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1.0 EXECUTIVE SUMMARY 

Laser and optical emission-based measurements have been developed and implemented for 
use on low-power hydrogen arcjet thrusters and xenon-propelled electric thrusters. In the 
case of low power hydrogen arcjets, these laser induce fluorescence measurements 
constitute the first complete set of data that characterize the velocity and temperature field of 
such a device. The research performed under the auspices of this NASA grant includes 
laser-based measurements of atomic hydrogen velocity and translational temperature, 
ultraviolet absorption measurements of ground state atomic hydrogen, Raman scattering 
measurements of the electronic ground state of molecular hydrogen, and optical emission 
based measurements of electronically excited atomic hydrogen, electron number density, 
and electron temperature. In addition, we have developed a collisional-radiative model of 
atomic hydrogen for use in conjunction with magnetohydrodynamic models to predict the 
plasma radiative spectrum, and near-electrode plasma models to better understand current 
transfer from the electrodes to the plasma. In the final year of the grant, a new program 
aimed at developing diagnostics for xenon plasma thrusters was initiated, and results on the 
use of diode lasers for interrogating Hall accelerator plasmas has been presented at recent 
conferences. 

The research is presented in detail in a total of 32 full length technical papers, 10 of which 
are in peer-reviewed journals. All of these reference the support provided by NASA 
Lewis, either through this grant, or, through donated equipment. 

2.0 PUBLICATIONS 

The following is a brief synopsis of the publications included in the Appendices of this 
report. The reader is also referred to the complete list of publications stemming from the 
research funded under this grant presented in Section 2.4. 

2.1 Near-Electrode Plasma Structure 

The need to better understand the structure and stability of a diffuse arc attachment at the 
anode of an arcjet thruster spawned a modelling effort on our part which focuses on the 
problem of current transfer between an electrode and a near-thermal plasma. This research 
culminated into two publications. In the first publication, we formulated a model for the 
case of a multi-fluid, single temperature, flowing thermal plasma impinging onto a plane 
anode, and identified two operating modes that were strongly influenced by the free-stream 
conditions. In the second paper (Meeks and Cappelli, 1993, see Appendix A), we 
extended the model to account for finite electron-ion recombination kinetics at the anode, 
and for the case of a multi-fluid, two temperature free stream condition. 

2.2 LIF Measurements of Temperature and Velocity 

We have used laser induced fluorescence to excite the B aimer-alpha transition in atomic 
hydrogen at 656 nm. By tuning the narrow bandwidth output of an argon-ion pumped 
ring-dye laser across the n = 3 -> n = 2 transition multiplet, we were able to resolve the 
Doppler-broadened and Doppler shifted lineshape to extract both translational temperature 
and velocity. The diagnostic strategy and first results were presented in Liebeskind et al., 
1993 (see Appendix B). In order to understand weather there is a significant slip between 
species comprising the flow in a hydrogen arcjet, LIF was performed on both atomic 
hydrogen and helium in a helium-seeded hydrogen arcjet. These results are presented in 
Liebeskind et al., 1995 (see also Appendix B). Finally, a comprehensive comparison of 
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the LEF data to MHD predictions of arcjet performance is presented in a paper co-authored 
by our laboratory, and engineers from Olin Aerospace Co., manufacturers of flight- 
qualified thrusters (Cappelli et al., 1995 in Appendix B). 

2.3 Raman Scattering Measurements of Molecular Hydrogen 

The Raman scattering diagnostic was motivated by the need to better understand 
performance limitations associated with frozen flow losses in arcjet thrusters. It provides 
the means of directly measuring the electronic ground state concentrations of molecular 
hydrogen. Measurements were performed for both cold-flow (see Boyd et al, Appendix 
C) and arc-heated flow conditions. In the cold-flow studies, experiments were in excellent 
agreement with the results of Monte Carlo simulations of Iain Boyd at Cornell University. 
Arc-heated flow studies were compared to both Monte Carlo and MHD predictions of the 
flow structure (see for example, AIAA 95-1956, Appendix C). The results support the 
conjecture that mass diffusion is an important transport mechanism which strongly 
influences the nozzle flow behavior. 

2.4 Complete List of Publications 


Peer-reviewed Publications 

1 . "The Non-Equilibrium Region of an Electrode in Contact with a Flowing Thermal 
Plasma," M.A. Cappelli, IEEE Trans. Plasma Sciences 21, pp.194-201, 1993. 

2 . "Two-Temperature Fluid Model for High Pressure Plasmas in Contact with Cooled 
Electrodes," E. Meeks and M.A. Cappelli, J. Appl. Phys. 73, pp. 3172-3182, 
1993. 

3. "LIF Diagnostic of Temperature and Velocity in a Hydrogen Arcjet," J.G. 
Liebeskind, R.K. Hanson, and M. A. Cappelli, Applied Optics 32, pp.6117- 
6127, 1993. 

4. "Numerical and Experimental Investigations of Low-Density Supersonic Jets of 
Hydrogen," I.D. Boyd, D.R. Beattie, and M.A. Cappelli, J. Fluid Mech. 280, pp. 
41-67, 1994. 

5. "Experimental Investigation of Velocity Slip Near an Arcjet Exit Plane," J.G. 
Liebeskind, R.K. Hanson, and M.A. Cappelli, AIAA Journal 33, pp. 373-375, 
1995. 

6. "Comparison of Hydrogen Arcjet Thruster Properties to Single-Fluid Model 
Predictions," M.A. Cappelli, J.G. Liebeskind, R.K. Hanson, G.W. Butler, and 
D.Q. King, submitted to Journal of Propulsion and Power, October, 1994. 

7 . "Interior Plasma Diagnostics of Arcjet Thrusters," M.A. Cappelli and P.V. Storm, 
Journal of Propulsion and Power, in press, 1995. 

8 . "Radiative Emission Analysis of an Expanding Hydrogen Arc Plasma I: Arc Region 
Diagnostics Through Axial Emission," P.V. Storm and M.A. Cappelli, submitted 
to JQSRT, 1995. 

9. "Radiative Emission Analysis of an Expanding Hydrogen Arc Plasma II: Plume 
Region Diagnostics Through Radial Emission," P.V. Storm and M.A. Cappelli, 
submitted to JQSRT, 1995. 
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10 . 


"Raman Scattering Measurements of Molecular Hydrogen in an Arc-Heated Flow, 
D.R. Beattie and M.A. Cappelli, submitted to Applied Optics. 


Conference Papers 

1 . "Velocity Measurements in a Hydrogen Arcjet using LIF," J.G. Liebeskind, R.K. 
Hanson, and M.A. Cappelli, AIAA 91-2112, 27th Joint Propulsion Conference, 
June 24-26, Sacramento, California, 1991. 

2. "Optical Diagnostics of a Low Power Hydrogen Arcjet," M.A. Cappelli, R.K. 
Hanson, J.G. Liebeskind, and D.H. Manzella, IEPC 91-091, 22nd International 
Electric Propulsion Conference, October 14-17, Viareggio, Italy, 1991. 

3. "Laser-Induced Fluorescence of Atomic Hydrogen in an Arcjet Thruster," J.G. 
Liebeskind, R.K. Hanson, and M.A. Cappelli, AIAA 92-0678, 30th Aerospace 
Sciences Meeting, January 6-9, Reno, NV 1992. 

4. "Vacuum Ultraviolet Absorption Measurements in a Hydrogen Arcjet," D.H. 
Manzella and M.A. Cappelli, AIAA 92-3564, 28th Joint Propulsion Conference, 
July 6-8, Nashville, TN, 1992. 

5 . "Modeling the Near-Electrode regions of Arcjets I: Coupling of the Flowfield to the 
Non-Equilibrium Boundary Layer," M.A. Cappelli, AIAA 92-3109, 28th Joint 
Propulsion Conference, July 6-8, Nashville, TN, 1992. 

6. "Flow Diagnostics of an Arcjet Using Laser-Induced Fluorescence,", J.G. 
Liebeskind, R.K. Hanson, and M.A. Cappelli, AIAA 92-3243, 28th Joint 
Propulsion Conference, July 6-8, Nashville, 'EM, 1992. 

7. "Molecular Hydrogen Raman Scattering in a Low Power Hydrogen Arcjet 
Thruster," D.R. Beattie and M.A. Cappelli, AIAA 92-3566, 28th Joint Propulsion 
Conference, July 6-8, Nashville, TN, 1992. 

8. Predictions," M.A. Cappelli, J.G. Liebeskind, R.K. Hanson, G.W. Butler, and 
D.Q. King, AIAA 93-0820, 31st Aerospace Sciences Meeting, January 11-14, 
Reno, NV, 1993. 

9 . "Monte Carlo and Experimental Studies of Nozzle Flow in a Low Power Hydrogen 
Arcjet," I.D. Boyd, M.A. Cappelli, and D.R. Beattie, AIAA 93-2529, 29th Joint 
Propulsion Conference, June 28-30, Monterey, CA, 1993. 

10. "Plume Characteristics of an Arcjet Thruster," J.G. Liebeskind, R.K. Hanson, and 
M.A. Cappelli, AIAA 93-2530, 29th Joint Propulsion Conference, June 28-30, 
Monterey, CA, 1993. 

11. "Axial Emission Diagnostics of a Low Power Hydrogen Arcjet Thruster," P.V. 
Storm and M.A. Cappelli, IEPC-93-219, Proceedings of the 23rd International 
Electric Propulsion Conference, Seattle, WA, September 13-16, 1993. 

12. "Chamber Effects on Plume Expansion for a Low- Power Hydrogen Arcjet," I.D. 
Boyd, D.R. Beattie, and M.A. Cappelli, IEPC-93-126, Proceedings of the 23rd 
International Electric Propulsion Conference, Seattle, WA, September 13-16, 1993. 

13. "A Direct Comparison of Hydrogen Arcjet Thruster Properties to Model 
Predictions," M.A. Cappelli, J.G. Liebeskind, R.K. Hanson, W.G. Butler, and 
D.Q. King, Proceedings of the IEPC-93-220, 23rd International Electric 
Propulsion Conference, Seattle, WA, September 13-16, 1993. 
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14. 


"LIF Measurements of Species Velocities in an Arcjet Plume," J.G. Liebeskind, 
R.K. Hanson, and M.A. Cappelli, IEPC-93-131, Proceedings of the 23rd 
International Electric Propulsion Conference, Seattle, WA, September 13-16, 1993. 

15. "Interior Plasma Diagnostics of Arcjet Thrusters," M.A. Cappelli and P.V. Storm , 
AIAA-94-2654, 25th AIAA Plasmadynamics and Lasers Conference, Colorado 
Springs, CO, June 20-23, 1994. 

16. "Semiconductor Laser Diagnostics Development for Xenon Plasmas," R.J. 
Cedolin, R.K. Hanson, and M.A. Cappelli, AIAA 94-2739, 30th Joint Propulsion 
Conference, Indianapolis, IN, June 27-29, 1994. 

17. "Axial Emission Measurements on a Medium Power Hydrogen Arcjet Thruster," 
P.V. Storm and M.A. Cappelli, AIAA 94-2743, 30th Joint Propulsion Conference, 
Indianapolis, IN, June 27-29, 1994. 

18. "Spectral Imaging of the Arcjet Electrode Region," D. Bems, P.V. Storm, and 
M.A. Cappelli, AIAA-95-1957, 26th AIAA Plasmadynamics and Lasers 
Conference, San Diego, CA, June 19-23, 1995. 

19. "High Spectral Resolution Emission Study of a Low Power Hydrogen Arcjet 
Plume," P.V. Storm and M.A. Cappelli, AIAA-95-1960, 26th AIAA 
Plasmadynamics and Lasers Conference, San Diego, CA, June 19-23, 1995. 

20. "Raman Scattering Measurements of Molecular Hydrogen in an Arcjet Thruster 
Plume," D.R. Beattie and M.A. Cappelli, AIAA-95-1956, 26th AIAA 
Plasmadynamics and Lasers Conference, San Diego, CA, June 19-23, 1995. 

21. "Laser-Induced Fluorescence Measurements within an Arcjet Thruster Nozzle," 
P.V. Storm and M.A. Cappelli, AIAA-95-2381, 31st AIAA Joint Propulsion 
Conference, San Diego, CA, July 10-12, 1995. 

22. "Pressure Measurements in the Plume of a Low power Arcjet Nozzle," W. Hargus 
Jr. and M.A. Cappelli, AIAA-95-2818, 31st AIAA Joint Propulsion Conference, 
San Diego, CA, July 10-12, 1995. 

3.0 GRADUATES STUDENTS FUNDED UNDER THIS GRANT 

John Liebeskind, "Temperature and Velocity Diagnostic for an Arcjet Plume Using Laser 
Induced Fluorescence." Ph.D. April, 1994. 

Douglas Beattie, "Raman Scattering Measurements of Molecular Hydrogen in an Arcjet 
Thruster Plume." Ph. D. April, 1995. 

Renato Cedolin, "Laser Diagnostics for Xenon Thrusters." Ph.D. expected 1996. 

David Manzella, "Laser and Emission Spectroscopy of Hall Accelerator Plumes." Ph.D. 
expected 1996. 
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Two-temperature fluid model for high-pressure plasmas in contact 
with cooled electrodes 

E Meeks PRECEDING PAGE BLANK NOT FILMED 

Sandia National Laboratories, Livermore, California 945551-0969, and Department of Mechanical 
Engineering, Stanford University, Stanford, California 94305 

M. A. Cappelli 

Department of Mechanical Engineering, Stanford University, Stanford, California 94305 

(Received 3 August 1992; accepted for publication 11 December 1992) 

A detailed continuum model of a weakly ionized plasma in stagnation flow against an 
electrically biased surface is presented. The model includes bulk fluid motion as well as electron 
and ion momentum and continuity equations in both general and quasineutral formulations. The 
inclusion of both electron and bulk energy equations allows for the determination of electron 
temperature and deviations from local thermodynamic equilibrium. Results for the quasineutral 
model demonstrate the importance of considering thermal nonequilibrium in the evaluation of 
plasma properties in the presheath region. The evaluation of transport properties is 
self-consistent with properties for the charged species depending on the calculated electric field 
and neutral number density. In addition to investigating thermal boundary layer and fluid flow 
effects, a finite electron-ion recombination rate at the electrode surface is employed and the 
consequences of varying this boundary condition are explored. 


I. INTRODUCTION 

The widespread use of plasma enhancement in materi- 
als processing has sparked new interest in the development 
of both glow-discharge and arc-jet plasma models. With 
few exceptions, most existing plasma models describe qui- 
escent systems with uniform or decoupled bulk (neutral) 
properties. In many plasma-processing reactors, however, 
ionized gases impinge against cooled substrates in a 
stagnation-point flow geometry. In such cases, there may 
be significant coupling between fluid-dynamic and thermal 
boundary-layer development, plasma electrical properties, 
and chemical kinetics. Unraveling these complex interac- 
tions is essential to understanding and predicting the ef- 
fects of plasma enhancement and electrical biasing of sub- 
strates. This article outlines the theoretical treatment of a 
weakly ionized, chemically reacting plasma in stagnation 
flow against a cooled, electrically biased substrate. The 
equations are then solved for a simplified case of a 
quasineutral argon plasma, to explore the effects of thermal 
and fluid boundary layers on plasma properties in the 
presheath region. 

Over the past 30 years there has been much work on 
the analysis and characterization of plasma systems with a 
variety of applications considered. Early studies were pri- 
marily concerned with the interaction of diagnostic probes 
and the plasma surrounding them (see, for example, Refs. 
1-5). Continuum probe models have usually focused on 
predicting current-voltage characteristics and have been re- 
viewed by Chung, Talbot, and Touryan. 6 Interest in arc 
modes and arc jet operation has also generated much work 
in plasma modeling near electrode surfaces. 7-10 These mod- 
els are directed toward understanding the principles behind 
diffuse and constricted arc modes and in predicting the 


potential fall in the vicinity of electrodes. Magnetohydro- 
dynamic (MHD) generators are often analyzed in a simi- 
lar manner with further attention to current-voltage char- 
acteristics and the effect of hydrodynamic boundary layers 
on plasma properties. 5,11,12 Numerous models have 
emerged recently to describe dc and rf glow discharges. 
These models attempt to capture the characteristics of a 
self-sustained plasma, predict voltage drops near the elec- 
trodes, and provide insight into various modes of discharge 
behavior. 11,13-18 While most models of plasma chemistry 
are limited to ionization and excitation reactions, 8,12,17,19,20 
Kushner 21 presented results for a quiescent, radio- 
frequency etching plasma that included an extensive chem- 
ical kinetic mechanism. We combine many of the physical 
details that have been introduced through these previous 
analyses to describe a comprehensive, consistent theory for 
a fluid continuum plasma in contact with a reactive elec- 
trode. In addition, we draw upon prior theoretical work in 
the area of neutral chemical-vapor-deposition modeling 
(see, for example, Refs. 22-28) which provides the 
groundwork for inclusion of multiple species in reacting, 
fluid systems. 

Among the many features common to the broad spec- 
trum of plasma analyses, it is generally recognized that the 
plasma region close to an electrode can be subdivided into 
several overlapping subregions in which different processes 
dominate. In particular, there usually exists 12 (a) an outer, 
quasineutral region in which net space charge is negligible 
but bulk properties such as velocity, temperature, number 
densities, and species fluxes may change significantly; (b) a 
sheath of length on the order of the plasma Debye length, 
where space charge dominates; and (c) a Knudsen sub- 
layer of thickness on the order of the electron or ion mean 
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free path, where continuum assumptions are no longer 
strictly valid. Several self-consistent models are reported 
that include the plasma sheath 15,16-18,19 * 20,29 " 32 with cou- 
pling of electric-field and space-charge effects. The 
quasineutral analysis presented here accounts for detailed 
nonequilibrium processes within the presheath region (a). 
This model can later be expanded to include solution of 
Poisson’s equation throughout the domain, coupling the 
space-charge-dominated sheath with presheath process in- 
teractions. 

One important drawback of previous plasma analyses 
is the almost exclusive use of a perfectly absorbing elec- 
trode boundary condition on the ions and electron number 
densities. The use of this boundary condition often intro- 
duces singularities causing numerical difficulties 6,12 and im- 
posing unrealistic constraints on the bulk plasma behavior. 
Even for a highly catalytic surface, a more appropriate 
condition would be a very high rate of recombination at the 
electrode, rather than a zero-density condition. Chung pro- 
vided some guidance in evaluating the effect of electron 
emission at an electrode, 33 but still maintained that the 
electrode is first perfectly absorbing. Graves 17 recently em- 
ployed a unity sticking coefficient, which provides a more 
appropriate formula for a surface with a high probability of 
recombination. In this work we include a finite reaction 
rate for electron-ion recombination and treat the rate co- 
efficient as a variable parameter. 


If. MODEL FORMULATION 
A. Governing equations 

The following governing equations provide a complete 
self-consistent formulation for a weakly ionized, contin- 
uum plasma in stagnation-point flow against an electrically 
biased substrate. Our numerical model includes solution of 
fully coupled electron and neutral energy equations. In 
addition we calculate net electron production rates by 
electron-impact ionization and electron third-body recom- 
bination with electron-temperature dependence. We derive 
the continuum equations by taking moments of the Boltz- 
mann equation and applying a stagnation-point flow simi- 
larity transformation to obtain the final one-dimensional 
set of continuity, momentum, and energy equations for 
both individual species and the bulk. The stagnation-flow 
geometry is illustrated in Fig. I. Further assumptions in- 
clude thermalization between positive ions and neutral spe- 
cies and the application of Chapman-Enskog theory 34 for 
determination of neutral transport properties. Transport 
properties of charged species derive from swarm data for 
drift velocities as functions of E/n, where E is the calcu- 
lated electric field, and n the neutral number density. The 
assumption of a quasineutral plasma then reduces the 
number of equations for the simplified calculation pre- 
sented in Sec. III. 

Equations (1)— (6) below describe the bulk fluid mo- 
tion with nonuniform temperature and pressure: 



FIG. 1. Schematic of stagnation-flow geometry. 
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These equations can be derived either from Navier-Stokes 
relations or from the summation of species conservation 
equations as derived from moments of the Boltzmann 
equation. In either case, a similarity transformation has 
been introduced for the stagnation-flow geometry. This 
transformation begins with the assumption of axisym- 
metry, where u and v are the mean-mass axial and radial 
velocities. Instead of the traditional stagnation flow formu- 
lation 35 in which a potential-flow solution determines the 
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pressure-gradient term in the radial momentum equation, 
we employ the finite-domain transformation described by 
Kee etal 36 This method introduces a stream function of 
the form ip(x,r) = U(x)r l , which exactly satisfies mass 
continuity. In neglecting edge effects, we assume that the 
radii of the substrate surface and flow inlet are infinitely 
large. Substituting the stream function into the axial and 
radial momentum equations results in a constant radial 
pressure gradient as indicated by Eq. (5). Also we note 
that u=2U/p and we can define V~{\/p)d V/dx = v/r, 
where p is the mixture mass density. In an ideal stagnation 
flow, temperature and concentration gradients occur only 
in the axial direction, precluding density variations in the 
radial direction. The dependent variables p, u, and V are 
therefore only functions of the axial coordinate x. For con- 
venience we choose to solve the continuity equation ( 1 ) 
along with the radial and axial momentum equations (2) 
and (3), rather than maintaining the stream-function no- 
tation. In the equations below, t] is the mixture viscosity. 
The average thermodynamic pressure P is assumed to be 
much larger than the radial contribution Ar*/2. The ideal 
gas law (4) relates pressure to density and temperature, 
with k B denoting the Boltzmann constant. Here n k is the 
number density and m k is the molecular mass of species k. 
The subscript e depicts electrons while t depicts electrop- 
ositive ions. There are a total of K, species including elec- 
trons and Ki positive ions. Although we are interested in 
the steady-state solution, we include transient terms in 
each equation to facilitate the description of time-stepping 
techniques in our numerical solution scheme. 

Equation (6) describes the bulk energy equation where 
neutrals and ions are in thermal equilibrium and their tem- 
perature is a function of the axial coordinate only. It is 
important to note that the bulk energy equation represents 
the summation of species energy equations over all species 
except electrons, which are allowed to depart from thermal 
equilibrium. The energy equation includes convection, dif- 
fusion, production, conduction, flow work, and Joule heat- 
ing terms. The last three terms describe energy gain due to 
elastic and inelastic collisions with electrons, followed by 
an electron production term. These three terms appear in 
the bulk energy equation due to the separate treatment of 
the electron energy and are described in more detail below. 
If the electrons were equilibrated at the bulk temperature 
T, and the electron energy equation added to Eq. (6), 
these terms would no longer appear in the bulk energy 
equation. In the ion Joule-heating term, the product Zp is 
the charge of the x'th electropositive ion, where — e denotes 
the charge of an electron. Also, h k is the specific enthalpy 
of species k, c pk is the species’ constant-pressure specific 
heat, and h k is the net species’ production rate. c p repre- 
sents the mixture-averaged specific heat, while A is the bulk 
thermal conductivity. 

Equations (7) and (8) describe continuity and conser- 
vation of momentum, respectively, for individual species 
within the plasma bulk: 

dn k d 
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Again, the similarity transformation results appear in the 
bulk convection terms. We have introduced the variable P k 
to describe the diffusive flux of species k, where r k =n k U k 
and U k is diffusion velocity. Equation (9) represents Pois- 
son’s equation, relating the electric field gradient to the 
local net space charge. Here e 0 is the free-space permittiv- 
ity. We can relate the total current density J to the ion and 
electron fluxes, as follows: 



z,r,-r. 



Zfii-ri' 


For neutral species the general form of the diffusive 
fluxes result from a multicomponent formulation including 
thermal diffusion as shown in Eq. (8a). 34 Here dj is the 
concentration gradient with respect to x. D kj is the ordi- 
nary diffusion coefficient while D k is the thermal diffusion 
coefficient. The mean molecular mass of the mixture is 
given by m. We note that Eq. (8a) is an explicit relation 
that could be substituted directly into Eq. (7) for the neu- 
tral species. However, we choose to keep T* as a separate 
dependent variable for consistency with the treatment of 
ions and electrons. For ions and electrons, the full species 
momentum equations must be solved as indicated by Eq. 
(8b) for ions and Eq. (8c) for electrons. The ion and 
electron momentum equations include balance between 
forces due to pressure, inertia, and electric field as well as 
momentum loss due to collisions with other species. Vis- 
cous forces on electrons and ions are neglected. For the 
ions, we assume that charge-neutral collisions dominate 
momentum transfer, consistent with the assumption of a 
weakly ionized plasma. The parameter v w represents the 
momentum-averaged momentum-transfer collision fre- 
quency summed over all neutral species. For electrons, 
we consider all heavy species (i.e., all species besides elec- 
trons) and introduce the parameter v eH = 1 h v eh to repre- 


•V * 100*5 


E. Meeks and M. A. CappelU 3174 


(12) 


sent the momentum-averaged momentum-transfer colli- 
sion frequency summed over all heavy species. The 
collision frequencies are derived from swarm data for drift 
velocities with electric-field dependence, as further dis- 
cussed below. 

Equation (10) is the electron energy equation in which 
\k B T e represents the average energy of an electron: 


3 dT e dn e 3 dT e dn e 

2 nJiB sr- k ‘ r - sr +2 “•**“ & 




A m h 


+|wr,-n« o. 


(10) 


We neglect all thermal diffusion effects on the electron 
energy. The steady-state equation includes convective 
terms, diffusion, conduction. Joule heating, energy loss due 
to both elastic and inelastic collisions, and a production 
term. From left- to right-hand sides, the first two terms in 
Eq. ( 10) are transient, followed by two terms which rep- 
resent a combination of convection and flow work analo- 
gous to the second and sixth terms in the bulk energy 
equation (6). The next term represents diffusion of elec- 
tron enthalpy analogous to the third term in Eq. (6), fol- 
lowed by a conduction term. Next is the electron Joule- 
heating contribution representing the increase in electron 
energy due to work by the electric field. Coupling the elec- 
tron energy equation to the bulk energy equation, the next 
term represents the transfer of energy between electrons 
and heavy species through elastic collisions. Inelastic col- 
lision losses result from energy lost by electrons participat- 
ing in electron third-body reactions, and take the form 
R imi =1/jAH n where q r is the rate of progress of the rth 
electron third-body reaction, and A H r is the heat of reac- 
tion. A positive sign for AH r indicates that the reaction is 
endothermic in its forward direction. Finally, the produc- 
tion term can be interpreted as the energy loss to the elec- 
tron swarm required to thermalize electrons produced 
through ionization of neutral species. The electron equa- 
tion is included here in its most general form, although 
some terms may prove to be insignificant under the condi- 
tions considered below. 


B. Quasineutral assumptions 

We now modify the set of governing equations for the 
case of a quasineutral plasma, 6 neglecting space-charge ef- 
fects in the thin-sheath limit: 

k, 

n e — X (ID 

i=i 


J ' 

r,+-- X z, r,=o, 

e i=i 

J ( du — \ 

-[2^2V +v ^ 

5i ( k B d e — \ 

- X 2, — — (n;T e ) +Z i — n i E-Z i T i v eH =0. 

,=i \ m e dx m e ) 

(13) 

This approach allows us to investigate thermal nonequilib- 
rium in the presheath region of a flowing plasma. Solution 
of the full set of equations described above can later be 
used to evaluate additional interactions within the sheath. 
The fundamental assumption of the quasineutral formula- 
tion is that, throughout the region of interest, the local net 
charge is approximately zero. To this end, the electron 
continuity equation is replaced by the charge balance con- 
dition given by Eq. (11). In addition we assume that the 
production rate of electrons is always balanced by the pro- 
duction of ions with n,=2,Z,nj, and that the gradient of 
the electron number density is balanced by gradients of 
ions, such that dn/dx ~ 2, Z,( dn/dx ) . These assumptions 
result in a constant total current density J. We therefore 
replace the electron momentum equation (8c) with the 
consequent quasineutral relation (12). Finally, in order to 
solve for the ambipolar electric field, we subtract Eq. (8c) 
from the sum of the ion momentum Eq. (8b) weighted by 
the integer charge Z; of each ion. The result is Eq. (13) 
which replaces Poisson’s equation (9). We now have ten 
equations with ten unknowns: u(x), V(x), P(x), T(x), 
p(x). A, n k (x), r*(x), E(x), and T e (x), with T* and n k 
representing unknowns for each species (neutral, ion or 
electron). 

C. Boundary conditions 

Solution of the above set of equations requires bound- 
ary conditions both at the reactor inlet and at the electrode 
surface. At the inlet (x=L) we specify the freestream con- 
ditions with u(L)=u ai , V(L)= 0, P{L) =P a> , T(L) 
= r„, n k (L)=n ka> , and T t =T ew . For the results pre- 
sented here, we additionally let T ew = T w and set the inlet 
species concentrations to their equilibrium values at that 
temperature, although this is not required for solution of 
the equations. 

At the electrode, the conditions are more complex. The 
radial velocity obeys the no-slip condition at the surface, 
with F(0) =0. The inclusion of heterogeneous reactions 
that may result in a net deposit on the wall as well as a net 
current into or out of the wall results in a nonzero value for 
w(0). The bulk velocity condition at the electrode becomes 



where s k is the net production rate of species k through 
heterogeneous reaction. We assume that the neutrals and 
ions thermally equilibrate with the electrode such that 
T(0) = T w> and T w is held constant The electrons, how- 
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TABLE I. Electron and ion transport properties. 


Symbol 

Transport 

property 

Formula 

v rt 

electron momentum- 
transfer 

collision frequency 
with species h 

v M =e£/m,ut rt 

A* 

binary electron 
diffusion 
coefficient in 
heavy species h 

D, h =k B T/m,v th 


electron thermal 
conductivity 

lWr - 


v ,r 

ion momentum- 
transfer 

collision frequency 
with species n 

— eE 

V W __ ... • — l m 

T7ij n Wj n /n,“T“ 

v iN 

ion momentum- 
averaged 

collision frequency 


D/„ 

binary ion diffusion 

coefficient 

in neutral species n 

D„ = k a T/m,v,' 


ion thermal 
conductivity 

/9y,-5\ , n, c„ 

Vr,-lJ c* 

Vi 

ion viscosity 

V, = 4k b T n/ TtVjf/ 


ever, may be sufficiently energetic to remain out of equi- 
librium at the electrode surface. The boundary condition 
on the electron energy is therefore a balance between en- 
thalpy loss from net current flow into the surface, net de- 
struction of electrons through surface recombination, con- 
duction due to electron temperature gradient, and flow 
work. This relation is stated as 



-A, 


dj_ , 

dx 


-(unJceT') 


jr=0 


=0 = 0 . 


U=0 


The boundary condition on species fluxes must also 
balance production rates of species through surface reac- 
tion. With a net current / specified in the electrically biased 
case, the electron, ion, and neutral flux conditions are as 
follows: 


r e (0)=s e -(J/e)-n e u\ x=0 , 
r,( 0 ) —Sj — fijU | x _ o, i— lj— , 


r*(0)=i*— n*u| x=0 , k=£e or /. 

K- 

Here we also require that s e = 2 f= l and assume that 
only electrons carry current into or out of the electrode. 
This completes the boundary-condition requirements for a 
well-posed set of governing equations. 


D. Numerical solution 

Through finite-difference approximations, we dis- 
cretize the set of differential equations over a one- 
dimensional computational grid which spans the distance 
between the electrode and the inlet or frees tream at x— L. 
Solution of the resulting algebraic equations employs a 
damped, modified Newton algorithm designed for 
boundary-value problems. 37 This numerical scheme, as first 
applied to premixed flame modeling, 38 has proven to be a 
successful method for convergence of stiff, nonlinear sys- 
tems of equations. Solution begins on a coarse mesh with 
trial guesses for the dependent variables. The coarse-mesh 
solutions then provide the initial guess for adaptively re- 
fined grids, until the solution achieves the desired accu- 
racy. In the event that the Newton algorithm fails to con- 
verge, a series of time integrations proceed from the trial 
estimate to gradually approach the steady-state solution. 
Since a successful Newton iteration is much faster than the 
time-integration procedure, time stepping only conditions 
the trial estimate until it is sufficiently good for the Newton 
iteration to converge. The implementation of our numeri- 
cal scheme and the overall program structure derives from 
the neutral-gas stagnation-flow software documented as 
the Sandia computer code SPIN. 26 


E. Incorporation of transport properties and 
chemical production terms 

In our general formulation for a reacting multispecies 
flow, we require transport properties for neutral species 
including thermal conductivity, diffusion coefficients, and 
mixture viscosity. Calculation of these properties employs 
a modified version of the Sandia transport package. 22 For 
interactions between neutral species, properties derive from 
standard kinetic theory where collision integrals are based 
on Stockmayer potentials using a Maxwellian neutral tem- 
perature and Leonard-Jones potential well depths. 39 To 
determine collision frequencies for charged species, this 
method is no longer valid and we must include electric- 
field dependence. In an effort designed to be parallel and 
analogous to the Sandia transport approach, we include the 
generalized capability of incorporating electron and ion 
transport dependence on electric field in a coupled manner. 
This is achieved by first mathematically fitting measured or 
calculated swarm parameters in the form of drift velocity, 
w ek vs E/n k for each charge-neutral species pair. We then 
access this information during numerical iteration for local 
values of E/n h The drift velocity w ek determined from 
swarm data is the velocity at which electrons (or ions) 
diffuse through a uniform medium in the presence of an 
applied electric field. This velocity is directly related to the 
momentum-averaged collision frequency between electrons 
and the background gas for a particular value of E/n h as 
indicated in Table I. The total momentum-average colli- 
sion frequency for electrons is determined by summing 
over all heavy species, while ion collision frequencies are 
summed over all neutral species. Charged species’ trans- 
port properties are then determined from or as 
shown in Table I. 40 
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Similarly, for inclusion of neutral and charged-species 
chemical production terms, we employ a modified version 
of the Sandia subroutine package CHEMKIN . 24,25 Here, as- 
sumption of Maxwellian distribution functions for neutral 
species allows the use of temperature-dependent rate for- 
mulas for chemical reactions that involve neutral species 
only. In general, electron reactions such as electron-impact 
ionization or electron-impact dissociation could be treated 
in a similar manner to electron and ion transport described 
above. In this case, ionization, dissociation, or recombina- 
tion rate coefficients as functions of E/n could be em- 
ployed. Alternatively, reaction cross sections as functions 
of electron energy could be integrated over an electron 
distribution function determined by a local field approxi- 
mation. For the high-pressure conditions considered here, 
however, we assume that the isotropic part of the local 
electron distribution is close to a Maxwellian function of 
In this way, rates for reactions that involve electron 
impact depend only on the local electron temperature and 
species concentrations. Also, the assumption of a near- 
Maxwell ian electron temperature simplifies the use of de- 
tailed balancing in determining reverse reaction rates. Fi- 
nally, use of the Sandia Surface CHEMKIN package 
facilitates incorporation of finite-rate chemistry at the elec- 
trode. 28 


III. RESULTS AND DISCUSSION 

For purposes of comparison, we begin by analyzing 
plasma conditions typical of arc-jet operation with stagna- 
tion flow against a perfectly catalytic electrode. We choose 
plasma free-stream conditions to be argon gas in local ther- 
modynamic equilibrium (LTE) at 6000 K and 1 atm. In 
addition, we specify a free-stream velocity of u co = 570 
cm/s at a distance of 5 cm above the electrode. Approxi- 
mating the perfectly catalytic condition requires a very 
high rate for electron-ion recombination, since we maintain 
a finite-rate formulation. To achieve this condition, we 
choose a recombination rate such that the electron gener- 
ation rate, mole fraction, and flux behavior compare well 
with analyses that employ zero electron density at the elec- 
trode. To compare with analysis by Cappelli, 9 we employ a 
reaction rate coefficient of approximately 10 15 cm 4 /s mol 
for the recombination reaction 

Ar + -f-e=>Ar 

occurring at the electrode surface. The effects of reducing 
this reaction rate to a more realistic value are discussed 
later in this section. In the gas phase, we allow both ion- 
ization and recombination through electron third-body im- 
pact. Thus, for the reaction stated as 

Ar + +e+e^>Ar-fe, 

we employ a forward reaction rate from Hinnov and Hir- 
shberg, 41 with 

*/= 3.959 X10 39 77 9/2 cmVsmol 2 , 

and the reverse rate derived from detailed balancing be- 
comes 



FIG. 2. Spatially varying electric field for various current values and 
uniform neutral temperature case. 
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k r = 1.906 X l0 i2 T e 3 exp/ — J cmVs mol. 


For the results presented here, we consider no other chem- 
ical reactions in either the gas or at the surface. 

We consider the three cases of electrode operation: un- 
biased or floating (/=0), anode (/> 0) and cathode (J 
< 0) . It is convenient to introduce a reference value for the 
electrode current density for scaling purposes. We there- 
fore let JR=n eiB ((i e +Hj)(.k B T ea> /l R ), where l R is the ap- 
proximate recombination length at the inlet conditions, 
while n e and /x, are characteristic values of the electron and 
ion mobilities at the inlet. 9 For the conditions described 
above, we use fi f =0.6 mVv s, /x,— 0.002 m 2 /V s, and ar- 
rive at a value of approximately 6.S mm for l R and 230 
A/m 2 for J R . Also, since we use the electrode current as a 
boundary condition, we implicitly assume that, in the case 
of J<0, the cathode is capable of providing the specified 
current. Here we are only interested in the effect of the net 
current on the plasma behavior. If we instead specified a 
fixed material for our electrode, we would have to restrict 
cathode currents to those consistent with the material 
work function at the specified electrode temperature. 


A. Uniform neutral temperature case 


If we first neglect the thermal boundary layer and set 
T u /T a> = 1.0, the results of our stagnation flow analysis 
are as shown in Figs. 2 and 3. These results underline the 
importance of including a separate electron energy equa- 
tion when a net current passes through the electrode, es- 
pecially in the cathode case. In all cases, a sharp electrical 
boundary layer develops due to ambipolar diffusion, as il- 
lustrated in Fig. 2, even in the unbiased case. The magni- 
tude of this presheath electric field is largely determined by 
the electron-ion recombination boundary condition, as fur- 
ther discussed below. Figure 3 shows the importance of 
electron temperature in determining electron and ion be- 
havior. In the cathode case (/<0), electrons diffuse with 
the electric field and are thereby energized through Joule 
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FIG. 3. Effect of varying current on plasma properties for uniform neutral temperature case. 


heating. In the anode case, however, electrons diffuse 
against the field and thereby lose energy through diffusion 
cooling. In either case, the presence of the electrode causes 
the electron temperature to depart significantly from equi- 
librium with the neutral/ion temperature, even when ne- 
glecting the neutral thermal boundary layer. 

The importance of thermal nonequilibrium is apparent 
in Figs. 3(b), 3(c), and 3(d). For the cathode case, the 
elevated electron temperature causes the plasma to ionize 
further and results in a source of electrons and ions near 
the surface. This source causes the ions to reverse direction 
midstream such that ions diffuse away from the generation 
location, despite the negative bias of the electrode. In the 
case of the anode, recombination occurs in the diffusion- 
cooling region, an effect that in turn reduces the electric- 
field value. For the anode, convection of the inlet gas mix- 
ture dominates electron dynamics, by limiting the effect of 
electron-ion recombination. We also note that, as the elec- 
trode is biased increasingly positive (i.e., behaves more and 
more like an anode), the ion flux toward the surface grad- 
ually approaches zero. As previously discussed by Self and 
Eskin 12 and Cappelli, 9 the point at which the ion flux 
reaches zero can be considered an ion-current saturation 
limit, since further increases in the drawn current would 
require ion emission from the electrode. 


B. Inclusion of the neutral thermal boundary layer 

Survival of a realistic electrode in the above conditions 
requires active cooling to a temperature well below the 
free-stream plasma temperature. The presence of a cooled 
wall results in the development of a thermal boundary 
layer for the bulk gas, as shown in Fig. 4(a). Here we have 
fixed the electrode temperature at T w = 1200 K or T JT m 
=0.2. If there were no net flow at x= L, the bulk temper- 
ature profile would be linear. The magnitude of the inlet 
velocity and the bulk plasma transport properties deter- 
mine the thermal boundary layer thickness. Although the 
neutral/ion temperature includes Joule heating of the ions 
and collisional exchange with electrons, these effects are 
negligible for the currents considered in Fig. 4. For the 
bulk temperature, convection and conduction terms dom- 
inate. 

Comparison of Fig. 4(b) with Fig. 2 above demon- 
strates a significant influence of electrode cooling on the 
electric field. This effect is more easily understood by con- 
sidering the results presented in Fig. 5. In Fig. 5(a) the 
presence of the bulk-flow thermal boundary layer sup- 
presses the effect of Joule heating in the cathode case 
through increased collisional energy transfer between the 
electrons and the bulk. In the floating case, collision energy 
exchange causes the electrons to cool below their free- 
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FIG. 4. (a) Neutral temperature profile for cooled electrode and (b) 
resulting electric field for cathode, anode, and floating cases. 


stream temperature. In the anode case, however, the 
boundary layer has a smaller effect on electron tempera- 
ture, since energy exchange with the bulk is now smaller 
than in Fig. 3(a). In all cases the electron temperature 
remains significantly out of equilibrium with the neutral 
species near the electrode. In the cathode and floating case, 
lower electron temperatures result in decreased ionization, 
increased recombination, and a shift of the electron pro- 
duction peak nearer to the electrode surface. The conse- 
quent reduction in ion flux away from the production 
source suppresses the tendency of the electric field to 
change sign midstream. In all cases, electrode cooling re- 
sults in significantly larger magnitudes of the electric field 
near the electrode. With decreased ionization, increased 
recombination in the bulk and a high recombination rate at 
the surface, the electric field compensates for the tendency 
of electrons to diffuse rapidly toward the surface, maintain- 
ing overall charge neutrality. 

An important issue in arc-jet modeling is the determi- 
nation of the potential fall near an electrode surface. Al- 
though we cannot provide a quantitative measure of the 
total potential drop through the plasma sheath without 


inclusion of Poisson’s equation, we can determine the qual- 
itative behavior of the potential fall in the presheath region 
that results from the presence of the surface. We define 
A<p quas as the integral of the electric field, excluding the 
field contribution arising from the net current flow, for the 
quasineutral formulation: 

A^quas = <P U=0 ~ <P I x= Z. 

= [ L [E{x)-E{L)]dx 
Jo 

= ^ E(x)dx-E(L)L. 

The results for A<p quis as a function of total current density 
are shown in Fig. 6. When the thermal boundary layer is 
present, a distinct transition in potential behavior occurs 
for both the anode and cathode cases. For anode current 
values less than J/Jr ~ 8, the plasma has positive resis- 
tance, corresponding to a diffuse mode. For higher cur- 
rents, however, the plasma behavior reverses, correspond- 
ing to a possible transition to a constricted or arc mode. 
For the cathode case, this transition occurs at much lower 
currents. Neglect of the neutral thermal boundary layer 
would result in no observed diffuse mode for the cathode 
and a transition at higher currents for the anode case. 

The ion fluxes at x=0 for the anode cases including 
the thermal boundary layer are shown in Fig. 7. The ion 
flux towards the electrode exhibits a minimum correspond- 
ing to the transition shown in Fig. 6. 


C. Varying the rate of surface recombination 

In the above cases we observe that the near-catalytic 
boundary condition imposed at the electrode surface 
largely determines the magnitude of the electric field. Fig- 
ure 8 illustrates the effect of reducing the surface recombi- 
nation rate on the plasma behavior for the case of a floating 
electrode. An important result indicated in Fig. 8(a) is 
that the electron number density at the electrode ap- 
proaches a finite nonzero value as the surface recombina- 
tion rate increases from a low value to the near-catalytic 
value. This would suggest that the extreme condition of 
n e = 0 employed by many authors (see, for example. Refs. 
9, 12, 14, and 19) is unrealistic. Also, the fact that the 
electron number density appears to asymptote at high sur- 
face recombination values indicates that our choice of 10 15 
cm Vs mol was appropriate in approximating a perfectly 
catalytic surface. In Fig. 8(b) the magnitude of the electric 
field at the surface becomes very large as the rate of catal- 
ysis increases. This is consistent with the electric field’s 
inverse relationship with the electron number density at 
the surface and has been the source of singularities when a 
zero number-density condition is employed. 12 Figure 8(c) 
illustrates the effect of the surface electron-loss rate to the 
local gas-phase production rate at the electrode. For low 
surface recombination, electrons recombine in the gas, con- 
sistent with the boundary-layer cooling of the electrons. At 
high surface-recombination rates, however, the loss of elec- 
trons drives the ionization/recombination balance toward 
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FIG. 5. Effect of varying current on plasma properties including thermal boundary layer at electrode. 


ionization. Although a lower surface recombination rate 
reduces the electron temperature near the surface, this ef- 
fect is small with the value of T JT x at *=0 remaining 
between 0.8 and 0.9 for the electrically floating case con- 
sidered here. 

The above observations indicate the strong dependence 
of plasma behavior on the choice of a surface recombina- 
tion rate. It is therefore difficult to predict quantitatively 
important plasma characteristics in this quasineutral for- 


mulation without first obtaining a better understanding of 
the surface third-body behavior in electron-ion recombina- 
tion. However, inclusion of the plasma sheath in the com- 
plete formulation may significantly reduce the consequence 
of the rate choice. In the presence of net space charge near 
the surface, electrons are repelled from the sheath bound- 
ary. In this way, the electron loss may become transport 
limited, rather than surface-kinetic limited. 



FIG. 6. Potential drop between the plasma free stream and the electrode 
as a function of total current density. 
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FIG. 7. Ion flux at x=0 as a function of total current density for the 
anode case. 
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FIG. 8. Effects of varying the surface recombination rate coefficient on 
(a) electron number density, (b) electric field, and (c) electron genera- 
tion rate near the surface for the case of 7=0. 

IV. CONCLUSIONS 

We have presented a general formulation for a chem- 
ically reacting, multicomponent, weakly ionized plasma in 
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stagnation flow against a cooled electrode. Through inves- 
tigation of the quasineutral or thin-sheath limiting case, we 
provide a broad understanding of various parameters on 
plasma behavior in the presheath region. In particular, we 
have demonstrated the importance of including the neutral 
thermal boundary layer that arises from electrode cooling 
and electron thermal nonequilibrium near the stagnation 
surface. These effects have important consequences on the 
determination of the presheath potential fall and prediction 
of a transition between diffuse and constricted modes of arc 
attachment. In addition, we have demonstrated the high 
sensitivity of the magnitude of the ambipolar electric field 
on the choice of boundary condition for ions and electrons 
near the electrode. In future analyses, inclusion of Pois- 
son’s equation in the general formulation can be compared 
with the thin-sheath approximation to better understand 
the consequences of the quasineiitrai assumptions. 
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Laser-induced fluorescence diagnostic 
for temperature and velocity measurements 
in a hydrogen arcjet plume 


John G. Liebeskind, Ronald K. Hanson, and Mark A. Cappelli 


A diagnostic has been developed to measure velocity and tr anslational temperature in the plume of a 
1-kW-class arcjet thruster operating on hydrogen. Laser-induced fluor escen ce with a narrow-band cw 
laser is used to probe the Balmer a transition of excited atnmir hydrogen. The velocity is determined 
from the Doppler shift of the fluorescence excitation spectrum, whereas the temperature is inferred from 
the lineshape. Analysis shows that although Doppler broadening is only significant broadening 
mec h a n is m , the fine structure of the transition must be into account. Near the exit plane, axial 
velocities vary from 4 to 14 km /s, radial velocities vary from 0 to 4 km /s, and swirl velocities are shown to 
be relatively s ma ll. Temperatures from 1000 to 5000 K indicate high diaanriatinn fractions. 

Key words: Laser-induced fluorescence, arcjet, atomic hydrogen, Doppler shift, Doppler broadening, 
velocity, temperature. 


Introduction 

Arcjets are effective devices for producing high- 
enthalpy gas flows for applications in materials pro- 
cessing, metal cutting, atmospheric reentry simula- 
tion, and spacecraft propulsion. This paper is 
motivated by the propulsion application, in which 
arcjet thrusters are targeted for satellite station 
keeping, vehicle orbit transfer, and space exploration. 1 
Arcjets’ value in this application stems from the high 
mean mass velocity that is generated at the nozzle 
exit. This figure of merit is referred to as the specific 
impulse, or I, p , and is determined from the ratio of 
thrust to propellant mass-flow rate. The 1^ of chemi- 
cal rockets is limited because the stagnation tempera- 
ture cannot exceed the adiabatic flame temperature. 
This limits the amount of energy available for conver- 
sion to kinetic energy. In arcjet thrusters, the en- 
ergy is coupled via a constricted arc, where much 
higher temperatures can be obtained, and in practice 
an I sp can be achieved that is 2 to 3 times that of the 
best chemical rockets. 

The interaction of the arc with the gasdynamic flow 
creates a complex environment that is difficult to 
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analyze. The cathode is located in a stagnation 
plenum in the converging section of the nozzle (see 
Fig. 1). The arc extends from the tip of the cathode 
through the throat and attaches diffusely to the 
nozzle-anode in the low-pressure diverging section. 
The submillimeter throat carries high currents (10 A) 
and sustains high radial temperature and velocity 
gradients. Simplified analytical models suggest cen- 
terline temperatures in excess of 30,000 K, 2 whereas 
the throat walls are below the tungsten nozzle’s 
melting point. The ionization and molecular disso- 
ciation followed by the large expansion and the high 
velocities contribute to nonequilibrium conditions at 
the exit plane. Diagnostics of the plume can lead to a 
better understanding of the plasma dynamics by 
determining the conditions of the flow field. Veloc- 
ity is important because it correlates with thrust, 
whereas kinetic temperature indicates the amount of 
enthalpy remaining in the gas. 

This paper describes the development of a diagnos- 
tic that uses laser-induced fluorescence (LIF) to mea- 
sure velocity and temperature in the plume of an 
arcjet thruster. The laser probe is focused to a small 
diameter in the arcjet plume. Subsequent fluores- 
cence is collected and imaged onto a detector. Flow 
properties are determined from the behavior of the 
fluorescence signal as a function of the laser fre- 
quency (wavelength). This nonintrusive spectro- 
scopic technique has several advantages. Measure- 
ments can be made in the vicinity of the exit plane 
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Fig. 1. Arcjet nozzle configuration. The throat diameter is 0.64 
mm, and the expansion area ratio is 225. 


where solid probes would fail because of the high 
enthalpy of the flow or would otherwise be difficult to 
interpret. In addition, the ability to focus the probe 
laser to small sizes allows for high spatial resolution. 
This is particularly important in arcjets, where large 
variations in velocity and temperature occur across 
the small, 1-cm diameter of the nozzle exit. 

Velocity is determined from the Doppler shift of the 
fluorescence excitation spectrum. This technique 
has been demonstrated in a variety of flows and has 
been adapted for use in arcjets. 3 - 4 - 5 The transla- 
tional temperature is determined from the shape of 
the fluorescence excitation spectrum. In contrast to 
velocity, the temperature is more difficult to infer 
because of competing mechanisms that act to broaden 
the measured profile and the complex nature of the 
overlapping line components. This paper uses a 
tunable narrow-band cw laser to resolve the line- 
absorption profile generated through fluorescence 
excitation. 

Theory 

LIF involves excitation of an atomic or molecular 
species from a lower to a higher quantum state by 
absorption of laser radiation, with subsequent emis- 
sion of radiation (fluorescence). To avoid use of 
additive tracer molecules, an inherent species with a 
suitable transition is preferred. The exhaust plume 
of a hydrogen-fueled arcjet consists primarily of 
neutral and singly ionized atomic and molecular 
hydrogen as well as electrons. The Balmer a transi- 
tion of atomic hydrogen was selected owing to its 
accessibility at 656.28 nm, the substantial population 
of atomic hydrogen in the absorbing state at arcrjet 
conditions, and the transition’s relatively large oscil- 
lator strength. A partial energy level diagram for 
hydrogen is shown in Fig. 2 (Ref. 6). The Balmer a 
transition gives rise to the first spectral line of the 
well-known Balmer series. In this series the lower 
energy level is the first exdted state with n (the 
principal quantum number) = 2. Note that the 
Lyman series, in which the lower level is the ground 
state (n = 1), involves transitions in the vacuum 
ultraviolet for which suitable cw laser sources are not 
yet available. Fluorescence is obtained when the 



Maine Frequency (cm* ') 

Fig. 2. Hydrogen atom energy level diagram with the fine struc- 
ture that is due to spin-orbit coupling. The corresponding transi- 
tions are shown below the diagram with their relative separations 
and intensities. 


laser is tuned to the frequency resonant with the 
energy difference between the two states of an al- 
lowed transition. The laser, which is focused to the 
point of interest in the flow, induces some of the 
atoms in the lower state to absorb a photon, thus 
pumping them into the upper state. Some of these 
atoms in the upper state will spontaneously decay 
back to the lower state by emitting a photon. If the 
atom undergoes an elastic collision in the interim 
between the absorption and re-emission processes, 
the photon will be emitted in an arbitrary direction 
with a random phase. The rate at which photons 
arrive at the detector is 7 

n 

n p = nJiJVt — > (1) 

where n u is the number density of the upper level, A ui 
is the Einstein spontaneous emission coefficient, V is 
the fluorescence probe volume, e is the efficiency of 
the collection optics, and ft is the solid angle sub- 
tended by the collection optics. Because of the high 
temperature, the population in the upper level, n u , 
may be significant in the absence of laser excitation. 
This can lead to a large background signal. For this 
paper, the number of atoms excited to n = 3 by the 
probe laser is 1 to 2 orders of magnitude smaller than 
the number of atoms in the n = 3 state without laser 
excitation. 

In considering LIF as a velocity diagnostic, we first 
recognize that each species comprising the flow will 
have a distribution of velocities. Let u represent the 
ith species’ velocity component along a specified axis 
and fi(u) be the distribution function of this species’ 
velocity component. The distribution function rep- 
resents the fraction of species i having a velocity 
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The mean velocity for a species i along the specified 
coordinate can now be defined by 

/* -*■ 30 

Vi = j ufi(u)du. (3) 

Mean velocities for different species may vary because 
of diffusion. The mean mass velocity of the flow is 


2 p iVi 


(4) 


where p* is the local density of species i. LIF tech- 
niques determine the mean velocity of the probed 
species ( ) along the direction of laser propagation 
from the Doppler shift of the fluorescence excitation 
spectrum. The Doppler shift, Av, is given by 


u i 

Av = v 0 - > (5) 

where v 0 is the laser frequency, i/; is the mean velocity 
in the direction of laser propagation, and c is the 
speed of light. A stationary hydrogen p lasma dis- 
charge is used to provide a reference (unshifted) 
frequency. 

The temperature is determined from the shape of 
the fluorescence excitation spectrum. However, be- 
cause the Balmer a transition actually consists of 
several overlapping transitions, any effort to analyze 
line shapes must account for this structure. The 
upper and lower states are split by spin-orbit coupling. 
The split states, referred to as the fine structure, are 
identified by the quantum number j, which indicates 
the total angular momentum. Selection rules allow 
for five different optical transitions between the two 
manifolds. The line shape function of the B alm er a 
transition can be modeled as the superposition of 
these five components with their appropriate broaden- 
ing mechanisms and relative intensities. The fine 
structure and the five transition components are 
shown in Fig. 2. The relative position and intensi- 
ties of each component are listed in Table l. 8 Note 
that the last two components each consist of two 
degenerate transitions. 

Given the relative position and intensities of the 
various components, the line-shape function of the 
Balmer a transition is determined by the broadening 
mechanisms. The predominant broadening mecha- 
nism associated with the Balmer a transition in 
low-power (1-5 kW) arcjet exhausts is Doppler broad- 
ening. Each velocity class u, where u is the velocity 
component in the direction of laser propagation, will 
have its own Doppler shift. Generally the distribu- 
tion function f{u) is taken to be Maxwellian centered 


about the mean velocity u,-. Although large gradi- 
ents in the temperature and velocity can distort the 
Maxwellian distribution, 9 we consider these distor- 
tions to be ne gligib le. The Doppler shifts are com- 
bined with the Maxwellian distribution to give a 
Gaussian fluorescence intensity profile 10 


^ 


/In 2\ 1/2 
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$ d(v ) is generally referred to as the normalized 
line-shape function and Av D is the FWHM (full width 
at half maximum or Doppler half-width) given by 


(8kT In 2\ 1/2 

’ v °'r^n ■ 


(7) 


where v 0 ' is the line center or unshifted frequency, m 
is the mass of the atom, k is Boltzmann’s constant, 
and T is the translational temperature. Because 
hydrogen is the lightest atom, wide profiles are 
encountered. Temperatures in the arcjet vary from 
1000 to 5000 K, corresponding to half-widths (for 
each component) of 0.35 to 0.77 cm -1 . Neglecting 
the fine structure at 4000 K leads to an inferred 
temperature of 5500 K. This error increases at 
lower temperatures. 

The broadening of the Balmer a transition that is 
due to other mechanisms is expected to be small in 
comparison with Doppler broadening and has been 
neglected. The most significant of these is Stark 
broadening. The Stark-broadened half-width is a 
monotonically increasing function of electron num- 
ber density. 11 Emission scans of the Balmer (3 tran- 
sition (n - 4 to n = 2 at 486 nm) were made with a 
monochromator having entrance slits imaged onto 
the c enterline of the plume near the exit plane. A 
FWHM of 0.38 A (1.6 cm -1 ) was measured, which is 
approximately the spectral resolution of the mono- 
chromator. This width can be used to determine an 
upper limit of 2 x 10 13 cm -3 for the electron number 
density assuming a temperature of 5000 K (the exit 
plane temperature at the plume’s center). 11 This 
estimated upper limit is in agreement with a study of 
a similar arcjet operating on simulated hydrazine. 12 
The B alm er a transition is much less affected by 
Stark broadening. At this electron number density, 


Table 1. Components of the Balmer a Transition and 
Their Properties i Ref. 8) 


Upper j 

Lower/ 

8v(cm _1 ) 

Transition 

Relative 

Intensity 

1/2 

3/2 

-0.144 

3Si,2 2 P 3/2 

0.05 

3/2 

3/2 

-0.036 

3 D 3/2 2 P 3/2 

0.28 

5/2 

3/2 

0.0 

3Ds/ 2 — 2 P 3/2 

2.50 

1/2 

1/2 

0.220 

3Sj. 2 2 P\ 2 

3Pi,2 2 Si,2 

0.31 

3/2 

1/2 

0.328 

3 P 3, 2 2 S 1 . 2 

3D 3 2 2Pi. 2 

1.96 


20 October 1993 / Vol. 32. No. 30 / APPLIED OPTICS 6119 


the Balmer a transition is broadened with a Lorentz- 
ian line shape having a FWHM of only 0.01 cm -1 . 11 

Resonance broadening occurs when atoms in the 
upper or lower levels of the transition can undergo a 
resonant exchange of energy with a perturbing species. 
The broadening is Lorentzian in shape with a half- 
width proportional to the number density of the 
perturbing species and the oscillator strength of the 
resonant transition. 13 In the arcjet plume, the most 
populated species with a resonant energy level is 
ground-state atomic hydrogen. To get an estimate 
of resonance broadening, an estimate of the ground- 
state H-atom density is needed. For this purpose, 
we assume that the density, p, of the flow can be 
estimated from a quasi-one-dimensional conservation 
of mass. The mean mass-flow rate is given by 

m = f>vA, (8) 

where A is the nozzle exit area and v is velocity. If a 
velocity of 5 x 10 3 m/s and a dissociation fraction of 
30% are used, the H-atom number density is 10 22 
m“ 3 . If interactions between both the upper and 
lower levels of the Balmer a transition are considered, 
the Lorentzian resonance-broadened half-width is 
0.005 cm -1 . Within the framework of the above 
assumptions, this gives an estimate that is negligible 
compared with Doppler broadening. Nonresonant 
pressure broadening is also considered negligible 
because of the low pressures and high temperatures 
in the arcjet plume. 

Lifetime broadening arises as a result of the depopu- 
lation of both levels of the transition by radiative 
decay. 10 This broadening leads to a Lorentzian line 
shape. The half-width, which is calculated from the 
sum of all downward Einstein spontaneous emission 
coefficients from both the upper and the lower states 
involved in the transition, is given by 

Avz * * 2 tt ( j?) Aui + ^ 

For the Balmer a transition this summation yields a 
linewidth of 0.003 cm -1 . 6 

Transit-time broadening occurs when high mean or 
thermal velocities are present in a s m all fluorescence 
probe volume. The short interaction time between 
the atoms and the laser acts to decrease the lifetime of 
the fluorescence state. A characteristic time is deter- 
mined from the ratio of a characteristic length to the 
velocity. The characteristic length for this experi- 
ment is the laser beam waist of 0.1 mm. The peak 
velocity is 14 km/s. This gives a characteristic time 
of 7 x 10 ~ 9 s, which acts as a limiting lifetime. This 
would add 0.005 cm -1 to the linewidth. Hence this 
mechanism can also be neglected. 

At high power levels, the laser can significantly 
perturb the populations of the atomic states and 
thereby distort the absorption line-shape function. 
This effect is known as saturation broadening or 
power broadening. The induced absorption reduces 
the number of atoms in the lower (n = 2) state while 


increasing the number of atoms in the upper state 
(n = 3). If the number of atoms involved in this 
process is small compared with the total number of 
atoms in these states, the fluorescence signal is linear 
with the laser intensity and the line shape is indepen- 
dent of the laser intensity. When the laser intensity 
becomes large, the populations of the two states 
change significantly. This change is limited to a 
condition known as saturation, where the maximum 
ratio of populations of the two states is equal to the 
ratio of the degeneracies. As saturation is ap- 
proached, the fluorescence signal becomes a nonlin- 
ear function of the laser intensity, causing the exci- 
tation line-shape profile to become broader. To 
determine precisely the line centers and line shapes, 
saturation has been quantified and avoided. 

Other artifacts of the experiment can act to distort 
the shape of the measured excitation spectrum. The 
most significant of these are due to reductions of the 
spatial and temporal resolutions. If the velocity at a 
point changes in a time period shorter than the 
measurement time (several minutes), the measured 
profile will be a convolution of the varying Doppler 
shift and the instantane ous line shape. Ripple in the 
power supply that drives the arc discharge might lead 
to a modulation of the arcjet velocity. Because the 
temperatures obtained burn vertical excitation (with 
no mean velocity) are the same as those obtained from 
axial excitation (with high mean velocity) it appears 
that variations in velocity are not a problem. Poor 
spatial resolution can lead to the examination of 
atoms with a finite range of mean velocities. The 
superposition of the differently shifted fluorescence 
si gn a l s would give rise to an apparent increased 
broadening. Good spatial resolution was main- 
tained by keeping the diameter of the fluorescence 
probe volume small so that the range of mean veloci- 
ties was negligible and by ensuring that the collection 
optics were well focused. 

Experiment 

The arcjet thruster used in this experiment is a 
1-kW-class radiatively cooled thruster designed and 
built at the National Aeronautics and Space Adminis- 
tration Lewis Research Center. The tungsten nozzle, 
as shown in Fig. 1, has a 0.635-mm diameter, a 
0.25-mm-long constrictor, and a conical (20° half- 
angle) diverging section with an area ratio of 225 
(9.53-mm exit diameter). A more complete descrip- 
tion of the arcjet is available elsewhere. 14 Operating 
conditions for this experiment are 13 mg/s of H 2 
propellant with a maximum power dissipation of 1.48 
kW (voltage and current levels of 139 V and 10.3 A, 
respectively). The arcjet is operated in a 0.56-m- 
diameter cylindrical stainless-steel chamber 1.09 m 
long. The 0.35-Torr background pressure is main- 
tained by two 1250-CFM blowers evacuating through 
15-cm pipe. 

A schematic of the experimental setup is shown in 
Fig. 3. The fluorescence excitation laser is a cw ring 
dye laser (Spectra Physics Model 380) pumped by a 
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6-W argon-ion laser and operating on DCM dye. 
The laser produces several hundred milliwatts of 
power at 15,233 cm -1 (656 nm) with a no minal 
bandwidth of 10 MHz. Because this bandwidth is 
significantly less than the half-width of the measured 
fluorescence excitation spectrum of the Balmer a 
transition, the shapes of the recorded line accurately 
depict the broadening mechanisms. Part of the beam 
is passed through a 2-GHz Fabry-Perot interferom- 
eter to monitor mode structure and ensure single- 
mode operation. Another portion of the beam is 
directed into a wavemeter (Burleigh WA-10) that uses 
a Michelson interferometer to measure laser fre- 
quency to within 0.01 cm -1 . A third beam splitter 
directs the laser to a photodiode detector that moni- 
tors laser power. 

A beam splitter divides the laser into two beams of 
approximately equal intensity. One beam is for axial 
excitation and the other is for radial excitation. The 
beams are expanded and then focused through win- 
dows to a waist of 0.1 mm in the arcjet plume. The 
radial excitation beam is passed vertically through 
the plume. The axial beam is passed through a 
window in the end of the tank, and the thruster 
nozzle is used as a beam dump. The fluorescence 
signal is collected on a horizontal axis normal to both 
excitation beams. It is collected through an 83-mm- 
diameter window with a lens operating at f/5.7 (76 
mm diameter, 280 mm focal length) and imaged onto 
a R928 photomultiplier tube (PMT) with a magnifica- 
tion of 2.7. A notch filter with a nominal 10 A (23 
cm -1 ) bandwidth is placed directly in front of the 
PMT to filter extraneous light. A 0.4-mm aperture 
is mounted directly in front of the filter at the focal 
plane of the imaging lens to act as a field stop and 
determine the measurement location in the plume as 
well as the axial and vertical spatial resolution. 
Phase-sensitive detection is used to discriminate the 
fluorescence signal from the intrinsic plasma emis- 
sion, which is 30 to 500 times greater. The lock-in 


amplifier (SRS SR530) is synchronized to a mechani- 
cal chopper at 3 kHz. The output from the lock-in is 
collected and stored in a laboratory computer along 
with the output from the wavemeter and laser power 
detector. 

At each spatial point investigated, the laser is 
automatically scanned in frequency by varying the 
voltage on a (piezo-controlled) 75-GHz intracavity 
etalon. Radial and axial excitation are obtained on 
separate scans by hlnnking one of the beams. The 
laser is able to sran over a range of up to 2 cm -1 . 
This is enough to capture most of the fluorescence 
excitation line. Typical scan times are a few min- 
utes. 

Both the arcjet and a stationary hydrogen dis- 
charge were used as a stationary (unshifted) reference. 
In the former case, it is assumed that the flow is 
axisymmetric, which implies that the radial velocity 
is zero at the centerline. The symmetry of both the 
axial and r adial velocity profiles confirms this assump- 
tion. The wavemeter provides enough accuracy and 
repeatability to forego simultaneous reference scans. 
The stationary discharge was produced in a low- 
pressure hydrogen cell mounted in a microwave 
cavity, which acted as a source of slow (essentially 
stationary) H atoms. 

The velocity and temperature are determined by a 
least-squares fit to the line-shape model with the 
Levenberg-Marquardt method. 15 The model con- 
sists of five components with specified relative inten- 
sities and separations. Doppler broadening is the 
only broadening mechanism considered. There are 
three parameters that are used to fit the data to the 
model. These are the Gaussian half-width, the line 
shift, and the integrated intensity. The tempera- 
ture is deter min ed from the Gaussian half- width. 
The velocity is determined from the shift of the 
fluorescence excitation profile relative to that of the 
reference spectrum. The integrated intensity was 
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necessary to fit the other parameters but can also be 
used to determine relative excited H atom densities. 

Results and Discussion 

The effect of saturation was evaluated by comparing 
UF line profiles for a range of laser-power levels. 
For the 0.1-mm laser spot size, the fluorescence was 
found to be linear below 60 mW. Thus all data were 
taken with laser powers well below this level. 

Two typical fluorescence excitation scans are shown 
in Fig. 4. The scan on the left is taken from the 
plume centerline 0.4 mm from the exit plane with 
axial excitation where the measured temperature is 
4500 K. The scan on the right is from the colder, 
600 K stationary reference discharge. At tempera- 
tures below approximately 1000 K, the Balmer a line 
shape begins to look more like a doublet. This is due 
to the fine structure discussed in the Experiment 
section. The solid curves in Fig. 4 correspond to the 
model’s fit to the data. The fluorescence intensity is 
the ratio of the output signal from the lock-in ampli- 
fier to the signal from the photodiode that monitors 
laser power. Accounting for laser power in this way 
preserves the line shape because the fluorescence is 
linear and not saturated. The power varies by up to 
10% as the laser scans in frequency. Both profiles 
were normalized by their areas to facilitate compari- 
sons of the line shapes between the hot and cold 
gases. These profiles consist of 70 to 200 data points 
each, with the data taken at 1-s intervals. Profiles 
from regions further from the centerline require 
longer scan times because of decreased signal-to- 
noise ratios. 



Laser Frequency (cm -1 ) 

Fig. 4. Sample data taken near the centerline 0.4 mm from the 
exit plane with the fluorescence excitation spectrum from axial 
excitation and from the (unshifted) reference celL The 0.63 cm -1 
shift corresponds to a 12.4 km/s velocity. Both profiles are 
normalized by their areas after subtraction of the baseline. 


Two laser probe b eams are focused to the same 
small volume on the plume centerline 0.4 mm from 
the exit plane. The two beams have respective propa- 
gation vectors aligned along the radial and axial 
directions of the arcjet plume. Data were taken as 
the arcjet was positioned at a variety of vertical 
locations. Figure 5 shows the radial profile of the 
axial velocity. The profile is highly peaked at 14 
km/s with a fairly linear gradient away from the 
peak. Note that the profile is shifted approximately 
0.7 mm from the indicated centerline (0 mm radial 
location). This is most likely due to the difficulty in 
establishing the plume centerline, rather than devia- 
tion of the arcjet plume from axisymmetric flow. 
In establishing the plume centerline the axial laser 
beam is aligned by centering it on the cathode tip, 
which is more than 13 mm beyond the exit plane. 
If the arcjet nozzle is not coaxial with the laser, the 
laser will pass through the exit plane at an off-center 
location. A 0.7-mm deviation at the exit plane corre- 
sponds to a 3° misalignment of the laser and arcjet. 

The radial profile of the radial velocity component 
is shown in Fig. 6. The profile has a fairly linear 
shape and, like the axial velocity, is symmetric about 
the centerline. This validation of our symmetry 
assumption means that radial excitation at the center- 
line provides a reliable unshifted reference fluores- 
cence profile. The magnitude of the radial velocity 
reaches a maximum of 4 km/s. This occurs near the 
edge of the nozzle, where the axial velocity is approxi- 
mately the same magnitude, giving a velocity direc- 
tion of 45° (from the arcjet axis). The flow direction 
is more clearly represented in Fig. 7, where the radial 
and axial velocities have been combined in a velocity 
vector format. Note that the axial velocity data from 
Fig. 5 were shifted to correct for the arcjet axis 
alignment problems. Figure 7 best shows the small 
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Fig. 5. Radial profile of axial velocity 0.4 mm from the exit 
plane. The exit plane radius is 4 .8 mm. 
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Radial Position (mm) 

Fig. 6. Radial profile of the velocity in the vertical direction along 
a vertical line through the arcjet center, 0.4 mm from the exit 
plane. 

asymmetry near the edges of the plume. This is 
most likely due to tank effects. The tank pump port 
is only 375 mm downstream from the arcjet exit plane 
and off to the side (280 mm from the centerline), 
allowing a nonaxisymmetric pressure field to develop. 

Measurements were made of the velocity compo- 
nent in the angular direction about the axis (circum- 
ferential or swirl velocity). The swirl velocity is 
normal to both the radial and axial directions. To 
perform these measurements, the vertical probe beam 
was translating away from the plume centerline such 
that it remained normal to the radial velocity compo- 
nent. Velocities up to 1 km/s were measured near 
the edge of the plume. However, this is inconclusive 
because a minor misalignment of the arcjet with the 
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Fig. 7. Radial and axial velocities, combined to form a radial 
profile of the velocity vectors 0.4 mm from the exit plane. 


laser axis can account for this result. 11 tne laser is 
not completely normal to the radial axis, there wifi be 
a nonzero component of the radial velocity in the 
direction of laser propagation. This problem can be 
overcome by taking velocity measurements through- 
out the fiowfield along a two-dimensional grid. The 
radial and swirl components can then be separated. 
The significance of the current results is that the 
swirl velocity components are apparently small com- 
pared with the radial or a vial velocity components. 

The temperature is determined by fitting the data 
to the model as discussed above. In applying the 
fitting procedure, it was necessary to account for a 
positive baseline, which was caused by the scattering 
of laser light from surfaces. The temperature ob- 
tained from the model is sensitive to small changes in 
this baseline. The precision in determiningthe base- 
line can therefore limit the uncertainty of the tempera- 
ture measurement The large linewidths at the 
higher temperatures and limited scanning range of 
the laser prevent a complete scan of the entire 
fluorescence profile. Therefore, to establish the base- 
line on one side of the profile, a portion of the profile 
on the other side is truncated. This forces a compro- 
mise in the ability to establish the baseline and the 
ability to capture the width of the profile. At lower 
temperatures (where the densities are also lower) 
uncertainties are increased by noise in the baseline 
combined with high baseline-to-peak fluorescence 
ratios. At the radial edges of the plume, this ratio is 
as high as 70%, whereas at the centerline it is as low 
as 1%. 

This sensitivity to the baseline is illustrated in Fig. 
8, which shows fluorescence data along with three 
profiles obtained by fits to the model for different 
baseline values. The program that fits the data to 
the model can fit to a specified baseline, or it can 
determine a baseline by a least-squares fit. When 
the program determines the baseline by a fit, it 
overpredicts the baseline to compensate for an under- 
prediction of the wings. That is, the program deter- 
mines a baseline that is dearly higher than the actual 
baseline. This produces a temperature on the low 
side. When the baseline is forced to zero, the fitting 
procedure produces a temperature that is signifi- 
cantly high er. The most appropriate baseline is 
somewhere between these limits. For the results 
presented here, the baselines were all determined 
manually using interactive graphics software. For 
the fluorescence profile in Fig. 8, the two extreme 
baselines produce temperatures that differ by 900 K, 
or 20% of the determined temperature. If any signifi- 
cant portion of the left-hand side of this profile were 
omitted, a baseline could not be established. Alterna- 
tively, if the laser were timed to begin scanning at a 
lower frequency, it could not scan enough of the 
profile to establish the width. All three of the fitted 
profiles underpredict the fluorescence in the wings. 
This error could be caused by Lorentzian broadening 
mechanisms, which are not included in the model. 
Lorentzian line shapes decay much more slowly in the 
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Fig. 8. Sample of raw data showing fits produced by the model. 
The three lines were fixed to different baselines and correspond to 
determined temperatures of 3880, 4306, and 4775 K_ This figure 
illustrates the sensitivity of temperature determination to correct 
baseline establishment. 


wings compared with Gaussian line shapes. There- 
fore, a small Lorentzian component to the broadening 
would first affect the wings. Because the effect is 
seen only in the wings, the overall error in the 
inferred temperature is small. To estimate the maxi- 
mum error this might cause, a profile was generated 
by convolving a Lorentzian profile with a half-width 
of 0.01 cm -1 corresponding to an electron density of 
lO 13 cm" 1 and a Gaussian profile with a 0.70 cm -1 
half-width corresponding to a temperature of 4200 K 
This profile (which included the effects of the fine 
structure) was then fitted to the model used for data 
reduction. The Lorentzian component led to an 
overprediction of the temperature by 300 K, which is 
less than the uncertainty that was due to baseline 
establishment. The misfit in the wings is only pre- 
sent near the plume centerline, where the electron 
number density is likely to be the highest. Toward 
the plume edges, the fluorescence signal is weaker 
and the baseline is noisier, making the baseline 
equally difficult to locate. Examples of both of these 
phenomena can be seen in Fig. 8. 

A radial profile of the temperature is shown in Fig. 
9. The temperature peaks between 4000 and 5000 K 
at the centerline and falls to 1000 K at the edges of 
the plume. The shape is similar to the radial profile 
of the axial velocity with the same shift in the 
centerline. The scatter is relatively high because of 
the large uncertainties. The data presented were 
obtained from a mix of radial and axial laser excitation. 



Vertical Position (mm) 

Fig. 9. Radial profile of translational temperature. 


The overlap in the data indicates that there are no 
anisotropies or fluctuations in velocity. 5 Note that 
fluctuations of temperature give rise to the measure- 
ment of an average temperature. Fluctuations of 
velocity give rise to the measurement of the convolu- 
tion of bulk velocity shifts and the instantaneous line 
shapes, leading to an erroneously high temperature. 
Because there is no radial velocity at the centerline, it 
is unlikely that there are any significant velocity 
fluctuations. Because the temperatures measured 
from radial and axial excitation are the same it 
follows that there are no significant fluctuations in 
the axial velocity. 

The integrated intensity of the fluorescence excita- 
tion spectrum is a function of the number density of 
hydrogen in the absorbing state (n = 2). The radial 
profile of the integrated signal from data near the exit 
plane is shown in Fig. 10. The data are plotted on a 
log scale to show that the fluorescence signal drops by 
2.5 orders of magnitude from the plume center to the 
edges. This is most likely due to decreases in both 
the Boltzmann fraction of the excited state as well as 
the dissociation fraction. This follows from the 3000 
K decrease in temperature. If the quenching rate is 
small compared with radiative relaxation, the inte- 
grated fluorescence signal is directly proportional to 
the number density of H in the absorbing state. The 
large differences in energy between the upper fluores- 
cence state and the adjacent states lead to small 
quenching cross sections, the low densities further 
decrease quenching rates, and the radiative lifetime 
of the fluorescence is short. However, without fur- 
ther knowledge of the collisional rates, we cannot 
make any quantitative conclusions about number 
densities. 

Measurements were taken 0.4 mm from the exit 
plane at the centerline for a variety of arcjet power 
levels and a mass flow rate of 13 mg/s. The axial 
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Vertical Position (mm) 

Fig. 10. Integrated fluorescence signal intensity profile, 0.4 mm 
from the exit plane. 


velocities and temperatures are shown in Fig. 11. 
Note that when the input power is more than doubled, 
the centerline axial velocity increases by less than 
30%. This is consistent with the variation of thrust 
with input power. 16 However, for the same change 
in power, the kinetic temperature increases approxi- 
mately 3 times. This indicates an increased level of 
enthalpy that is converted out of the nozzle. This 
represents lost energy that is not converted to kinetic 
energy or thrust. This is consistent with a drop in 
the efficiency as calculated from thrust measure- 
ments. 16 The efficiency is defined as the ratio of I SP 
to the sum of the specific power and I SP that is 
produced with the arc off. 
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Fig. 1 1. Axial velocity and temperature at the centerline. 0.4 mm 
from the exit plane, for a variety of power levels. 



When the arcjet is operating with argon propellant, 
a shock structure that is due to the mismatch be- 
tween exit and background pressure is clearly visible 
in the varying emission intensity. However, when 
the arcjet is operating on hydrogen propellant, no 
such structure is visible. There could be several 
reasons for this, including a higher mean free path 
leading to relatively thick or diffuse shock waves; a 
more closely matched pressure at the exit plane; or 
more oblique shock angles, making radial viewing of 
the change in emission (because of the temperature 
jump) difficult. Measurements of the axial velocity 
and temperature were taken at a variety of locations 
along the centerline. Fig. 12 shows these results 
along with the Mach number as determined from the 
data (assuming single-component -flow, see conclu- 
sions). The Mach number rises steadily from 1.4 at 
the exit plane to 3.0 at approximately 18 mm (2 exit 
diameters) downstream. After going through the 
shock, the Mach number drops. Note that the flow 
does not appear to be accelerating in the expansion 
region of the first 18 mm. The velocity is slowly 
decreasing until the Mnrh number peaks where the 
velocity begins to decrease more rapidly. The de- 
creasing velocity in the expansion section is in part 
due to an off-axis probe volume. The measurements 
were m ade along the indicated 0-mm radial position. 
As shown in Fig. 5, this is not the true location of the 
plume centerline. This deviation from the true cen- 
terline position is likely to increase further from the 
exit plane. Thus the velocity shown in Fig. 12 might 
be a component of and not the full centerline velocity. 
This could be corrected by taking data at various 
radial positions for each axial location. Although 
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Fig. 12. Axial velocity and temperature taken from a scan along 
the centerline with constant 1.43-kW input power. The Mach 
number in the lower graph is derived from the velocity and 
temperature data and shows the shock at 18 mm. 
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these gasdynamic property variations are more a 
function of the background pressure than the plasma 
dynamics in the nozzle, they do suggest that the flow 
is underexpanded. Under these circumstances, the 
exit pressure must be higher than the 0.35-Torr 
chamber pressure. 

A few simple experiments were conducted to deter- 
mine the presence of polarization of the fluorescence. 
The fluorescence process is an induced electric dipole 
phenomenon. 10 Unless there are a significant num- 
ber of collisions to randomize the orientation of the 
fluorescing atoms, the fluorescence will be emitted in 
directions and with polarization reflecting the polar- 
ization and orientation of the excitation laser. The 
laser polarization was rotated while monitoring the 
variation in fluorescence intensity and line shape. 
In addition, a polarizing filter was placed in front of 
the PMT to analyze the fluorescence for preferred 
polarization. No polarization effects were detected. 
This indicates that the time between collisions involv- 
ing excited H atoms is small compared with the 
lifetime of the transition. 

Conclusions 

The mean mass exit velocity of the hydrogen arcjet as 
determined from previous thrust measurements (at 
the 1.48-kW power level) is 8.3 km/s. 16 This appears 
consistent with the measured axial atomic velocities. 
Note that although the 14 km/s peak velocity is 
significantly higher than the mean mass velocity, the 
axial velocity drops to 4 km/s at the plume edges. 
It is likely that the density is not the highest at the 
centerline, where the temperature and dissociation 
fractions are hi gh . This interpretation is supported 
by mass flux profiles determined by sampling probes 
in similar devices. 17 To fully correlate thrust mea- 
surements to velocity data, a measurement of the 
total density is required as well as a measurement of 
the (possibly different) velocity distribution of molecu- 
lar hydrogen. The high radial velocities at the plume 
edge may represent a significant thrust loss. This 
effect also would be better quantified by density 
measurements. 

Velocity measurements far from the exit plane (120 
mm) have been made with molecular beam sampling 
mass spectrometry. 18 These measurements suggest 
the presence of high diffusi on velocities for atomic 
hydrogen. It would be valuable to measure directly 
the velocity of molecular hydrogen to determine if 
this diffusion is significant at the exit plane where the 
pressure ( > 0.35 Torr) is several orders of magnitude 
higher than the 5 x 10“ 4 Torr near the molecular 
beam sampling skimmer. There are a multitude of 
transitions between the excited states of molecular 
hydrogen that are accessible by cw dye lasers. 
Future work will include a feasibility study of these 
transitions for velocity measurements. 

At the high peak temperatures measured, the 
central portion of the plume should consist predomi- 
nantly of atomic hydrogen. Because the flow chemis- 


try may be partially frozen, equilibrium calculations 
will give lower bound values of the dissociation 
fraction. For hydrogen at 8 Torr, the equilibrium 
composition will vary from one-fourth atomic hydro- 
gen to one-fourth molecular hydrogen over the range 
2500 to 3000 K. (Although the exit plane pressure is 
unknown, an order-of-magnitude change in pressure 
shifts this dissociation transition range by approxi- 
mately 300 K.) We conclude that the central area of 
the plume consists predominantly of atomic hydrogen. 
The dissociation fraction outside this core is more 
dependent on the extent to which the flow is frozen. 
Because the gas at the centerline consists predomi- 
nantly of atomic hydrogen, the peak atomic hydrogen 
velocities measured are the same as the peak mean 
mass velocity. 

This paper has demonstrated the utility of spec- 
trally resolved LIF to measure velocity and transla- 
tional temperature in an arcjet plume. The measure- 
ments have provided a variety of information on the 
nature of the arcjet plume. In summary, we find 
high peak temperatures, indicating a high dissocia- 
tion fraction, low swirl velocities, and the presence of 
diffuse shock structure within 2 nozzle diameters 
from the exit plane. Even though the arqjet is 
axisymmetric it would be desirable to make measure- 
ments of the complete flow field. If measurements 
are taken throughout a plane at each axial location, 
problems that are due to laser— arcjet ali gnm ent are 
eliminated. This would be important for a more 
rigorous characterization of the arqjet plume. 

This research was supported by a grant from the 
National Aeronautics and Space Administration Lewis 
Research Center, and Frank Curran was the grant 
monitor. 

References 

1. P. J. Wilbur, R. G. Jabs, and F. M. Curran, “Space electric 
propulsion plasmas,” IEEE Trans, on Plasma ScL 19, 1167- 
1179(1991). 

2. M. A. Cappeili, J. G. Liebeskind, R. K. Hanson, G. W. Butler, 
and D. Q. King, “A comparison of arqet plume properties to 
model predictions,” paper AIAA 93-0820 in 31st Aerospace 
Sciences Meeting (American Institute of Aeronautics and 
Astronautics, New York, 1993). 

3. D. A. Erwin, G. C. Pham-Van-Diep, and W. D. Deininger, 
“Laser-induced fluorescence measurements of flow velocity in 
high-power arqet thruster plumes,” AIAA J. 29, 1298-1303 
(1991). 

4. J. G. Liebeskind, R. K. Hanson, and M. A. Cappeili, “Velocity 
measurements in a arqet using LIF," paper AIAA 91-2112 in 
27th Joint Propulsion Conference (American Institute of Aero- 
nautics and Astronautics, New York, 199 1 ). 

5. W. M. Ruyten and D. Keefer, “Characterization of electric 
thruster plumes using multiplexed laser induced fluorescence 
measurements” paper AIAA 92-2965 in 28th Joint Propulsion 
Conference (American Institute of Aeronautics and Astronau- 
tics, New York, 1991). 

6. C. E. Moore, “Atomic energy levels,” Natl. Bur. Stand. (U.S.) 
Circ. 467, 1-8 (1949). 

7. R. P. Lucht, “Applications of laser-induced fluorescence spec- 
troscopy for combustion and plasma diagnostics,” in Laser 


6126 APPUED OPTICS / Vol. 32, No. 30 / 20 October 1993 


Spectroscopy and Its Applications, L. J. Radziemski, R W. 
Solarz, and J. A. Paisner, eds. (Dekker, New York, 1987), 
Chap- 9, p. 627. 

8. H. A. Bethe and E. E. Salpeter, Quantum Mechanics of One 
and Two-Electron Atoms (Springer- Verlag, Berlin, 1957), Sec. 
64, p. 274. 

9. W. G. Vincenti and C. H. Kruger, Introduction to Physical Gas 
Dynamics (Krieger, Malabar, Fla., 1961), Chap. 9, p. 383. 

10. W. Demtroder, Laser Spectroscopy (Springer-Verlag, Berlin 
1981), Chap. 3, pp. 85-87. 

11. C. R Vidal, J. Cooper, and E. W. Smith, "Hydrogen Stark- 
broadening tables,’’ Astrophys. J. Suppl. Ser. 214, 37-136 
(1973). 

12. D. H. Manzella, F. M. Curran, R. M. Myers, and D. M. Zube, 
“Preliminary plume characteristics of an arcjet thruster,” 
paper ALAA 90-2645 in 21st International Electric Propulsion 
Conference (American Institute of Aeronautics and Astronau- 
tics, New York, 1990). 

13. H. R Griem, Plasma Spectroscopy (McGraw-Hill, New York, 
1964), Chap. 4, pp. 95-98. 

14. F. M. Curran and T. W. Haag, “An extended life and perfor- 


mance test of a low power arget,” paper AIAA 88-3 106 in A 

Joint Propulsion Conference (American Institute of Aeronau- 
tics and Astronautics, New York, 1988). 

15. W. H. Press, B. P. Flannery, S. A. Teukolaky, and W. T. 
Vettling, Numerical Recipes in C (Cambridge U. Press, Cam- 
bridge, 1988) Chap. 14, pp. 528-547. 

16. F. M. Curran, S. R. Bullock, T. W. Haag, C. J. Sarmiento, and 
J. M. Sankovic, "Medium power hydrogen arcjet operation,” 
paper AIAA 91-2227 in 27th Joint Propulsion Conference 
(American Institute of Aeronautics and Astronautics, New 
York, 1991). 

17. W. M. Van Camp, D. W. Esker, R J. Checkley, W. G. Duke, 
J. C. Kroutii, S. B. Merrifield, and R A. Williamson, Study of 
Arc-Jet Propulsion Devices, Rep. E368 (McDonnell Aircraft 
Corporation, St Louis, Mo., 1966); Rep. CR-54691 (National 
Aer onautic s and Space Administration, Cleveland, Ohio, 1966). 

18. R. P. Welle, J. E. Pollard, S. W. Janson, M. W. Crofton, and 
R- B. Cohen, “One kilowatt h ydr ogen and helium arcjet 
performance,” paper AIAA 91-2229 in 27th Joint Propulsion 
Conference (American Institute of Aeronautics and Astronau- 
tics, New York, 1991). 


20 October 1993 / Vol. 32. No. 30 / APPLIED OPTICS 6127 


Experimental Investigation of Velocity Slip near 
an Arcjet Exit Plane 

J. G. Liebeskind, R. K. Hanson, M. A. Cappelli 



A publication of the 

American Institute ot Aeronautics and Astronautics. Inc. 
370 L'Enfant Promenade. SW 
Washington. DC 20024-2518 


AlAA JOURNAL. VOL. 33. NO. 2: TECHNICAL NOTES 


373 


PRECEDING PAGE BLANK NOT FILMED 


Experimental Investigation 
of Velocity Slip near an 
Arcjet Exit Plane 


J. G. Liebeskind.* R. K. Hanson,' and M. A. Cappellr 
Stanford University. Stanford. California 94305 

Introduction 

V ELOCITY slip is a phenomenon associated with multicompo- 
nent flows in which species of different mass have different 
mean velocities. This phenomenon has been studied in molecu- 
lar beam nozzles' and has been exploited for isotope separation. 2 
Significant velocity slip has been measured in the plume of low- 
power arcjet thrusters.’ These recent measurements, obtained by 
mass spectrometry of a molecular beam sampling probe, were lim- 
ited to locations more than 10 diameters from the exit plane. 

The determination of whether slip develops in the nozzle, near 
the exit plane, or in the expansion plume of an arcjet thruster, is 
important to both modeling efforts and to the interpretation of mea- 
surements. If velocity slip develops inside the nozzle, numerical 
models will need to include separate momentum equations for each 
neutral species or. at the very least, an appropriate model for multi- 
component transport. These additions increase the complexity of 
an already difficult task. Measurements of exit-plane velocity in 
these devices have recently been made by laser-induced fluores- 
cence (LIF). 4 a species specific technique. If significant velocity 
slip exists, then separate measurements for each species are neces- 
sary to completely characterize the fiowfield. 

We have previously reported LIF-based measurements of atomic 
hydrogen velocity and kinetic temperature in a hydrogen-fueled 
arcjet. 4 5 We recognize, however, that the flow consists of both 
atomic and molecular hydrogen (along with small fractions of ions 
and electrons). The most direct method to investigate velocity slip 
in the arcjet fiowfield would be to measure, in addition to the veloc- 
ity of atomic hydrogen, the velocity of molecular hydrogen under 
the same conditions. As with atomic hydrogen, absorption transi- 
tions from the ground state require vacuum-ultraviolet wavelengths. 
Owing to the low densities of excited-state molecules and the ro- 


vibrational distribution which reduces the density of any particular 
state, these transitions are difficult to probe. 

A simpler approach, for the purpose of evaluating slip, is to seed 
the flow with a species that is accessible with the same laser used 
to probe atomic hydrogen. In the present study, helium was chosen 
as the seed species owing to its inertness (with respect to the arcjet 
nozzle), its relative mass (four times the mass of atomic hydrogen), 
and its convenient electronic transitions in the visible wavelength 
region. Velocity and temperature are measured by LIF of both he- 
lium and atomic hydrogen at the same arcjet operating condition. 
Absence of slip between helium and atomic hydrogen would sug- 
gest that slip is not a dominant mechanism in the nozzle or exit plane 
vicinity of our hydrogen arcjet. 

Theory 

In a multicomponent mixture, each species (denoted by i). may 
have its own mean velocity v,-, and the mean mass velocity v m is 
defined by 


Vn, 


YL n ' m < v i 


(i) 


Here. /?,- is the number density of species i and m, is the mass. The 
simplified steady-state momentum equation for a single species in 
the axial direction c can be expressed as 


dvi Bp, 
+ — 


mm,- 




( 2 ) 


where p, is the partial pressure. The term on the right side of Eq. (2) 
represents the momentum exchange between species of different 
types due to collisions. The momentum exchange term formulated 
in this way assumes that the force per unit volume exerted on par- 
ticles of species i. due to collisions with species of type j. is pro- 
portional to the difference between their respective mean velocities. 
The proportionality constant k,j. which accounts for differences in 
mass of the colliding panicles, can be considered a collision fre- 
quency for momentum transfer. This constant can be determined 
from kinetic theory. 6 - 7 The slip velocity t)j, p is usually defined for a 
binary mixture as the difference in velocities. • 


Received March 12. 1994: revision received Aug. 10. 1994; accepted 
for publication Aug. 19. 1994. Copyright © 1994 by the authors. Pub- 
lished by the American Institute of Aeronautics and Astronautics. Inc., with 
permission. 

■Research Assistant. High Temperature Gasdynamics Laboratory. De- 
partment of Mechanical Engineering. Student Member A1AA. 

^Chairman. High Temperature Gasdynamics Laboratory. Department of 
Mechanical Engineering. Associate Fellow AlAA. 

•Assistant Professor. High Temperature Gasdynamics Laboratory. 
Department of Mechanical Engineering. Member AlAA. 


ujip = vj - v, (3) 

If the density is high such that flow is in a continuum regime, the 
collision rate will be high, allowing efficient momentum exchange 
between species. Thus, the different species will be in equilibrium, 
and the slip velocity will be negligible. However, in an expanding 
flow, conditions can occur in the transition from continuum to free 
molecular flow, where velocity and pressure gradients are large and 
the collision rates are low. Insufficient collisions for momentum 
exchange between species allows a slip velocity to develop. The 
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proportionality constant k y s is inversely proportional to the Knudsen 
number Kn (Ref. 7) 

k 0 oc \/Kn (4) 

As the Knudsen number increases within the continuum regime, 
momentum transfer between species decreases. However, as the flow 
becomes rarefied to the extent that collisions become negligible, the 
slip approaches a constant value (which may be zero) and Eq. (2) 
is no longer valid. The Knudsen number for this flow has been 
determined, based on experimental results, to be just under unity. 5 
This puts the flow regime in the transition between continuum and 
rarefied flow. 

Experiment 

The arcjet thruster used in this experiment is a 1-kW class radia- 
tively cooled thruster designed and built at NASA Lewis Research 
Center. The tungsten nozzle has a 0.64-mm-diam throat. 0-25-mm- 
long constrictor, and a conical (20-deg half-angle) diverging section 
with an area ratio of 225 (9.53-mm exit diameter). A more complete 
description of the arcjet is available elsewhere. 8 - 9 When operating on 
the hydrogen-helium mixture, the current and voltage are 9.5 A and 
143 V corresponding to 1 .36 kW arc power. The arcjet is operated in 
a 0.56-m-diam cylindrical stainless steel chamber 1 .09 m long. The 
0.5-Torr background pressure is maintained by two 35.4 m 3 /min 
coupled to the chamber with 152-mm-diam pipe. 

The temperature and velocity measurements were made using 
LIF. A continuous wave ring dye laser scans across the atomic 
transition. Since the laser bandwidth is much smaller than the spec- 
tral feature, the fluorescence excitation spectrum accurately depicts 
the broadening mechanisms of the probed species. The velocity is 
determined from the Doppler shift whereas the temperature is in- 
ferred from the shape of the fluorescence excitation spectrum using 
a Doppler-broadened line-shape model.' 1 Since Doppler broadening 
arises from the velocity distribution of the species, this diagnostic 
technique probes the kinetic or translational temperature which may 
differ from the electronic, rotational, and vibrational temperatures of 
the molecular species in the flow. 1011 Details of the implementation 
and analysis have been published previously 4 ; the only difference in 
this experiment is the use of optogal vanic detection in a dc discharge 
for the helium reference (unshifted) signal. The optogalvanic cell 
provided a reference for helium with a high signal-to-noise ratio. 
Fluorescence and absorption in a microwave discharge were used 
for a hydrogen reference calibration because it provided a stronger 
signal-to-noise ratio than the optogalvanic technique. 

The atomic hydrogen transition probed is the Baimer or transition 
at 656 nm (15.233 cm -1 ), an electronic transition between the first 
and second excited states of H. i.e.. n = 2 — *• n = 3. This transition 
actually consists of the overlap of several fine structure components. 
Since the separation of these components is of the same magnitude 
as the broadening, the fine structure must be accounted for when 
analyzing the line shape. The helium transition probed is the 2 1 P — *■ 
3’D transition at 668 nm (1 4.970 cm -1 ). 

The relatively high energy of helium's first excited state. 
166.278 cm -1 (20.6 eV) above the ground state can lead to weak 
fluorescence signals. To obtain a significantly strong fluorescence 
signal (relative to the noise), a high helium flow rate of 6.2 standard 
liters/min (SLPM ) was necessary. This is significant when compared 
to the 9.0 SLPM of hydrogen, and hence altered the normal oper- 
ating conditions of the arcjet. In particular, the addition of helium 
to the propellant increases the arc voltage and reduces the veloc- 
ities. The voltage increase of less than 3 V can be explained by a 
combination of the increased mass flow (due to the helium addition) 
and helium's higher ionization energy. The 30% velocity decrease 
(at the centerline) can be attributed to helium's higher atomic mass 
and the reduced specific power. Also, the peak temperature drops 
20% from 5000 K to 4000 K. Although the addition of the helium 
modifies the flow parameters, the conditions that may lead to slip are 
expected to be similar to those encountered in the unseeded g arcjet. 

Results and Conclusions 

The raw data from the helium and hydrogen fluorescence are 
shown in Figs. 1 and 2. respectively. The Doppler shifts corre- 
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Fig. 1 Doppler-shifted fluorescence, excitation spectrum of helium 
along with the unshifted reference signal used to calibrate the wave 
meter, probe volume centered in the plume, 0.5 mm from the exit plane. 
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Fig. 2 Doppler shift of atomic hydrogen measured relative to its cali- 
bration, measurement location same as the helium data; note the smaller 
atomic mass leads to broader line widths even though the temperature 
is the same. 

spond to axial velocities of 10.4 km/s for helium and 10.3 km/s 
for hydrogen. The magnitude of these measurements is consistent 
with previous measurements, 4 - 5 thrust measurements. 8 and model 
predictions. 12 The hydrogen and helium velocities agree within the 
0.2 km/s uncertainty of the measurement. This indicates that the slip 
velocity is less than the uncertainty of the measurement; i.e.. the up- 
per limit of the slip velocity is about 2% of either species’ velocity. 

The inferred temperature for helium and hydrogen are 3990 K 
and 4260 K, respectively. The 6% difference is less than the 20% 
uncertainty of the measurement. These results allow us to conclude 
that all species share a common kinetic temperature. Because of the 
higher sensitivity to noise, the uncertainties of the temperature are 
higher than the uncertainties of the velocity. 

The results indicate that the velocities of atomic hydrogen and he- 
lium are the same, suggesting that in the present device, the velocity 
of any single species, as measured in the vicinity of the exit plane, 
closely represents the mean mass velocity (within the uncertainty of 
the measurement). This conclusion is significant because it implies 
that species specific velocity-measurement techniques such as LIF 
need only be applied to a single species. In a hydrogen arcjet, it is 
much easier to probe the atomic hydrogen than the molecular hydro- 
gen. furthermore, it is often easier to probe electronic excited states. 

For further investigation of the effect of collisions, it would be 
interesting to extend these slip measurements further along the axis 
of the plume to see where slip develops. This extention. however, 
would require a test chamber capable of maintaining high vacuum 
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back pressures that would better reproduce the thruster plume struc- 
tures expected in flight conditions. 
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Abstract 

This paper describes a detailed comparison of experimentally measured arcjet plume properties to 
model predictions. The experiments focus on the measurement of the exit plane vector velocity field 
and kinetic temperature by laser induced fluorescence. The model is based on a single-fluid 
magnetohydrodynamic description of the flow. Reasonable agreement between model predictions and 
experimental exit plane velocity profiles is obtained over a wide range of specific power, despite the 
complex nature of the flow. At higher powers however, there are significant differences between the 
measured and predicted exit plane temperature. At the highest specific power studied, the model 
overpredicts the temperature by as much as 2000K and the specific impulse by approximately 20%. 
These differences cannot be attributed to the present treatment of the collisionless anode sheath or 
finite-rate molecular recombination (dominant frozen flow loss) and is believed to be more closely 
related to the treatment of other physical processes such as finite-rate ionization, ambipolar diffusion 
or possibly the selection of exit plane boundary conditions. These findings suggests that further 
refinements in the treatment of these processes may be necessary in order to extend the predictive 
capabilities of the model. 


L Introduction 

Arcjet thrusters are high specific impulse 
electric propulsion devices that are presently 
targeted for use in stationkeeping and orbit 
transfer.^ In comparison to other propulsion 
alternatives such as low-thrust chemical rockets, 
the use of arcjets will enable smaller launch 
vehicles for payload delivery because of the 
significant savings in propulsion system mass. 

Much emphasis has been placed recently on the 
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development and flight qualification of low power 
(<2 kW) hydrazine arcjet thrusters.^ However, in 
order to meet a wider range of mission possibilities, 
arcjet operating ranges must be extended to higher 
powers (2 - 30 kW). Although much progress has 
been made in this direction through extensive 
laboratory studies, 5-7 it is now evident that scaling 
to higher powers will require a better understanding 
of the physical processes that control arcjet 
performance. These include: (i) the physics of arc 
root attachment at both the cathode and the anode; 
(ii) the processes that establish the local plasma 
conductivity (i.e., ionization, ambipolar diffusion of 
charged species); (iii) the dissipation of energy via 
Joule heating, the subsequent transport of this 
energy to the internal surfaces of the arcjet. and the 
conduction of energy away from these internal 
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surfaces; (iv) the magneto-hydrodynamic coupling 
of the flowfield to the electrical nature of the arc, 
and; (v) processes such as dissociation and 
vibrational excitation, which give rise to frozen 
flow losses and reduce thrust efficiencies. 

To better understand these processes, studies 
that include a combination of modeling and 
laboratory measurements are necessary. The data 
base generated by the experiments must extend 
beyond thrust and energy balances, and should be 
capable of resolving detail at the level predicted by 
the computations (i.e., capable of measuring 
velocity components, translational, rotational or 
vibrational temperatures, neutral species and 
electron number densities). 

This paper describes our continued effort to 
develop and validate computational models of 
arcjet performance by way of direct comparison to 
detailed measurements of flow properties. We have 
selected a nominally 1 kW class hydrogen arcjet 
thruster as a benchmark for this comparison. 
Unlike ammonia arcjet flows, the chemical 
composition is relatively simple, and the transport 
properties and chemistry are better known. 
Furthermore, hydrogen is a strong contender as a 
propellant for future arcjet demonstration 
missions. 1 

The numerical results presented here are based 
on a single-fluid model that has been developed at 
Olin-Rocket Research Company. 8,9 The 
experimental data are obtained from facilities in the 
High Temperature Gasdynamics Laboratory at 
Stanford University. Details of both the numerical 
model and experiments are available in other 
papers. 8-12 j n this paper, we shall only discuss the 
sensitivity of the model to parameters which 
influence some of the physical processes described 
above, such as its sensitivity to changes in electron- 
ion and atomic hydrogen recombination rate 
coefficients. In a previous paper, ^ we performed 
a detailed comparison of experimental 
measurements of flow properties taken at a nominal 
power (1.5 kW) to properties predicted by the 
numerical model. Our focus here is to extend this 
comparison over a range of specific power. 

A brief description of the single-fluid model is 
given in Section II. An overview of the 
experimental facility is described in Section HI. A 
detailed comparison between experimental results 
and model predictions is presented in Section IV. 
A discussion of this comparison and an analysis of 


the model's sensitivity to chemical rate coefficients 
is given in Section V. A summary of this study and 
conclusions are given in Section VI. 


II Single- Fluid Model 

The single-fluid model is based on a solution of 
the axisymmetric Navier - Stokes equations 
coupled to relations which govern the behavior of 
the magnetic and electric fields, as formulated 
within the magnetohydrodynamic (MHD) 
approximation. 14 The model tracks the relative 
mole fractions of the various species which 
comprise the mixture by solving individual species 
conservation equations that include source terms 
for the net production of species via chemical 
reactions. The model considers a mixture that 
includes atomic and molecular hydrogen (H, H 2 ), 
electrons (e), and protons (H + ). Table I 
summarizes the chemical reaction set and nominal 
values of the reaction rate coefficients (in modified 
Arrhenius form) used in the calculations. Reverse 
reactions are computed as required by using 
detailed balance and curve fits (in JANNAF 
format 1 5) to the high temperature thermodynamic 
properties derived from spectroscopic data. 16 

In the present study, diffusive transport is 
neglected. This simplification is expected to 
significantly affect the calculations of the electrical 
conductivity, a , in colder regions of the flow and 
is discussed in more detail later in the paper. 
Neglecting the diffusion of neutral species, such as 
the slip velocity between atomic hydrogen and 
molecular hydrogen may be a more reasonable 
approximation, especially in the low density nozzle. 
Recent velocity measurements of atomic hydrogen 
and seeded helium at the exit plane of a 1 kW arcjet 
verified that slip velocities between species in the 
flow may indeed be negligible. 1 ^ 

In the formulation, the global momentum 
includes the Lorentz force associated with self- 
generated magnetic fields. Body forces associated 
with excess free charge (space charge) are 
neglected within the MHD approximation. Quasi- 
neutrality is assumed throughout the flow. The 
global energy equation includes Joule heating. The 
required transport coefficients (thermal 
conductivity, viscosity, electrical conductivity) are 
computed from the species mole fractions, collision 
integrals and mixture rules as provided by Yos. 1 ^ 

The magnetic field. 5, is computed from the 
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magnetic induction equation in the steady-state 
approximation 

Vx5x5-Vx^7?(Vxij)j = 0 (1) 

In writing Equation (1), we have neglected the 
contributions associated with ion-slip and the Hall 
effect. In the above equation, u represents the 
mean mass velocity of the flow and q = \/n 0 a is 
a measure of the plasma resistivity. A detailed 
derivation of Equation (1) can be found in 
Reference 14. The current density is derived from 
the computed magnetic field using Ampere's law. 

The importance of the Hall effect is determined 
by the Hall parameter 14 (product of electron 
cyclotron frequency and mean collision time for 
electrons and ions and electrons and neutral 
species). The Hall parameter is less than unity for 
conditions anticipated in the arcjet. Ion-slip is 
generally unimportant in highly-ionized plasmas, 
and is only significant in weakly-ionized plasmas 
when the Hall parameter is much greater than unity. 

For boundary conditions, we assume a no slip 
condition for velocity and zero flux condition for 
species number densities at the electrode surfaces. 
For temperature, we couple the heat flux at the 
electrode surface to a model of heat transfer within 
the anode. Heat transfer to the cathode is 
neglected. The model systematically calculates the 
temperature distribution and flux condition on the 
internal and external surface of the anode. Only the 
nozzle of the complete anode structure is modeled. 
A specified heat flux or temperature distribution is 
required at the boundary where the anode/nozzle 
joins with the rear arcjet housing. In general, this 
boundary condition is found to have a small effect 
on the predicted arcjet performance. Constraining 
the current to a direction normal to the anode at the 
anode surface acts as a condition on the the 
gradient of the magnetic field on this boundary. At 
the exit plane, we assume that the magnetic field is 
zero. The flow is assumed to be axisymmetric, and 
so we solve only for the azimuthal component in 
the magnetic field. 

As mentioned above, we believe that the single 
fluid model is not capable of accurately predicting 
the near-electrode plasma properties and therefore 
the near-electrode conductivity. In the vicinity of 
the cathode or anode, property gradients are high 
and ambipolar diffusion will create an 
overpopulation of electrons and ions in this region. 
Furthermore, we expect an elevated electron 


temperature, since electron thermal conduction is 
very efficient and thermalization between electrons 
and neutral species occurs at a finite rate. Since the 
single-fluid model used here does not account for 
these effects, the model will underpredict the 
electron number density, and in tum, over predict 
the Joule heating rate (per unit volume), 
P Q =j J/o. In order to compensate for this, a 
minimum value of the near-electrode conductivity 
is imposed. This conductivity "floor", is 
found to have a dramatic effect on the computed 
arcjet behavior. In particular, we find that the value 
selected for the conductivity floor greatly 
influences the computed arc voltage. For this 
reason, a value of <7^,, is selected to provide 
reasonable agreement between the computed and 
measured discharge voltage over a wide range of 
power to mass flow rate ratio, P/m. Typical 
values of cr^,, correspond to electron number 
densities of approximately lO 1 ^ cm‘^. This value 
is approximately ten orders of magnitude lower 
than the electron number densities predicted in the 
core of the arc. The fact that the computed arcjet 
performance depends so critically on this value 
further emphasizes the need to include a self- 
consistent treatment of the plasma behavior in the 
near-electrode region of the arcjet. A proposed 
self-consistent treatment of the quasi-neutral, near- 
electrode region in arcjet thrusters has been recently 
described in the literature. 19,20 

The computed arcjet voltage reflects the 
resistive contribution to the plasma conductivity, as 
constrained within the quasi-neutral assumption. 
As a result, the single-fluid quasi-neutral model 
does not formally treat and therefore capture the 
structure of the non-neutral sheath. The space 
charge generated in the immediate vicinity of the 
electrodes (within a distance of the order of a 
Debye length) gives rise to an additional voltage 
drop. In the case of a collisionless anode sheath, 
this may represent an additional mechanism for 
heat transfer to the anode. We use a simple model 
to arrive at an estimate of the anode sheath 
potential , <5* . We shall assume that the sheath is 
collisionless (valid in low pressure regions of the 
nozzle, when the mean free path, X, is much 
greater than the Debye length, X D ), and that the 
anode fall is strongly positive relative to the 
potential of the distant plasma. Furthermore, we 
shall assume that the total current, J , is carried by 
the electrons within the space charge region, and 
that it is approximately equal to the one-way 
random electron flux at the sheath-plasma 
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boundary. Under these conditions, the voltage drop 
across the sheath is described by Child's law, ^ and 
is given here as 


the single-fluid and a extended two-fluid model in 
describing the performance of a 10 kW hydrogen 
arcjet thruster has been previously described.^ 


$ 3 / 4=1 
5 2 




( 2 ) 


Here, e 0 is the permitivity of free space, e is the 
magnitude of the elementary charge, m e is the 
electron mass, and d s is the sheath thickness, 
which for a collisionless sheath, can be as large as 
the electron mean free path, X e . For conditions 
near the arcjet constrictor, we estimate X e = 3 x 
105 m. If we assume that the current density to the 
anode is of order 10-* A/m^, equation (2) above 
gives a maximum value of = 13V. For the 
results presented here, the sheath voltage is taken to 
be constant over the entire anode surface and the 
total power , P = I <t> s , is assumed to be driven 
entirely into the anode. The power lost in the 
anode sheath is added to the power dissipated 
through Joule heating to arrive at the total electrical 
power invested into the device. 

We note that this simple treatment of the plasma 
sheath is invalid when the densities are sufficiendy 
high to render the sheath collisional. Under these 
conditions, a self-consistent multi-fluid treatment of 
a collisional sheath in a flowing plasma may be 
necessary. 

The governing single-fluid equations are 
integrated using a finite-volume algorithm 
originally developed for magnetohydrodynamic 
flows at Los Alamos National Laboratory. ^ * A 
discussion of the numerical method may be found 
in References 8 and 9. Its original capabilities 
have been extended by Olin-Rocket Research 
Company to include finite rate chemistry, high 
temperature transport properties, radiative losses, 
and thermal conduction. The computational grid 
includes the region within the arcjet plenum, throat, 
and nozzle. The domain extends to the nozzle exit 
plane. The model therefore accurately predicts the 
internal flowfield of the thruster, but is incapable of 
modeling expansion waves immediately 
downstream of the nozzle exit plane. 

Although there are some obvious limitations to 
this model described here, it offers a framework on 
which to construct multi-fluid models. ^ 
Furthermore, it provides an opportunity to guide the 
design of higher power thrusters. The usefulness of 


m. Experiment 

The arcjet thruster used in the experiment is a 1 
kW class radiatively-cooled laboratory-type 
thruster designed and built at NASA Lewis 
Research Center. The tungsten nozzle has a 0.635 
mm diameter throat and a conical (20° half angle) 
diverging section with an area ratio of 225 (9.53 
mm exit diameter). A more complete description 
of the arcjet design is available elsewhere. 2 
Baseline operating conditions are 13.7 mg/s of H 2 
propellant with a maximum power dissipation of 
1 .5 kW (voltage and current levels of 148 V and 
10.1 A respectively). Some measurements were 
also performed at slightly lower mass flow rates 
(13.1 mg/s). At these propellant mass flow rates, it 
is possible to achieve stable operation at powers as 
low as 650 W. The arcjet is operated in a 0.56 m 
diameter cylindrical stainless steel chamber 1.09 m 
long. The 0.35 Torr background pressure is 
maintained by two 1250 cfm blowers evacuating 
through 6 inch diameter pipe. 

The experimental data of plume velocity and 
translational temperature is obtained by laser- 
induced fluorescence (LIF) excitation of the Balmer 
-a transition of atomic hydrogen. In employing 
LIF, a laser probe is focused to a small diameter in 
the arcjet plume. Subsequent fluorescence from 
excited atomic hydrogen is collected and imaged 
onto a detector. Flow properties are determined 
from the way in which the fluorescence signal 
varies with the laser frequency (wavelength). This 
spectroscopic technique has several advantages. It 
is non-intrusive, and measurements can be made in 
the vicinity of the exit plane where solid probes 
would fail due to the high-enthalpy nature of the 
flow. In addition, the ability to focus the probe 
laser to small sizes allows for high spatial 
resolution. Having high spatial resolution 
(approximately 100 mm) is particularly important 
in low power arcjets where large variations in 
velocity and temperature occur across the small, 1 
cm diameter of the nozzle exit. 

Velocity is determined from the Doppler shift of 
the fluorescence excitation spectrum. The 
translational temperature is determined from the 
shape of the fluorescence excitation spectrum. In 
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contrast to velocity, the temperature is more 
difficult to infer due to competing mechanisms 
which act to broaden the measured profile as well 
as the complex nature of the overlapping line 
components. For conditions present in the arcjet 
plume, the Balmer-a transition is primarily Doppler 
broadened and has a Doppler width that is much 
larger than the laser bandwidth. As a result, the 
shape of the recorded fluorescence spectrum 
accurately depicts the excited state velocity 
distribution function. The excited state atomic 
hydrogen is assumed to be in translational 
equilibrium with the ground state species. Details 
of the diagnostics have been presented 
previously. ^ 1 1 

A schematic of the experimental setup is shown 
in Figure 1 . The fluorescence excitation laser is a 
CW ring dye laser (Spectra Physics Model 380) 
pumped by a 6 W argon-ion laser and operating on 
DCM dye. It produces several hundred milliwatts 
of power at 15,233 cm' 1 (656 nm) with a nominal 
bandwidth of 10 MHz. Part of the laser beam is 
passed through a 2 GHz Fabry-Perot interferometer 
to monitor mode structure and ensure single mode 
operation. Another portion of the beam is directed 
into a wavemeter (Burleigh WA-10) which uses a 
Michelson interferometer to measure laser 
frequency to within 0.01 cm' 1 . A third beamsplitter 
directs the laser to a photodiode detector which 
monitors laser power. 

A beamsplitter divides the laser into two beams 
of approximately equal intensity. One beam is for 
axial excitation and the other is for radial 
excitation. The beams are expanded and then 
focused through windows to a waist of 0. 1 mm in 
the arcjet plume. The radial excitation beam is 
passed vertically through the plume. This allows 
the fluorescence signal to be collected on a 
horizontal axis normal to both excitation beams. It 
is collected through an 83 mm diameter window 
with a lens operating at 175.7 (76 mm diameter, 280 
mm focal length) and imaged onto a R928 
photomultiplier tube (PMT) with a magnification of 
2.7. A notch filter with a nominal 10 A (23 cm -1 ) 
bandwidth is placed directly in front of the PMT to 
filter extraneous light other than that from the 
transition of interest. A 0.4 mm aperture is 
mounted directly in front of the filter at the focal 
plane of the imaging lens to act as a field stop and 
determine the measurement location in the plume 
as well as the axial and vertical spatial resolution. 
Phase sensitive detection is used to discriminate the 


fluorescence signal from the background emission 
which is 30 to 500 times greater. The lock-in 
amplifier (SRS SR530) is synchronized to a 
mechanical chopper at 3 kHz. The output from the 
lock-in is collected and stored in a laboratory 
computer along with the output from the 
wavemeter and laser power detector. 

At each spatial point investigated the laser is 
automatically scanned in frequency by varying the 
voltage on a (piezo-controlled) 75 GHz intracavity 
etaion. Radial and axial excitation scans are 
performed on separate scans by blocking one of the 
beams. The laser is able to scan over a range of up 
to 2 cm' 1 , enough to capture most of the 
fluorescence excitation line. For adequate signal to 
noise, scan times are a few minutes in duration. 

For a stationary (unshifted) reference, both the 
arcjet itself and a stationary hydrogen discharge 
were used. It is assumed that the flow is 
axisymmetric which implies that the radial velocity 
is zero at the centerline. The symmetry of both the 
axial and radial velocity profiles confirms this 
assumption. The wavemeter provides enough 
accuracy and repeatability to forego simultaneous 
reference scans. The stationary discharge was 
produced in a low pressure hydrogen cell mounted 
in a microwave cavity which acted as a source of H 
atoms with essentially zero average velocity. 

The velocity and temperature are determined by 
a least squares fit to the lineshape model using the 
Levenberg-Marquardt method.^ 2 The model 
consists of the five components associated with the 
fine structure of the probed transition, with 
specified relative intensities and separations. 11 
Doppler broadening is the only broadening 
mechanism considered. There are three parameters 
that are used to fit the data to the model: Gaussian 
halfwidth, line shift, and magnitude. The 
temperature is determined from the Gaussian 
halfwidth. The velocity is determined from the 
shift. The magnitude was necessary to fit the other 
parameters but can also be used to determine 
relative excited H atom densities. 

The effect of saturation was evaluated by 
comparing LIF line profiles for a range of laser- 
power levels. For the 0. 1 mm laser spot size, the 
fluorescence was found to be linear below 60 mW. 
Thus all data was taken with laser powers well 
below this level. 

Two typical fluorescence excitation scans are 
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shown in Figure 2. The scan on the left is taken 
from the plume centerline 0.4 mm from the exit 
plane with axial excitation where the measured 
temperature is 4500 K. The scan on the right is 
from the colder, 600 K reference cell. At 
temperatures below about 1000 K, the Balmer-a 
lineshape begins to look more like a doublet. This 
is due to the fine structure discussed briefly in the 
previous section and in more detail elsewhere. * ^ 
The solid lines correspond to the model's fit to the 
data. The fluorescence intensity is the ratio of the 
output signal from the lock-in amplifier to the 
signal from the photodiode which monitors laser 
power. Accounting for laser power in this way 
preserves the lineshape since the fluorescence is 
linear and not saturated. The power varies by up to 
10% as the laser scans in frequency. Both profiles 
were normalized by their areas to facilitate 
comparisons of the lineshapes between the hot and 
cold gases. These profiles consist of 70 to 200 data 
points each, with the data taken at one second 
intervals. Profiles from regions further from the 
centerline require longer scan times due to 
decreased signal-to-noise ratios. 


IV. Results 

As described in Section II above, a-priori 
knowledge of the current-power operating 
characteristics of the arcjet was necessary in order 
to establish a value of the conductivity "floor", 
CTjnjn . Figure 3 shows the measured and predicted 
variation in operating current versus operating 
power for a specific power of 90 MJ/kg and three 
values of the conductivity floor. For a single 
conductivity floor, we find that the predicted 
current-power characteristics are best modelled by 
a value of cr^,, = 4.07 mho/m. We note however, 
that use of a single value of <7^,, may be 
inappropriate, as it is evident in the figure that this 
value of a mn over predicts the arc current at low 
arc powers, and underpredicts the current at high 
arc powers, over the range of powers studied here. 
In general, use of a single value of <7^,, = 4.07 
moh/m can result in a 25-30% overprediction or 
underprediction of the arcjet operating voltage for a 
given specific power. In order to more accurately 
model the current-power characteristics, we use a 
variable cr,^,, that is quadratically dependent on 
the arcjet power: 

<7™. =-3.969 + 11. 529 P- 4.099 P 1 moh/m 


where P is the arcjet power in kilowatts. 

We also note that although a power dependent 
value of is employed in the computation for 
the range of specific powers studied here, a more 
appropriate model for the current-power 
characteristics would include a mass flow 
dependence on <7^,, (i.e., a dependence on specific 
power). The dependence on specific power is not 
included in this study, and we assume that over the 
range of specific powers investigated here, 
adjustments in the value of <7^,, are insignificant. 

In the simulations, a value of <1> 3 =13V is taken 
for the anode sheath voltage. This value is found 
to accurately reproduce the performance 
characteristics of a lOkW hydrogen arcjet^ and is 
approximately that estimated by Child’s Law when 
the collisionless sheath thickness approaches the 
electron mean free path. The sensitivity of the 
computed arcjet performance to values of the 
sheath voltage drop is discussed in Section V. 

For the calculations presented here, we have 
assumed an adiabatic boundary condition for the 
surface where the nozzle assembly is attached to 
the remaining arcjet structure. Although this 
assignment is somewhat arbitrary, we find that on 
the basis of a comparison to calculations that 
employ an isothermal condition of 1000K, the 
computed flow properties are only weakly sensitive 
to this boundary condition. 

Unless otherwise stated, the experimental 
measurements presented here are for hydrogen 
mass flow rates of 13.7 mg/s. This is 
approximately 7% higher than the mass flow rate 
used for earlier comparisons to model 
predictions. * 3 Figures 4-6 compare the computed 
and measured exit plane velocities and kinetic 
temperature for values of specific power, 
PI m =64, 73, and 1 10 MJ/kg. In all cases, there 
seems to be good agreement in magnitude and in 
profile shape (near the jet core), between the model 
predictions and experimental measurements of exit 
plane axial velocity (see Figure 4). The systematic 
differences between the model predictions and 
experimental results in the regions close to the 
nozzle wall may be indicative of additional post- 
nozzle flow expansion that is not simulated by the 
model. We note that the measurements are made 
with a finite spatial resolution (0.1 mm) 
approximately 0.4 mm from the true exit plane and 
may capture the gasdynamic behavior in the region 
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of the nozzle lip. The differences however, are 
within the experimental uncertainties. The 
increased uncertainty in the measured velocity as 
one moves outwards from the central arcjet axis 
arises from the concomitant drop in signal as a 
result of the dramatic drop in excited state number 
density. 

As can be seen from Figure 5, there is 
reasonable agreement between the simulated and 
measured spatial variation in the radial velocity 
component for the range of specific powers studied. 
Positive and negative values in Figure 5 represent 
data taken at positions extending across the entire 
arcjet plume. Since the laser is directed downwards 
with respect to the laboratory reference frame, a 
negative velocity implies downward motion and is 
not to be confused with motion of the flow along a 
direction towards the arcjet axis. The measured 
radial velocities are found to vary nearly linearly 
with increasing distance from the arcjet axis. The 
measured velocities are in satisfactory agreement 
with model predictions near the arcjet exit plane 
core. However, we find that there is a noticeable 
difference between model predictions and 
experimental measurements made at positions far 
from the plume axis. We emphasize again 
however, that the model forces the radial velocity 
component to zero for the internal flowfield 
calculation in order to invoke a no-slip boundary 
condition. In fact, additional gasdynamic effects 
near the lip of the nozzle may be captured in the 
measurements of the radial velocities. It is 
noteworthy that both the measurements and 
calculations indicate that the radial velocities are 
only weak functions of the arcjet specific power. 
Also, the magnitude of the radial velocity reaches 
approximately 4 km/s for the highest specific 
power studied. This occurs near the nozzle lip, 
where the measured axial velocity is of the same 
magnitude thus implying a 45° flow direction at this 
spatial location. 

As can be seen from Figure 6, a significant 
difference between model predictions and 
experimental measurements of exit plane 
temperature is evident, both in magnitude and in 
profile shape. In particular, we find that the 
discrepancy increases with increasing specific 
power. 

These results are consistent with earlier 
comparisons^ of the measured and predicted 
centerline axial velocity and temperature variation 
with specific power. These earlier data are 


presented in Figure 7. As can be seen in Figure 7, 
good agreement between simulations and 
experiment was obtained for centerline values of 
temperature at relatively low specific powers (50 
MJ/kg), however, rather poor agreement (factor of 
two difference) was obtained at the higher values of 
specific power (110 MJ/kg). Both the 
measurements and simulation clearly indicate that 
when the input power is more than doubled, the 
centerline velocity increases by only 30%. This is 
consistent with the observed variation in thrust with 
changes in operating power^ for medium power 
hydrogen arcjet thrusters. The experimental results 
also indicate that the kinetic temperature increases 
by a factor of three over the range of powers 
studied. At higher powers, there is an increased 
fraction of the enthalpy convected out the nozzle. 
The increased specific impulse therefore comes at 
the expense of thrust efficiency (defined here as 
g 2 l]pPI 4m). Figure 8 shows the computed 
variation in both the specific impulse and efficiency 
with arcjet power at a constant mass flow rate of 
13.1 mg/s. Although the computed thrust 
efficiency is seen to decrease with increasing 
specific power, the experimental measurements 
suggest that the efficiency must fall more 
precipitously to account for the considerable rise in 
exit plane temperature. 

The measurements of temperature are based on 
the measurements of Doppler broadened 
fluorescence excitation profiles and are found to be 
very sensitive to the baseline fitting procedure. 1 ^ 
The experimental uncertainty is greater for 
temperature than it is for velocities, and is reflected 
in the scatter in the experimental data. Despite the 
increased uncertainty, we find that the discrepancy 
between measurements and model predictions lies 
outside of the expected experimental uncertainty in 
temperature for the higher specific power cases. 
For the highest specific power investigated here, we 
find that the model underpredicts the temperature 
by at least 50% over a good portion of the arcjet 
radial domain (r < 2 mm). The reason for this 
difference is speculated on in Section V. There has 
been some question raised^ as to wether the 
excited electronic state that is probed in these 
experiments is produced with excess kinetic energy 
resulting from dissociative recombination. This is 
certainly possible, although it should be noted that 
for the highest specific power case, the model 
significantly overpredicts the specific impulse 
(about 1000s), which, for conditions close to that 
of the high specific impulse case, has been 
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measured to be approximately 830 s.^ An 
overprediction of 20% in the I S p is consistent with 
the observed difference in the measured and 
computed temperatures. 

Although measurements of exit plane electron 
number density and ground state atomic hydrogen 
number densities under the above mentioned 
conditions have not been made, it is noteworthy 
that the predicted centerline values, which can be 
inferred from the predicted radial variation as 
shown in Figure 9, are consistent with 
measurements previously reported in the literature. 
For a specific power of 1 13 MJ/kg and a mass flow 
rate of 13.1 mg/s, a centerline exit plane electron 
number density of lO 7 ^ m'^ falls within the 
centerline value of 10^ m' 7 reported by 
Manzella 7 ^ at a downstream location of 1 cm, and 
a centerline value of 10 7 * m"’ reported by Zube, 7 ^ 
1 cm upstream of the exit plane inside the 
expansion nozzle. More recent measurements 7 ^ of 
the electron number density in the arcjet throat 
indicate that the axial electron number densities are 
2-5 times higher than that predicted for the highest 
specific powers studied. Centerline values of the 
ground state atomic hydrogen 2 cm from the exit 
plane as high as 5 x 1<)18 m -3 have been 
reported, ^ 7 indicating that additional plume 
expansion or hydrogen recombination may be 
taking place within a very short distance of the exit 
plane. 


V. Discussion 

Although the single fluid model predictions and 
measurements of arcjet velocities compare 
favorably, the differences observed between 
predictions and measurements of exit plane 
temperatures at higher specific powers (> 75 
MJ/kg) suggest that some features of the model 
require further refinement. It is apparent that there 
is a systematic underprediction of the exit plane 
temperature. This result is consistent with the 
overprediction of specific impulse, based on 
measurements made on. similar thrusters at NASA 
Lewis. We find that these differences cannot be 
accounted for by refinements in the anode (internal) 
heat transfer model. Instead, we believe that there 
are physical processes that are either not included 
in the model (i.e., ambipolar diffusion), or are 
included but not accurately represented (i.e., 
ionization), which must account for the observed 
differences in predicted and measured arcjet 


performance. It is apparent that there is an 
underprediction in the dissipated power, and 
therefore an underprediction in the exit plane 
temperature. Part of the discrepancy may be 
accounted for by the simplified treatment of the 
anode sheath. We used here an upper limit to the 
sheath potential, <J> Jt of 13V, which represents a 
10% loss of power to the arcjet anode. According to 
the experimental data, increasing the dissipated 
power by this amount (i.e., assuming no electrical 
power loss to the anode) can lead to as much as a 
1000 K increase in the exit plane centerline 
temperature. It is important to note however, that it 
the model predictions that does not reflect such a 
trend in that a 10% increase in power raises the 
temperature by a mere 300 - 400 K. 

The electron third-body recombination rate used 
here is representative of that estimated for a 
strongly recombining low temperature (i.e., 3000- 
4000 K) plasma 7 ®- The conditions in the 
constrictor of the arcjet are that of a high 
temperature (> 30000 K) ionizing plasma. A 
proper treatment of the ionization/recombination 
processes would require the coupling of a 
collisional-radiative model to the flow calculations. 
At this time, such an addition to the model would 
be prohibitively expensive from a computational 
standpoint. If the rate presently employed is too 
low, then the ionization rate calculated from 
detailed balance will also be low. This will result 
in too low a constrictor electron number density (as 
is observed 7 ^), and hence too low a conductivity. 
A low conductivity will give rise to a relatively 
high voltage drop over a shorter distance (for the 
same current density) and hence a shorter arc 
length. In order to increase the flow temperature 
(decrease the thrust efficiency) we must therefore 
increase the recombination rate that is employed in 
the simulation so that the arc penetrates further into 
the supersonic region of the expanding flow. 
Simulations were performed for values of 10 and 
100 times that of the nominal recombination rate 
coefficient that is used, for the highest arcjet 
specific power studied. The results of this 
computational study indicate that changes of such 
magnitude have only a marginal influence on the 
exit plane temperature. More dramatic differences 
between the true electron-ion recombination rate 
coefficient and that employed here would be 
necessary to account for the differences in the 
measured and predicted exit temperature. 

Perhaps a more important process not accounted 
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for in the model is species diffusion. In the single- 
fluid model, this can be accounted for to first order 
by the addition of diffusive fluxes for both neutral 
and charged species (ambipolar diffusion). 
Because of the significant temperature gradients 
that can exist in the arcjet constrictor (10° K/m), it 
is necessary to include thermal diffusion. A more 
rigorous treatment of the decoupled electron, ion, 
and neutral species dynamics would require the 
addition of an appropriate momentum equation for 
each plasma constituent. Such multi-fluid models 
would permit the diffusion of electrons and ions 
into colder regions of the constrictor and nozzle, 
thereby elevating the local plasma conductivity, 
and removing the need to artificially impose a 
conductivity "floor" in the region near the nozzle 
wall. 

The extension of the single fluid model to two 
fluids has been previously described. 9 A two fluid 
model, with the electrons distinguished from both 
the neutrals and ions, would allow for the possible 
departure from local thermal equilibrium. In the 
two fluid model, energy is added to the electron 
pool through Joule heating, and is transferred to the 
heavier species by both elastic and nonelastic 
processes (ionization, dissociation). An elevated 
electron temperature is established, especially near 
the anode, due to the more efficient conduction of 
electron energy across regions of the plasma where 
steep temperature gradients may exist. This two- 
temperature nonequilibrium would also help in 
maintaining an elevated plasma conductivity, and 
would therefore also result in an increase in arc 
length and therefore decrease in nozzle 
performance or thrust efficiency. 

In general, we find that the differences in 
measured and predicted exit plane temperature 
cannot be accounted for by uncertainties in atomic 
hydrogen recombination rate coefficients. We find 
that a large (i.e., order of magnitude) change in the 
rate coefficients for third-body enhanced 
recombination of atomic hydrogen (H + H + M -> 
H 2 + M) has only a minor influence on the 
calculated exit-plane temperature field and arcjet 
performance. This indicates that frozen flow losses 
associated with finite-rate molecular recombinauon 
processes are difficult to recover under present 
arcjet design constraints. 

At present, the computational domain is limited 
to the exit plane of the arcjet. We impose a zero 
magnetic field constraint on this boundary. This 
boundary condition is expected to influence the 


current distribution within the nozzle. It is well 
known that the arc attachment points extend well 
beyond the constrictor and into the expanding 
portion of the anode/nozzle in these low power 
arcjet thrusters.^^ Constraining the magnetic field 
to be zero at the exit plane may limit the extent to 
which the arc penetrates out from the constrictor 
and therefore constrain the ohmic dissipation to the 
subsonic portion of the flow thereby increasing the 
nozzle performance. 


VI Summary and Conclusions 

In this paper, we have described a detailed 
comparison of experimentally measured low power 
hydrogen arcjet plume properties to model 
predictions. The experiments focus on the 
measurement of the exit plane vector velocity field 
and kinetic temperature by laser induced 
fluorescence. The model is based on a single-fluid 
magnetohydrodynamic description of the flow. 
Excellent agreement between model predictions 
and experimental exit plane velocity profiles is 
obtained over a wide range of specific power, 
despite the complex nature of the flow. At higher 
powers however, there is significant difference 
between the measured and predicted exit plane 
temperature. The model overpredicts the 
temperature by as much as 2000K. These 
differences cannot be attributed to the present 
treatment of the collisionless anode sheath or finite- 
rate molecular recombination (dominant frozen 
flow loss) and is believed to be more closely related 
to the treatment of other physical processes such as 
finite-rate ionization, ambipolar diffusion or 
possibly the selection of exit plane boundary 
conditions. 

The simulations overpredict performance (Isp) 
at the higher specific powers studied. This 
overprediction in performance is consistent with the 
underprediction in temperature and implies an 
overprediction in the thrust efficiency or nozzle 
performance. The thrust efficiency is shown to be 
only marginally sensitive to order of magnitude 
changes in the value of the electron third-body 
recombination rate coefficient that is employed in 
the simulation. However, little data is available as 
to what value of this rate coefficient is appropriate. 
This suggests that further refinements in the 
treatment of finite-rate ionization and species 
diffusion in these thrusters may be necessary in 
order to extend the predictive capabilities of the 
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model. 
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Reaction 

A 

B 

E 

H + H + M^H 2 + M 

6.4xl0 17 

-1 

0 

H + +e + M H + M 

5.26x1 0 26 

-2.5 

0 

H 2 + e — > H + H + e 

1.91xlO n 

1 

203,000 

H + + e + e — > H + e 

7.08X10 39 

-4.5 

0 


Table 1. Chemical reaction mechanisms. 



Figure 1 . Schematic of experimental facility. 
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Figure 2. Typical fluorescence excitation scans 
taken with the laser focused along the 
centerline 0.4 mm from the exit plane. 
The more narrow profile with resolved 
features is taken from a 600 K low 
pressure stationary discharge ceil. 



Arc power [kW] 

Figure 3. Comparison of experimental cuurent- 
power characteristics at a specific 
power of 90 MJ/kg to that predicted 
for three different values of the 
conductivity floor, a mn . 



Radial Position [mm] 

Figure 4. Comparison between measured 
(symbols) and predicted (lines) 
spatial variation in exit plane axial 
velocity, m = 13.7 mg/s. 



Radial Position [mm] 

Figure 5. Comparison between measured 
(symbols) and predicted (lines) exit 
plane radial velocity component, m = 
13.7 mg/s. 
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Figure 6. Comparison between measured 
(symbols) and predicted (lines) 
variation in the exit plane temperature, 
m = 13.7 mg/s. 
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Figure 7. Computed (solid line) and measured 
(open symbols) variation in centerline 
axial velocity and temperature with 
power, m = 13.1 mg/s. 



efficiency and specific impulse with 


power, m = 13.1 mg/s. 



Figure 9. Computed radial variation in the exit 
plane electron and H-atom number 
densities, m = 13.1 mg/s. P/m =113 
MJ/Kg. 
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Low-density flow of molecular hydrogen from a small nozzle is studied using numerical 
and experimental techniques. The conditions in the nozzle indicate that non- 
equilibrium effects will significantly influence the flow. Therefore, the numerical 
analysis is undertaken using a Monte Carlo approach. The experimental studies 
employ spontaneous Raman scattering. Comparisons of the measured data and 
computed results are made for total number density, rotational temperature, and for 
the number density of the first rotational level. The numerical results are found to be 
quite sensitive to the rotational relaxation rate, and a strong degree of thermal non- 
equilibrium is observed at the exit plane of the nozzle. Comparisons between 
experiment and analysis permit estimation of the rotational relaxation rate for 
hydrogen. Investigations are also conducted for expansion of the supersonic jet into a 
finite back pressure. The interaction of the plume with the chamber background gas is 
found to form shock waves in both the simulations and experiments. This phenomenon 
is investigated further by increasing the background pressure. Direct comparison of the 
simulation results and experimental measurements is very favourable. 


1. Introduction 

Low-power hydrogen arcjet propulsion engines are under development for use on 
spacecraft for orbit transfer and station-keeping manoeuvres. Potentially, these 
engines offer significant improvements in specific impulse over current chemical 
rockets. However, the thrust efficiencies achieved to date with these rockets have only 
reached about 30%. Many of the factors that cause this relatively low efficiency are 
concerned with the fluid mechanics and chemical kinetics of hydrogen. The purpose of 
the present study is to examine some of these effects at a fundamental level. Detailed 
studies are performed for flow of hydrogen in a small nozzle. The geometric size and 
mass flow are representative of a spacecraft thruster. These aspects combine to give a 
low Reynolds number in the nozzle throat that yields flow in strong thermodynamic 
non-equilibrium. In this paper, both numerical and experimental studies are reported 
for an unheated hydrogen flow. These represent the first phase of a continuing effort 
to investigate arc-heated expanding hydrogen flows. The experimental investigation 
employs Raman scattering to obtain number density and rotational temperature. 

Sophisticated experimental techniques are being applied to investigate flow 
properties in the nozzle and plume of various arcjet thrusters. In these investigations, 
the thruster plume expands into a vacuum tank. Of course, owing to the finite pumping 
speeds of ground-based experimental facilities, the back pressures obtained are orders 
of magnitude higher than orbital flight conditions. There is an interest in understanding 
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how measurements taken in such facilities may be affected by the interaction of the 
expanding thruster plume with this residual pressure in the test chamber. This study 
addresses the issue for the hydrogen thruster operated in the arcjet facility at Stanford 
University. 

Unfortunately, it is quite difficult to estimate such effects. The plume expanding from 
a typical arcjet is complicated by viscous and chemical kinetic effects. The flow field is 
two-dimensional such that densities at the nozzle exit centreline may be an order of 
magnitude higher than those at the nozzle lip. In addition, the physical behaviour that 
occurs in the interaction region between the plume and the residual gas is not simple. 
Depending on flow conditions, shock waves may be formed, and the viscous and 
kinetic effects again become significant. Therefore, a detailed study of the interaction 
phenomenon has been undertaken. Experimental investigations again employ Raman 
scattering to measure number density along the axis of the plume. The measurements 
are carried through the interaction region where the thruster plume encounters the 
residual gas. In addition, Monte Carlo calculations of the same flow conditions are 
performed. 

Details of the experimental and numerical methods are provided in the following 
sections. In the numerical study, particular attention is given to the modelling of 
rotational relaxation for H,. A description of the nozzle and plume flows is given, and 
comparison of the numerical results and experimental data is discussed. 


2. Experimental investigations 

This section describes experiments performed using spontaneous Raman scattering 
to determine the rotational-level populations of H, for the flow in the expansion plume 
of a small rocket. The theory of this technique is described first. Discussions then 
follow on the experiments performed, and the methods employed for interpretation of 
the results. 

2.1. Theory 

Spontaneous Raman scattering is an inelastic, linear, two-photon scattering process. In 
this case, an incident photon scatters off an H 2 molecule causing a change in the 
molecular quantum state and a corresponding change in the energy of the scattered 
photon. The intensity of the signal is given by Valentini (1987): 


N u = N L ln^(%) di2. 


\*0)u 


( 1 ) 


where N i f is the number of Raman photons collected for transition from initial state / 
to final state /; N L is the number of incident laser photons; / is the sample length 
imaged; « 0 is the number density of scattering species; is the population 

fraction in state i; v i is the initial vibrational state; J t is the initial rotational state; 
(dcr/dP),. f is the differential Raman cross-section for transition from state i to state /; 
and c IQ is the collection solid angle. 

The number of scattered photons is linearly proportional to the number of incident 
photons and the density of the scattering species. Spontaneous Raman scattering has 
the advantage of allowing spatially resolved measurements of ground-state densities 
with an easily interpretable signal. The primary disadvantage is the small scattering 
cross-section (typically 10~ 33 m a ). Measurements here are based on the purely 
vibrational Raman transitions from v = 0 to i> = 1 with AJ = 0 for J = 0. 1,2 and 3 
(i.e. the Q-branch). Table 1 shows a summary of initial and final states, wavelength. 


Low-density supersonic jets of hydrogen 


43 



*7 

J, 

J, 

A 

(nm) 

<r 

10- M nr 


Relative 

intensity 

0 

1 

0 

0 

612.325 

1.544 

0.128 

1.982 

0 

1 

1 

1 

612.103 

1.648 

0.657 

10.828 

0 

1 

2 

2 

611.662 

1.645 

0.117 

1.930 

0 

1 

3 

3 

611.004 

1.665 

0.092 

1.529 


Table 1. Summary of Raman transitions 


cross-section. Boltzmann fraction, and relative intensity for these transitions using a 
laser wavelength of 487.986 nm and a gas temperature of 300 K. 

The density can be determined by using the absolute intensity from all transitions 
combined in conjunction with the signal from a reference cell. The sum of the intensity 
of each transition divided by its cross-section is proportional to the density and a 
constant determined by the optical system : 


y 

f (dtr/dO),., 




y 

i 


= N L ln 0 dQ 


= 2 « 0 . 


( 2 ) 


A spectrum recorded from a known gas density can be used to determine a and hence 
the density from subsequent spectra. 

The relative intensity of the various transitions can be used to determine a rotational 
temperature, assuming partial thermodynamic equilbrium exists among the rotational 
states. The intensity of each transition divided by its cross-section is proportional to the 
density and the population fraction of the initial state: 




(dcr/dQ) if 


n n 


= 


«o 


= zn 0 f[T,i). 


( 3 ) 


When two or more transitions are observed, their intensities can be fitted to a function 
of temperature and density (using the predetermined optical system constant). 


2.2. Experiments 

The nozzle considered in the present study is based on a 1 kW class NASA hydrogen 
thruster. The converging-diverging nozzle geometry is shown in figure 1. Note the 
cathode that is present in the device. The nozzle has a 0.64 mm throat diameter with 
a 20° half-angle, and an exit-to-throat area ratio of 225:1. The conditions chosen for 
examination represent a cold-flow configuration in which the arc is not ignited. The 
specific flow conditions are a stagnation temperature of 295 K and a mass flow rate of 
14.1 mgs -1 . These conditions correspond to a Reynolds number at the nozzle throat 
of about 2800. This value indicates that the flow is dominated by the effects of thermal 
relaxation, and by viscosity in the nozzle boundary layer. The thruster is operated in 
a vacuum chamber 1.09 m long and 0.56 m in diameter with optical access through 
75 mm diameter ports. It is mounted on a two-axis translation stage (axial and radial). 
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Figure I. Nozzle geometric configuration (dimensions in mm). 



Figure 2. Experimental apparatus. 


The tank is evacuated by two Roots blowers backed by mechanical roughing pumps 
with a total capacity of 1.2 m 3 s -1 . This system is capable of maintaining a back 
pressure of 43 ±2 Pa with a flow rate of 14.1 ±0.1 mg s _1 . 

The Raman excitation source is an Ar + laser with 2.5 W in the 487.986 nm line used 
for this experiment. The laser beam is focused axially into the plume with a beam waist 
under 100 pm (see figure 2). Scattered light is collected and focused onto a 
monochromator entrance slit with an achromatic two-lens system. The collection 
optical axis is perpendicular to the laser polarization and the direction of laser 
propagation. The spatial resolution is 4 mm in the axial direction (determined by the 
slit height) and 0. 1 mm in the radial direction (determined by the laser beam waist). 
The monochromator is a 1.0 m focal length f/8.7 single-pass Czemy-Tumer design 
with a 1 180 line mm -1 grating blazed at 500 nm. A long-pass filter is placed before the 
entrance slit to further reduce the intensity of collected laser light to negligible levels. 
The monochromator slits are typically 300 pm giving a spectral resolution of 0.2 nm. 

Light is detected at the exit slits by a cooled low-noise photomultiplier tube (PMT) 
with dark noise of 0.4 counts s" 1 . The PMT output is amplified and sent to a 
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Figure 3. Schematic of arcjet in vacuum chamber (dimensions in m). 


continuously gated photon counter. The monochromator is scanned over the spectral 
ration of interest and the photon counter output is recorded on a computer. The signal 
intensity is calibrated by filling the chamber with 10 kPa of H 2 . A spectral scan is 
recorded with the operating conditions unchanged from those used for the data scans. 

The overall experimental configuration is shown schematically in figure 3. Of 
particular interest is the position of the vacuum pumps with respect to the expanding 
flow. The location of the pumps will clearly give a non-symmetric background pressure 
field. It is presumed that such effects are small given the relatively large distance of 
separation between the pumps and the arcjet. 

2.3. Data analysis 

Spectral scans of the Raman Q-branch are recorded at various radial positions from 
the centreline to the nozzle edge. The collection segment starts at the nozzle exit plane 
and extends 4 mm downstream. Scans are also made on the centreline at various axial 
positions. The total number of scans taken is limited by the somewhat long scan 
duration of 40 minutes at 0.05 nm min -1 . In order to provide a more detailed view of 
the flow profile in the exit plane, data are also collected at a fixed spectral position while 
scanning the axial position. 

The spectral data are analysed by fitting to theoretical spectra with the unknown 
quantities such as temperature and density as fitting parameters. The theoretical 
spectra are obtained by calculating the predicted position and intensity of the various 
Q-branch lines, then convoluting with the instrument response function. The line 
positions are taken from the experimental values of Veirs & Rosenblatt (1987) which 
are in excellent agreement with the quantum mechanical calculations of Wolniewicz 
(1983). The line intensities are proportional to the Boltzmann fraction in the initial 
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state and the transition Raman cross-section. The cross-sections are taken from the 
dynamic polarizability calculations of Ford & Brown (1973). Since these values are 
averaged over all solid angles for unpolarized light, it is assumed that any variations 
with solid angle or polarization are the same for all Q-branch transitions. Note that the 
Ford & Brown cross-sections are defined in terms of scattered energy rather than 
scattered photons, as defined in (1). The two definitions differ by a factor of V'V 
where A s is the scattered wavelength and A 0 is the incident wavelength. 

As expected, the observed spectral features are principally instrument broadened. 
The Raman transitions are primarily Doppler broadened at the low pressures involved. 
According to Weber (1973), the Doppler width (FWHM) is 

2/ m i/2 r e t /2 

Av D = -[2\n2—J j^4(vj + »’ 0 v R ) sin 2 -+ , (4) 

where T is the gas temperature = 300 K; m is the hydrogen molecule mass; i< 0 is the 
incident photon frequency; v R is the Raman frequency shift; and 6 is the collection 
angle = 90°. 

Using (4), a value of 3.57 GHz is obtained for Av D . This Doppler width of 4.5 pm and 
the laser linewidth of 5 pm are both small in comparison to the 0.2 nm monochromator 
resolution. Hence the monochromator slit function is used to represent the instrument 
function for fitting purposes. Since the monochromator slits are equal in width and the 
resolution is large compared to the diffraction limit (0.01 nm), a triangular function is 
utilized. 

Two fits are performed. The first fit, NTE (Non-Thermodynamic Equilibrium), does 
not assume that the flow is in thermodynamic equilibrium. The rotational population 
fractions (and thereby relative intensity of the transitions) are freely varied to best fit 
the data. Additional fit parameters include a wavelength offset to account for 
monochromator error, an intensity offset to allow for PMT dark noise, and the width 
of the slit function. 

The second fit (LTE) assumes that the flow is in local thermodynamic equilibrium 
(LTE) with a rotational temperature dictating the rotational population fractions. The 
two parameters in this fit are the temperature and number density. The wavelength 
offset, intensity offset and spectral resolution have the same values determined by the 
NTE fit. The rotational energy levels used to compute the Boltzmann fractions are 
calculated with a six-parameter fit from Jennings, Rahn & Owyoung (1985). This fit is 
in excellent agreement with the energy level calculations for J=0, 1,...,5 of 
Wolniewicz (1983). The vibrational partition function is calculated assuming a 
harmonic oscillator using the energy difference from v = 0 to v — 1 also provided by 
Wolniewicz. The harmonic approximation introduces negligible error in absolute 
intensity since the population fraction in the ground state is close to unity at the low 
temperatures of this study. 

At the low temperatures involved here, the symmetric/antisymmetric nuclear spin 
states of H., do not reach equilibrium. This equilibration process can be extremely slow 
in the absence of a catalyst (McQuarrie 1973). Owing to the coupling of nuclear and 
rotational spin states, this non-equilibrium results in the rotational populations 
deviating from a Boltzmann distribution. Since H nuclei have spin h (i.e they are 
fermions), H., must have an antisymmetric total wave function: 

w — a/ m w w \p i 

total trans rot 1 rib elec 1 nuclear' 

The only components of the total wave function which can be antisymmetric are the 
nuclear and rotational wave functions (see table 2). Therefore, an H., molecule must 


Wave function Symmetry Degeneracy 

Translational Symmetric 1 

Rotational 

J even Symmetric 2J + 1 

/odd Antisymmetric 2J+1 

Vibrational Symmetric 1 

Electronic: '£+ Symmetric 1 

Nuclear 

Ortho Symmetric 3 

Para Antisymmetric 1 


Table 2. H 2 wave function 

be in a symmetric nuclear state and antisymmetric rotational state or vice versa. These 
two combinations are often referred to as ortho and para respectively. The nuclear 
degeneracy ratio of 3 : 1 (ortho-para) results in a 3 : 1 statistical weighting of odd-even 
J states. Above 300 K, the upper rotational states are sufficiently populated in a 
Boltzmann distribution that the ortho-para ratio is in fact 3:1. However, at lower 
temperatures where the population increasingly fills J = 0, the ortho-para ratio 
decreases. In the limit as T approaches 0, the equilibrium ortho-para ratio also 
approaches zero. 

The implication of this phenomenon is that equilibrium dictates an ortho-para ratio 
below 3 : 1 at the temperatures in this flow. However, since the ortho-para ratio is 
frozen at the 3 : 1 high-temperature limit, the rotational population fractions will not 
be Boltzmann, even if the rotational states are in LTE. This effect is accounted for in 
the fitting process by upholding the 3 : 1 ortho-para ratio, while at the same time 
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Figure 5. Boltzmann plot of rotational population under reference conditions. T — 323 ± 1 K, 

n = (2.27±0.02) x 10 24 m' 3 . 

requiring that the odd and even states independently satisfy a Boltzmann distribution 
with the same rotational temperature. On a Boltzmann plot, the odd and even states 
each fall on a line with the same slope (i.e. same temperature) but a different intercept. 
As the temperature is decreased the intercepts move further apart. 

In order to measure absolute density, reference scans are taken with the tank 
containing 10 kPa of hydrogen. A sample reference scan showing the raw data, filtered 
data, NTE fit, and LTE fit is shown in figure 4. The raw data are low-pass filtered to 
allow the spectral features to be seen more clearly. Fitting to the raw data and filtered 
data is found to give identical results. The two fits agree quite well, as expected, since 
the gas should be in equilibrium. A Boltzmann plot of the results from the two fits 
shows that the data match the equilibrium fit to within experimental error (figure 5). 

A scan taken in the exit plane of nozzle is shown in figure 6. The intensity is lower 
than the reference intensity and the shot noise is much more evident. It is also clear that 
the population has shifted almost entirely into the 7 = 0 and 7 = 1 states, indicating a 
lower temperature. The LTE fit temperature was 1 13 ±5 K. From the Boltzmann plot 
shown in figure 7, one can see the equilibrium fit lines match the 7 = 0, 1, and 2 
intensities quite well. The discrepancy between the equilibrium and non-equilibrium fits 
at = 3 is not surprising since the signal-to-noise ratio of that data point is very poor 
(< 1 ). 

The calculation of number density is slightly more accurate than temperature, 
especially at the lower signal levels. The cross-sections of interest are almost equal, so 
it is not critical that the fitting process determine exactly what portion of the 
overlapping peaks corresponds to 7 = 0 and 7=1. However, the error in density of the 
data scan is the sum of the density parameter fitting error and the reference density 
error. 
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A more detailed view of the exit-plane flow profile is obtained by setting the 
monochromator at the J = 1 peak and scanning the flow position radially. For this 
step, spectral resolution is less important, so the slits are opened to 400 pm to maximize 
the signal. This signal is proportional to the population in the J = 1 state, so it is 
normalized by the J = 1 signal from the spectral scans to give the density of state 
7=1. 

For the experimental data taken along the plume centreline, again the signal from 
the 7 = 1 state is measured. In this study, data are taken along the arcjet axis from the 
exit plane to 60 mm downstream. The back pressure in the vacuum tank is varied from 
43 to 67 Pa by varying the number of pumps used. Data are also taken with the 
chamber evacuated below 1 Pa to give a background signal, and filled with 10 kPa 
hydrogen for a reference signal. The plume density is found by subtracting the 
background signal from the measured signal, then scaling by the known density of the 
reference scan. The error in density is estimated as one standard deviation, assuming 
that the counts follow a Poisson distribution. The errors range from 6% at the highest 
densities to 36 % at the lowest. 


3. Numerical investigations 

This section discusses the numerical and physical aspects of simulating the flows of 
hydrogen using a Monte Carlo method. The general methodology for simulating the 
nozzle flow is described, details of the kinetics models employed are provided, and the 
rocket plume simulations are described. 

3.1. Nozzle flow simulation 

The low Reynolds number at the nozzle throat indicates that thermal non-equilibrium 
effects may be significant in the nozzle flow. The local value of the Reynolds number 
at the nozzle exit will be reduced by a factor close to the ratio of the exit and throat 
areas. Under such conditions, it is well established that the Navier-Stokes equations 
are invalid, owing to the failure of the linear forms for viscosity and thermal 
conductivity. A previous investigation of a small nitrogen thruster with a throat 
Reynolds number of 850 was reported by Boyd et al. (1992). In that investigation, 
experimental measurements of Pitot pressure were compared to continuum calculations 
at the nozzle exit. It was found that the continuum method overpredicted the pressure 
by 25 %. An alternative numerical methodology for low-Reynolds-number flows is the 
direct simulation Monte Carlo method (DSMC) developed by Bird (1976). Instead of 
solving an equation set, this method attempts to provide a direct simulation of the flow 
by employing a large number of model particles. Boyd et al. (1992) employed this 
method and found it to give excellent agreement with the experimental data. Further 
studies by Penko et al. (1993) and Boyd, Jafry & Vanden Benkel (1994) have shown 
that the DSMC technique gives good agreement with experimental data taken in low- 
Reynolds-number plumes. 

In a similar way to the previous study on N z , the DSMC computation begins at the 
nozzle throat and extends out to the near field of the plume. Initially, throat profiles 
obtained from a Navier-Stokes solution of the flow were employed to begin the DSMC 
simulation. However, it was subsequently found that simple isentropic theory could be 
used without affecting the final DSMC solution far downstream of the throat. The 
interaction of the gas with the nozzle wall is assumed to be diffuse with full 
accommodation to a temperature of 300 K. In the study by Boyd et al. (1992) it was 
found that a temperature gradient existed along the nozzle wall. Such data were not 
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Figure 8. Temperature dependence of the viscosity coefficient for molecular hydrogen. 

available for the present study and the assumption of a constant wall temperature may 
introduce a small error in the numerical results. For the nozzle simulation, the 
expansion of the plume is assumed to occur in a perfect vacuum. The computational 
grid employed in the present studies consists of 781 cells in the axial direction and 60 
cells in the radial direction. The cell dimensions are varied with local flow conditions 
to maintain a maximum dimension of one local mean free path everywhere in the 
domain. In the important regions next to the nozzle wall and at the nozzle lip the cell 
sizes employed are much smaller than the local mean free path to obtain better 
resolution. An example of the type of grid employed in the present study is shown in 
Boyd et al. (1994). At the steady stage of the calculation there are about 700000 
particles in a typical simulation. For this configuration, a typical run time is about 
1.5 CPU hours on a Cray C-90. 

3.2. Hydrogen kinetics models 

Model parameters for translational and rotational relaxation must be chosen to apply 
the DSMC technique to low-temperature flows of molecular hydrogen. The stagnation 
conditions investigated indicate that it is only necessary to consider temperatures close 
to 300 K. Relaxation of the translational mode occurs in the simulation through 
collisions between particles. These collisions are computed by applying collision 
probabilities to pairs of particles in each cell in the computational domain. In the 
current study, the form of the collision probability is determined by the Variable Hard 
Sphere (VHS) model of Bird (1981). In this model, the important parameters are a 
reference collision diameter d ref and a viscosity temperature coefficient w. These 
parameters are obtained simultaneously by considering the temperature dependence of 
the coefficient of viscosity for hydrogen. In figure 8 experimental data for viscosity are 
compared to the VHS model using parameters d ref = 2.88 x 10 _1 ° m and w = 0.67. 
These are the values given by Bird (198 1 ) for hydrogen and clearly give an excellent fit 
to the measured data. 

In the DSMC method, relaxation of the rotational mode occurs by applying a 
probability of exchanging energy between the translational and rotational modes of 
particles that collide. For many diatomic species, a suitable energy-dependent form for 
this probability has been developed by Boyd (1990). The form of this probability 
reflects the observations made experimentally that the probability of rotational 
equilibration increases as temperature decreases. However, this trend appears to be 
reversed for diatomic hydrogen. In fact, the number of collisions required to equilibrate 
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Figure 9. Temperature dependence of the rotational collision number for various forms of 

molecular hydrogen. 

the rotational mode of this species increases as the temperature is decreased. This point 
is illustrated in figure 9 where the rotational collision number for three different forms 
of hydrogen are shown as a function of temperature. The experimental data are taken 
from literature cited by Lambert (1977). As discussed previously, the rotational state 
of hydrogen occurs in two different forms: para (p-H 2 ) and ortho (0-H 2 ). Under 
normal conditions (n-H 2 ), ortho-H 2 is present as 0.75 and para-H 2 is 0.25 by mole 
fraction. 

Because the behaviour of H 2 differs from other diatomics, it is necessary to develop 
a new rotational energy exchange probability for the DSMC technique. In addition to 
the experimental data shown in figure 9, the available literature (Boitnott & Warder 
1971; Lensch & Gronig 1977; Gallagher & Fenn 1974) indicates that the rotational 
relaxation time of hydrogen at temperatures below 300 K is within a factor of 2 of 
r r = 10~ 8 s at a pressure of 1 atmosphere. The definition for the rotational collision 
number is 

Zr = 7 r /7 c , (6) 

where r c is the inverse of the molecular collision rate. Evaluating r c using the VHS 
model, and using a slightly higher value for the rotational relaxation time of t t — 
1.4 x 10 -8 s, it may be shown that the rotational collision number for hydrogen is 

Z r = 10480/7^. (7) 

This analytical form is compared to the available experimental data in figure 9. The 
simple model gives good agreement with the experimental data for normal hydrogen. 
By multiplying (7) by a factor of 2, good agreement is obtained for the experimental 
data for o-H.,. These comparisons give an idea of the spread in the available 
experimental data. Of course, this rather simple approach only gives an estimate of the 
rotational relaxation rate. Future studies may attempt to include more detailed 
analysis. 

By following a similar analysis to that presented in Boyd (1990) for other diatomic 
species, it may be shown that the appropriate energy-dependent form of the probability 
of rotational energy exchange to be used in the DSMC technique based on (7) is 

, _ 1 _ 1 /-(£+! -») (E C \ 

Vr z T 10480 nc+ 1 ) 


( 8 ) 
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In (8), E c is the total collision energy, k is the Boltzmann constant, and £ is the average 
number of degrees of freedom participating in the energy exchange process (£ = 2 for 
a collision between two diatomic molecules that both exchange rotational energy in the 
collision). Because of its approximate nature, the leading constant in (8) is considered 
a parameter to be identified by comparison with the experimental data generated in this 
investigation. 

Two separate adjustments must be made to (8) before it can be applied in the DSMC 
technique. The first is due to the implementation of collision mechanics. It has been 
shown by Lumpkin, Haas & Boyd (1991) that if the rotational energies of two particles 
are changed simultaneously in one collision then the probability must be multiplied by 
a factor 

5^ + 2^ 

5 — 2w ' 

The second adjustment is due to the fact that the collision time given by the VHS 
collision model is not generally the same as that employed in converting experimentally 
measured rotational relaxation times to rotational collision numbers. In this conversion 
process, it is common (Parker 1964) to employ a collision time of the form 

t c = (10) 

where /t is the coefficient of viscosity, and p is the pressure. Accounting for the viscosity 
and collision time produced by the VHS model (Bird 1981) the following correction 
must be applied to (8): 

8/TS-6*) 208-12a, 

2 15trr(f— w)211 — 12w(i— w)' 11 ; 

In the current application, with to = 0.67, and £ = 2, (9) and (11) are evaluated as 
Fj = 2.093 and F 2 = 0.870 respectively. The probability of exchanging translational 
and rotational energy for each collision is then given by 

<j> r = F l x F z x equation (8). (12) 

3.3. Plume flow simulation 

The code employed in the nozzle study is based on a numerically efficient program 
designed to be executed on Cray supercomputers (Boyd 1991). With minor 
modifications, the same code is used in the plume study on a computer work-station. 
The rate of execution of the code on the work-station is ten times slower than on a 
Cray YMP. 

Application of the DSMC technique to model the nozzle and near-field plume 
expansion of this unignited hydrogen arcjet is discussed above. The hydrogen kinetics 
model described above is again employed in the plume study. The flow field considered 
begins at the exit plane of the nozzle and extends to distances of 10 cm axially and 6 cm 
radially. The results computed for the nozzle flow are employed as an input condition 
at the nozzle exit plane for the present study. Decoupling the nozzle and plume 
calculations saves a great deal of computational expense. The densities in the nozzle 
flow are orders of magnitude higher than those occurring in the plume. The 
computational expense of the DSMC method is approximately proportional to the 
density of the flow. Therefore, it takes much less computer resources to perform a 
plume calculation. The smaller size of the plume problem allows the simulation to be 
performed in reasonable times on the work-station. The decoupled approach for 
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simulating nozzle and plume flows using the DSMC technique has been successfully 
applied in previous investigations (Boyd et al. 1992; Penko et al. 1993). It is observed 
that only at the highest chamber pressure considered below is there a very small 
upstream influence. The decoupled approach appears to be valid at least for the 
conditions considered in this study. 

Boundary conditions are applied along the radial line directly above the nozzle lip, 
and along the maximum axial and radial coordinates of the computational domain (at 
Z — 10 cm and R = 6 cm respectively). In each case, these are specified in terms of 
pressure and temperature. It is very difficult to perform an exact simulation of the 
chamber residual pressure. As illustrated in figure 3, the pressure is not distributed 
homogeneously throughout the chamber. Also, it is not clear what the temperature of 
the residual gas will be. For the purposes of the present investigation, a constant 
pressure ( p b ) for the background gas is assumed at each of the boundaries. In addition, 
the residual gas is assumed to be at room temperature (300 K). Particles are introduced 
into the simulation across each boundary with properties corresponding to a gas at rest 
at the conditions of the residual gas defined in this manner. This technique was applied 
in a previous investigation of helium micro-thrusters (Boyd et al. 1994), and diffuse 
shock-like structures were observed in the solutions. However, no experimental data 
were available to assess the validity of those calculated results. 


4. Results 

In this section, direct comparison is provided of the experimental data and numerical 
results obtained in this study. The comparisons are divided into those for nozzle and 
those for plume flow. 

4.1. Nozzle flow 

Flow-field properties predicted by the DSMC method are first considered before 
comparison with the experimental data is made. Contours of Mach number for the 
nozzle flow are shown in figure 10. It is noted that a very large viscous boundary layer 
develops in the nozzle flow such that only about 30% of the nozzle exit plane is 
isentropic. The DSMC method is found to successfully predict the important 
phenomenon at the nozzle lip which is intersected by the sonic line. This behaviour is 
not generally predicted by continuum methods and is due in part to the need to include 
slip effects along the nozzle wall. An idea of the degree of thermal non-equilibrium in 
the flow field is obtained from figures 1 1 (a) and 1 1 ( b ) where contours of temperatures 
obtained for the translational and rotational modes are shown. It is clear that the 
rotational mode relaxes at a significantly slower rate than the translational mode. 
These results are obtained using (12). It is concluded that rotational non-equilibrium 
plays an important role in the expansion of molecular hydrogen under conditions of 
low Reynolds number. 

Comparison of numerical results and experimental data for total number density in 
the nozzle exit plane is made in figure 12(a). Three different DSMC solutions are 
shown corresponding to Case 1 using (12), Case 2 using half of (12), and Case 3 using 
twice (12) for the rotational relaxation rate of H 2 . Note that Case 1 corresponds to the 
solid line given in figure 9 and Case 2 corresponds to the dashed line. Generally, the 
computational results for density are relatively insensitive to the exact form of the 
rotational relaxation rate, and quite good agreement with experiment is obtained. In 
both experiment and simulation a small drop in density close to the nozzle axis is 
observed. This is believed to be due to the influence of the thick viscous boundary layer 
that causes the gas on the centreline to expand more rapidly. 
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Figure 12. Comparison of experimental data and DSMC results for (a) number density 
and ( b ) rotational temperature in the nozzle exit plane. 


Rotational temperature profiles in the nozzle exit plane are compared in figure \2{b). 
As might be expected, this quantity is much more sensitive to the probability employed 
for rotational equilibration. At the axis, there is a difference in rotational temperature 
of a factor of about 2 between Cases 2 and 3, corresponding to a difference in 
relaxation rate of a factor of 4. The best agreement with experiment is obtained for 
Case 2, the slowest rate considered. In agreement with the experimental data there is 
a small rise in rotational temperature at the nozzle axis. This makes physical sense 
considering the small decrease in density described previously. To display the 
significant degree of thermal non-equilibrium predicted by the DSMC method in the 
flow, the translational and rotational temperatures obtained in Case 2 are compared 
with the experimental data in figure 13. At conditions of thermal equilibrium, the 
rotational temperature could be combined with the number-density measurements to 
estimate the pressure in the nozzle exit plane using the equation of state. The DSMC 
analysis reveals that this procedure leads to errors in pressure on the axis of a factor 
of about 5. This observation provides a strong argument in favour of the need to 
perform detailed numerical simulations and to directly measure the gas kinetic 
temperature for these non-equilibrium flows. 

The final comparison of results and data is for the number density of the first 
rotational level in the nozzle exit plane. These measurements are made at a finer spatial 
resolution than the total number density and rotational temperature measurements 
considered previously. The comparison is shown in figure 14. The DSMC results are 
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Figure 13. Degree of thermal non-equilibrium computed for Case 2 in the nozzle exit plane. 



r (mm) 


Figure 14. Comparison of experimental data and DSMC results for density of the J = 1 state in 

the nozzle exit plane. 


computed from the values for total number density and rotational temperature in the 
exit plane. Note the small degree of asymmetry in the experimental data. The data on 
the left wing are in better agreement with those shown in figure 12. For this reason, it 
is concluded that the best comparison of theory and measurement occurs for Case 2. 
The good agreement obtained using this probability function for the three different 
data sets allows a value to be estimated for the rotational relaxation rate of hydrogen 
at temperatures below 300 K. When the probability used in Case 2 is converted into a 
relaxation time the value is t T p = 2.8 x 10 -s s atm. This value is higher than those 
generally quoted for the relaxation of normal hydrogen. Recalling that a constant 
rotational relaxation time was employed in the analysis for temperatures below 300 K. 
this higher value may indicate that there is an increase in t, at low temperatures. On 
the other hand, this value is consistent with the experimental measurements for ortho- 
hydrogen. It is concluded that the relaxation time obtained in the present study is at 
least within the range of experimental data reported previously in the literature. The 
value obtained in the present study is a more complete test of the relaxation model due 
to the two-dimensional, viscous, and rapid-expansion processes experienced in the 
flow. 
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Figure 15. Contours of log 10 number density in molecules m' 3 at a background pressure of 0 Pa. 


4.2. Plume flow 

First, the DSMC method is applied to model the hydrogen-plume expansion process 
into a perfect vacuum. This is similar to the in-flight operational mode of the device 
(although here the arc is unignited). The boundary conditions are therefore specified 
as p b = 0 Pa. Contours of number density are shown in figure 15. Note that the nozzle 
exit is about 4.75 mm long, and the origin is at the nozzle exit. It is clear that the flow 
is highly directed. The variation of density along the plume centreline is shown in figure 
16(u). The density decreases rapidly by a factor of about 50 at a distance of 10 cm from 
the nozzle exit. A significant decrease in the translational temperature is shown in 
figure 16(6). The translational temperature shown in figure 16(6) is a weighted average 
of the single mode parallel to, and the two modes perpendicular to, the flow direction. 
While the former mode is frozen due to lack of collisions, the latter modes continue to 
decrease due to geometric effects. Specifically, as the distance away from the source of 
a highly expanded, low-density jet increases, the molecules located on the centreline 
tend to have smaller and smaller velocity components perpendicular to the flow 
direction. This is because only those molecules with such velocity components have not 
moved far away from the centreline during the expansion process. Thus, when the 
frozen parallel velocity component is combined with the slowly decreasing per- 
pendicular components, the overall translational temperature continues to decrease as 
shown in the computations. Also included in figure 16(6) is the rotational temperature 
which is frozen. Note the degree of thermal non-equilibrium that exists at the nozzle 
exit plane which occurred during the nozzle expansion process. 

The lowest value of the residual chamber pressure that can be achieved in the 
experimental facility is p b = 43 Pa. Contours of number density in particles m -3 
obtained using the DSMC technique for this condition are shown in figure 17. The 
finite chamber pressure creates a very thin jet with a diameter of about 1 mm. The 
diagram illustrates the important features that are typical of those observed in all 
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Figure 16. ( a ) Number density and ( b ) translational and rotational temperature, along the plume 
centreline computed with the DSMC method at a background pressure of 0 Pa. 



Figure 17. Contours of number density in molecules itt 3 at a background pressure of 43 Pa. 

(Contour levels multiplied by 10' 22 .) 
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Figure 18. (a) Number density and ( b ) translational and rotational temperature, along the plume 
centreline computed with the DSMC method at a background pressure of 43 Pa. 

calculations which include the chamber pressure. Careful examination of figure 17 
reveals that the density along the jet centreline falls then rises sharply before falling 
again. This behaviour is considered in greater detail in figures 18 (a) and 18(6) which 
show the density and temperatures obtained along the centreline. The profiles indicate 
that complex fluid dynamics processes are at work. A very strong shock wave is formed 
at a distance of about 10 mm from the nozzle exit plane. Behind this shock, the gas 
expands rapidly again before a second compression is reached. This secondary shock- 
like structure is much weaker, and the gas eventually starts to expand again further 
downstream. Although not clear in figure 18(a), a further weaker compression occurs 
at about 70 mm from the nozzle exit plane. These results are consistent with the 
diamond shock patterns observed in larger supersonic nozzles. The pressures in the 
nozzle exit plane are less than the residual pressure of 43 Pa, thus the nozzle is 
overexpanded. Actual diamond structure is not observed in the calculations. This is 
presumably due to viscous dissipation at the low densities encountered. The 
compression and expansion wave behaviour is also found in the translational 
temperature profile shown in figure 18(6). The translational mode reaches a local 
maximum after the shock formation at iOmm from the nozzle exit plane. The 
subsequent expansion is followed by further compression of the temperature up to and 
beyond the temperature of the background gas. Extension of the flow field beyond 
10 cm axially would presumably allow the translational mode to equilibrate with the 
residual gas. By comparison, the rotational mode undergoes a much slower rate of 
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Figure 19. Contours of the ratio of pressure obtained at p b = 0 Pa to that obtained at p b = 43 Pa. 

relaxation. The rotational temperatures slowly rise towards the equilibrium value. This 
illustrates the degree of thermal non-equilibrium that accompanies the complex flow- 
field structure of these interaction regions. 

The results shown in figures 18(a) and 18(b) indicate that the effective test section 
of flow undisturbed by the residual gas only extends along the axis to a distance of 
about 10 mm along the plume centreline. Radially, the influence of the chamber 
pressure is much stronger owing to the reduction in nozzle exit pressure in the radial 
direction. In figure 19 contours are shown of the ratio of pressures computed with 
p„ = 0 to those obtained at p b = 43 Pa. The contour value just less than unity 
indicates where the residual gas first affects the plume expansion process. While the flow 
along the axis is unaffected to a distance of about 10 mm from the nozzle, it is clear 
that in the radial direction interference effects are felt much sooner. Indeed, it appears 
that there is some interference close to the nozzle lip in the nozzle exit plane. 
This is confirmed in figures 20(a) and 20(b) which compare the number-density and 
translational-temperature profiles computed just downstream of the nozzle exit plane 
for the vacuum and p b = 43 Pa cases. Figure 20(a) indicates that the number density 
at the nozzle exit is only slightly affected by the back pressure beyond about 4 mm. The 
translational temperature, however, is much more sensitive as indicated by figure 
20(b). The effect of the residual gas is first felt at a radial distance of 3 mm from the 
centreline. 

Consideration is now given to comparison of the numerical results with the 
measured data. To facilitate taking high-spatial-resolution experimental data, only the 
density in the J = 1 rotational state is measured. These measurements are taken along 
the plume centreline. In the present study, the corresponding numerical prediction is 
calculated using the DSMC results obtained for the total number density and the 
rotational temperature. These are combined in the usual equilibrium relation to 
determine the number density in the first rotational level. The numerical and 
experimental results are compared in figure 21. It should be noted that each 
experimental data point is obtained over a volume that is 4 mm in length axially. The 
measurements resemble the variation in total number density shown in figure 18(a). 
The agreement between the two data sets is remarkably good. Although the computed 
shock location is further from the nozzle exit, the peak value is well predicted. While 
difficult to discern in the experimental data, the secondary compression predicted 
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Figure 20. Comparison of radial profiles for (a) number density and ( b ) temperature, near the 
nozzle exit plane obtained in vacuum and at finite back pressure. 
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Figure 21. Comparison of measured and calculated data for the density of the J — 1 state along 
the plume centreline at a background pressure of 43 Pa. 


numerically is not inconsistent with the measurements. Given the complicated fluid 
mechanics and physics of these flows, and the uncertainty in the spatial distribution of 
the residual pressure, this degree of agreement is very satisfying. The most obvious 
explanation for the difference in shock location is that the back pressure is higher at 
the DSMC boundaries than at the measuring point, which is further downstream. This 
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Figure 22. (a) Number density and ( b ) translational and rotational temperature, along the plume 
centreline computed with the DSMC method at a background pressure of 56 Pa. 

difference would move the shock position back in the simulation slightly downstream 
of the measured position. It is certainly possible that, locally along the axis, the back 
pressure may be higher than the value measured experimentally, thus pushing the 
shock closer to the nozzle exit. It would be difficult to include this type of behaviour 
in the numerical analysis. 

The results obtained at higher residual pressures in the test chamber are now 
considered. The density and temperature profiles for p b = 56 Pa computed along the 
plume centreline are shown in figures 22(a) and 22(b). The results show the same trends 
as those obtained for the lower value of p b . The shock is moved closer to the nozzle exit 
plane, and is now stronger. The peak in translational temperature is reached sooner. 
Comparison with the experimental data shown in figure 23 is once again very good. As 
before, the precise shock location is predicted slightly further downstream than that 
measured. Finally, the results obtained at the highest back pressure ( p b = 67 Pa) are 
shown in figures 24(a). 24(b), and 25. Comparison to the previous figures clearly 
illustrates the progressive increase in shock strength and movement of the shock closer 
to the nozzle exit plane as the back pressure is increased. In all these simulations, the 
DSMC shock-location prediction is slightly further from the nozzle than that observed 
experimentally. This gives a trend that is consistent with the idea proposed earlier that 
the pressure set at the DSMC boundaries may be too low since the chamber back 
pressure is actually measured further downstream of the nozzle. The general agreement 
between calculation and experiment is maintained. For the experimental data shown 
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Figure 23. Comparison of measured and calculated data for the density of the J = 1 state along 
the plume centreline at a background pressure of 56 Pa. 




Figure 24. (a) Number density and ( b ) translational and rotational temperature, along the plume 
centreline computed with the DSMC method at a background pressure of 67 Pa. 


in figure 25 it appears that the measurements at the nozzle exit plane are affected by 
the back pressure. Clearly, operation of the facility at higher pressures will produce a 
flow that is not anywhere representative of the plume expanding from the thrusier 
under flight conditions. Simulation of flows at higher values of p b would require 
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Figure 25. Comparison of measured and calculated data for the density of the / = 1 state along 
the plume centreline at a background pressure of 67 Pa. 


Rankine-Hugoniot p b - 43 Pa p b = 56 Pa p b = 67 Pa 

5.74 4.47 4.79 5.19 

Table 3. Density ratio across first shock wave in jet expansion 


simultaneous computation of the nozzle and plume flows. At p b = 67 Pa, there is a 
small amount of flow of the residual gas into the nozzle. This is considered a second- 
order effect for the results presented. However, at higher pressure, this flow would 
become more significant and could lead to a region of separated flow inside the nozzle. 

It is observed in the density profiles along the jet centreline that when the chamber 
pressure is increased the shock strength also increases. It is informative to compare the 
rise in density across the first shock wave predicted by the DSMC method for each 
value of p b with the continuum Rankine-Hugoniot theory. At the point where the 
density first starts to rise in figures 18(a), 22(a), and 24(a), the local Mach number is 
about 10. The Rankine-Hugoniot theory predicts a density ratio (pjp ^ ) across a 
normal shock for a diatomic species of 5.74 at this condition. Each of the density rises 
obtained from the three DSMC simulations is listed in table 3. As the residual pressure 
is increased, the density ratio increases and approaches the theoretical value. The 
density rise in the expanding flow is lower than that predicted by theory owing to non- 
equilibrium and viscous effects. At the onset of the density rise, there is strong thermal 
non-equilibrium between the translational and rotational modes. Through the shock 
wave, low-density collisional effects give rise to non-continuum behaviour. While 
similar in nature, these flow-field structures are not Mach disks: the flow remains 
supersonic downstream of the maximum density. 

A final note concerns the computational resources required to perform the plume 
simulations reported here. The solution times required by the DSMC method for the 
calculations with finite back pressures increase with p b . For the most expensive case of 
p b = 67 Pa. a total of 15 h on a work-station is required. The simulation of expansion 
into vacuum is the least expensive and requires about 3 h. 
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5 . Conclusions 

Advanced numerical and experimental methods were applied to investigate the low- 
density supersonic flow of a hydrogen jet. For the nozzle flow, it was found that the 
flow is in strong rotational non-equilibrium. A new rotational model was developed for 
the Monte Carlo method that was employed to compute the flow field. Computation 
of rotational temperature in the nozzle exit plane was found to be sensitive to the rate 
of rotational relaxation employed in the simulation. Comparison of the computed 
results with the experimental data was quite good and permitted the rotational 
relaxation time for hydrogen to be estimated. Fine details in the flow field, such as a 
small decrease in density and corresponding small increase in rotational temperature 
on the axis in the nozzle exit plane, were observed in both simulation and experiment. 
It was concluded that the Monte Carlo technique offers accurate prediction of the 
non-equilibrium flows produced in small nozzles used for low-thrust propulsion 
applications. 

Detailed experimental and numerical results were also obtained for the interaction 
region of the hydrogen thruster plume expanding into a finite back pressure. Despite 
the low densities involved, the characteristic phenomena involving formation of shock 
and expansion waves were observed. Direct comparison of the experimental data and 
numerical results gave very good agreement along the plume centreline. There were two 
significant findings from this portion of the present study. First, for this particular 
facility, there was only a small region of the plume centreline that was unaffected by 
the back pressure. At the lowest pressures achieved, this region extended to 10 mm 
along the centreline : this was less than 3 nozzle diameters. It is likely that interference 
effects first occurred in the radial direction at a much shorter distance from the nozzle. 
The second important conclusion was that the DSMC method is a very useful 
technique for calculation of low-density expanding flows. The complex fluid mechanics 
combined with strong viscous and non-equilibrium effects makes modelling of these 
flows challenging. The direct comparisons made in this study with experimental 
measurements indicated that the DSMC technique is capable of accurately simulating 
such complicated phenomena. It therefore appears that the method is appropriate for 
assessing the important question of plume interference effects for similar experimental 
vacuum facilities. 

Support by NASA Lewis Research Center for I.D.B. (Grant NAG3-1451) with Dr 
P. F. Penko as technical monitor is gratefully acknowledged. The work of D. R.B. and 
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Abstract 

Raman scattering measurements of molecular hydrogen density and temperature were made in the 
plume of a 1 kW-class hydrogen arcjet thruster. Pulsed-laser excitation was used to improve the signal-to-noise 
ratio in this high-background environment. Quantum-limited detection was achieved through the use of gated 
photon counting and a high-power Nd:YAG laser. Radial profiles of rotational temperature and density at the 
exit plane were measured for 5 power levels. In all cases the profiles were asymmetric about the arcjet center- 
line. The rotational temperatures were compared with the translational temperatures of atomic hydrogen from 
LIF studies and found to be significantly lower. This result suggests that this flow, like the cold flow, is not in 
translational-rotational equilibrium. 


Introduction 

Arcjet thrusters are a type of electric-powered 
rocket which can be used for propulsion of space vehi- 
cles. They offer higher exhaust velocities than conven- 
tional chemical rockets; hence they require less 
propellant mass to deliver a given impulse. For certain 
missions, arcjets offer a good compromise between the 
often conflicting requirements of small size and low mass 
for the thruster and propellant combined. They are now 
being commercially employed in station-keeping roles on 
communication satellites. 

There is a desire to increase both the thermal 
efficiency and exhaust velocities of current designs in 
order improve cost effectiveness and allow for use in a 
wider range of missions. Furthermore, there is interest in 
developing higher thrust versions which could be used in 
applications such as orbit transfer. Such improvements 
will depend partly on a better understanding of the 
plasma and gas dynamic processes undergone in the arc- 
jet nozzle. 

Currently both numerical analysis and experi- 
mental studies are being used to investigate flow proper- 
ties in low-power arcjet thrusters. Attempts at modelling 
the plasma flow in hydrogen arcjet thrusters use both 
continuum (Butler et al. 1994, Miller and Martinez- 


Sanchez 1993) and Monte Carlo techniques (Boyd et al. 
1994). These models are useful in understanding the var- 
ious loss mechanisms in an arcjet, including anode radia- 
tion, ionization, dissociation and other frozen-flow losses. 
Numerical models with accurate predictive capabilities 
will also be very helpful by allowing parametric studies to 
be performed much more easily than by testing proto- 
types in the laboratory. However, experimental measure- 
ments of relevant physical parameters in both the nozzle 
and plume are needed to validate these models. Of partic- 
ular interest are density, temperature and velocity mea- 
surements of atomic hydrogen (H) and molecular 
hydrogen (H 2 ). 

Recent studies have relied primarily on optical 
diagnostics such as emission, absorption, and fluores- 
cence. Optical techniques have the potential advantages 
of high spatial resolution and minimal perturbation to the 
flow field and most efforts have probed atomic hydrogen 
because of the large number of transitions in the visible 
spectrum. These experiments have included measure- 
ments of: H velocity and translational temperature using 
laser induced fluorescence (LIF) of the H a Balmer line 
(Liebeskind et al. 1993), H density using vacuum ultravi- 
olet absorption (Manzella and Cappelli 1992) and multi- 
photon LIF (Pobst et al. 1995), and H temperature and 
electron density using emission (Storm and Cappelli 
1994). 
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Molecular hydrogen (like other homonuclear, 
diatomic molecules) does not undergo radiative transi- 
tions between rotational or vibrational levels in the 
ground electronic state since it does not have a dipole 
moment. As a result, one can not use emission, absorp- 
tion or fluorescence techniques to study these energy 
states. Until now, measurements of H 2 properties have 
been limited to spatially averaged absorption measure- 
ments (Pollard 1992) and density and temperature mea- 
surements in the distant plume using mass sampling 
spectroscopy (Pollard 1993). However, numerical simula- 
tions (Butler et al. 1994) have predicted that -95% of the 
mass flow at the exit plane is molecular hydrogen so it is 
important to have experimental data on H 2 exit-plane 
properties. 

The research described in this paper involved 
making measurements of molecular hydrogen properties 
in the plume of a hydrogen arcjet thruster. Spontaneous 
Raman scattering was chosen as a suitable technique for 
making spatially resolved measurements of the ground 
electronic state of H 2 . Raman scattering is a linear, inelas- 
tic, two-photon scattering process. It can be used to mea- 
sure the density of various rotational and vibrational 
levels, allowing calculation of the total H 2 density and the 
rotational and vibrational temperatures. Measurements 
of the Doppler shift of the scattered light can provide the 
velocity distribution function, which gives both the bulk 
velocity and the translational temperature. The primary 
disadvantage of Raman scattering is the very low signal 
levels involved, 'typical scattering cross sections are 10 -33 
m 2 , or about 1000 times smaller than Rayleigh (elastic) 
scattering cross sections. As a result, the data collection 
process is generally very slow and it can be difficult to 
detect the signal in a luminous background, such as an 
arcjet thruster plume. 

In this study, the population of various rotational 
levels was measured, giving the rotational temperature 
and density. Initial measurements (Beattie and Cappelli 
1992) were made in the plume of a- cold-flowing arcjet 
thruster. These flow conditions were well suited to devel- 
oping both the experimental and numerical methods. 
Exit-plane profiles showed strong two-dimensionality in 
the flow, and also a moderate degree of asymmetry about 
the axis. The measured rotational temperatures on cen- 
terline were significantly higher than expected for isen- 
tropic flow, suggesting that the flow was not in rotational- 
translational equilibrium. Comparison of the results with 
Monte Carlo numerical simulations (Boyd et al. 1994) 
confirmed this result. Good agreement was obtained 
between experimental and modeling results for density 
and rotational temperature profiles by using a slightly 
lower value of the rotational relaxation rate than the 
accepted values. Axial profiles of density were also mea- 
sured at various back pressures. Shock structure was vis- 


ible downstream of the exit plane, and the shock behavior 
with changing pressure was consistent with gas dynamic 
theory and the simulation results. 

Theory 

The Raman effect is a linear, inelastic, two-pho- 
ton scattering process. An incident photon of frequency 
v 0 scatters off a particle, causing a transition in the parti- 
cle from the initial state i to a final state /. Since the sys- 
tem energy is conserved, the scattered photon must be 
shifted in energy correspondingly: 

£ moL i + £ v, i ~ £ mol / + £ v,/ 

Af v = ~A£ mol , 

where £^ ol and £ v are the energy of the molecule and 
photon respectively. In the harmonic oscillator-rigid 
rotor approximation, the molecular energy can be writ- 
ten as: 

^mol — £ trans + £ rot + £ vib + £ elec + £ nucl • (3) 

The Raman transitions under consideration here all 
begin and end in the ground electronic state. Further- 
more, changes in the translational and nuclear energy 
components as a result of the Raman scattering process 
are small. Hence the changes in can be approxi- 
mated as changes in the rotational and vibrational states, 
allowing Eq. (1) to be written as: 

A£ v = 

For a diatomic molecule, the selection rules are 
A/ = 0, ±2 and Av = 0, ±1, where A J = Jj- J i and 
Av - Vj-v t . The case of AJ = 0 and Av = 0 is actually 
Rayleigh scattering. Scattering is designated as purely 
rotational when AJ = ±2 and Av = 0, purely vibrational 
when AJ = 0 and Av = ±1, and rotational- vibrational 
when AJ = ±2 and Av = ±1. 

Theoretical Raman spectra were calculated for 
H 2 in equilibrium at a given temperature and density. The 
Raman energy shifts were taken from the experimentally- 
measured values of Veirs and Rosenblatt (1987) for pure 
rotational, vibrational and rotational- vibrational transi- 
tions from v i = 0 and J t - = 0 through J* = 5 (6 for Q 
branch). These measurements are in excellent agree- 
ment with the quantum mechanical calculations of 
Woiniewicz (1983), differing by less than one part in 10°. 

The H 2 Raman line shape is primarily determined 
by Doppler and collisional effects (Rahn et al. 1991). At 
the low pressures in these experiments, the transitions 
will be primarily Doppler broadened. The Doppler width 
for Raman scattering is different than for emission or 


absorption, since two photons are involved. By solving 
the equations of conservation of energy and momentum, 
Weber (1973) shows that the Doppler width is: 

1 

*1/ s.K. / \ — r* a -i 

Av fwhm 


- K 2ta ^H 4<v » + v ° AVR,s " 2 f +4v *] s 


which implies: 

^ total — ^trans + ^vib + ^elec + ^rot.nucl ‘ 

Hence the partition function can be factored as follows: 

® ~ ^trans^elec^vib^rot.nucl 


( 4 ) 

where 9 is the angle between the incident and scattered 
light (0 = 180° implies forward scattering), Av R is the 
Raman frequency shift, and FWHM denotes the full- 
width at half maximum. 


Raman cross sections for transitions from v i = 0 
and J t - 0 through J t = 3 were calculated by Ford and 
Browne (1973) using quantum mechanical evaluations of 
the dynamic polarizability. Their results are in good 
agreement with previous experimental measurements 
and calculations. The cross sections were obtained at the 
excitation wavelength of interest using their 8 th order 
expansion in wavelength. Cross sections for states J t = 4 
and higher were not available, so they were estimated 
from the J 2 = 3 value taking into account the v| fre- 
quency dependence. 

In a typical Raman experiment, the incident pho- 
tons are provided by a laser beam at a rate N L and light 
scattered into a solid angle di 2 is collected from a beam 
segment of length l. The rate of Raman photons scattered 
for a transition from the initial state i to the final state /is: 

n-(da\ 

N « • "‘‘“0^1 32 V" 3 - <» 

where N L is the incident laser photon rate, n 0 is the 
molecular number density, n i /n§ is the population frac- 
tion in the initial state, and (da /dQ) t j is the differen- 
tial Raman cross section for the transition from i to/. 


The population fractions n i /n 0 were calculated 
assuming the gas was in local thermodynamic equilibrium 
(LTE) with the molecules obeying Boltzmann statistics, 
in which case: 
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where q is the molecular partition function. Neglecting 
interactions between all energy modes except rotational 
and nuclear, the total molecular wavefunction for H 2 may 
be expressed as follows: 

V l ot al — ^tmris *^vib ^ elec *^rot .nucl ^ 1 ^ 


allowing the population fraction for a particular state of 
an energy mode to be calculated independently of the 
other modes. For example, the population fraction in the 
vibrational state v is: 


n ( v ) _ 9 yjb ( v ) 6 


-e^(v)/kT 


Vjr, 


2>vi b( y ) e 


-e^(v)/kT 


(10) 


which only requires knowledge of the vibrational energy 
levels and degeneracies. The vibrational populations 
were calculated using the harmonic-oscillator approxi- 
mation with As^/y = 0 — » v = 1) from Wolniewicz 
(1983) used as the vibrational energy constant. The 
error introduced by making the harmonic approximation 
is small at the temperatures involved in this study, since 
most of the molecules are in the ground vibrational state 
(>95% at T < 2000 K). The rotational energies were 
taken from a 6 parameter fit to experimental data pro- 
vided by Jennings et al. (1985). These energy levels 
agree to better than 1 part in 10 4 with the previously 
mentioned results of Veirs and Rosenblatt (1987) and 
Wolniewicz (1983). The electronic partition function was 
taken as unity. Furthermore, all measurements taken 
were integrated over essentially all translational veloci- 
ties, so the translational partition function can be 
neglected. 


Unlike the other energy modes, the rotational 
and nuclear modes interact, so their populations may not 
be calculated independently. A hydrogen nucleus (pro- 
ton) is a fermion with spin of / = 1/2. As a result, the 
wavefunction of an H 2 molecule must be antisymmetric 
with respect to interchange of the two nuclei. This prop- 
erty results in a restriction on the symmetry of the rota- 
tional wavefunction. For H 2 in the ground electronic 
state, the electronic wavefunction is symmetric. In addi- 
tion, the translational and vibrational wavefunctions are 
both symmetric. The only components of v'cotai which can 
be antisymmetric are lynuci and i// rot . There are 4 possible 
nuclear spin states of the 2 nuclei: 3 symmetric states 
with parallel spins (often referred to as ortho-hydrogen), 
and 1 antisymmetric state with opposed spins (para- 
hydrogen). The rotational wavefunctions are symmetric 
for even J and antisymmetric for odd J. Hence there are 
only two ways of satisfying the antisymmetric constraint 
on the total wavefunction: (i) coupling a symmetric 
nuclear spin state (ortho) with an antisymmetric rota- 
tional state (odd J) or (ii) coupling the antisymmetric 
nuclear spin state (para) with a symmetric rotational 
state (even J). 


The rotational-nuclear partition function is cal- 
culated as follows: 

~ -e m /kT 

^rot.nucl = X ^nucl^rot e 
J = 0 

where g TOl is the rotational degeneracy (2 J + 1), p nucl is 
the nuclear degeneracy (3 for odd J, 1 for even J), and 
the nuclear state energy £n UC i is neglected. The nuclear 
degeneracy ratio of 3:1 (ortho-para) results in a 3:1 sta- 
tistical weighting of the odd-even J states. At tempera- 
tures above 300 K, the upper rotational states are 
sufficiently populated in equilibrium that the ortho-para 
ratio is in fact 3:1. At lower temperatures, the population 
increasingly fills the J = 0 state, and the equilibrium 
ortho-para ratio decreases. In the limit as T approaches 
0, the equilibrium ortho-para ratio also approaches 0. 

Using the previously outlined techniques, the 
predicted signal from H 2 at reference conditions (T - 
300 K, P = 101 325 Pa) is shown in Figure 1 for various 
rotational and/or vibrational transitions from v i = 0. The 
plot shows the intensity versus wavelength for excitation 
by a 2 W laser at 532 ran, with l = 2 mm, di 2 = 0.024 sr, 
and a detection system efficiency (detected photons per 
collected photon) of q sys = 0.005. The Rayleigh line is 
shown for comparison, though it is -1000 times more 
intense than the strongest Raman transition. The lowest 
4 rotational states are shown, since they are the only ones 
with a significant population at this temperature. The 
pure rotational transitions, near the Rayleigh line, are the 
most intense at up to 6000 photon s -1 . However the 
slightly weaker pure vibrational transitions near 680 nm 
were used in these experiments for two reasons. First, 
they are well separated from the laser wavelength, which 
makes it easier to separate them from elastically scat- 
tered light due to both Rayleigh scattering and reflections 
off the arcjet and chamber walls. Second, since the vari- 
ous transitions are close together they can be measured 
in a single, short spectral scan. These transitions are 
shown in more detail in Figure 2. To indicate how these 
features might appear with typical instrument resolution, 
the spectrum was also convolved with a triangular func- 
tion of 0.2 nm FWHM. Under these conditions, all the 
peaks are resolved except J = 0 and J = 1 . However, only 
the wings of these two lines overlap and the centerline 
intensity is unaffected by the adjacent transition. Table 1 
shows a summary of these calculations, including the 
transition wavelengths, cross sections and intensities. 
The Q-branch transitions are highlighted and Rayleigh 
scattering information is included for comparison. 

In order to determine if Raman scattering is a 
feasible technique for use in the arcjet plume, Raman 
spectra were calculated for the estimated flow condi- 
tions. The average exit-plane density was estimated at 
lxlO 22 nr 5 using 6000 m s' 1 as a typical mass-averaged 
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Figure 1 Raman intensity versus wavelength for various 
rotational and/or vibrational transitions from H 2 at refer- 
ence conditions. 
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Figure 2 Raman intensity versus wavelength for pure 
vibrational transitions under the same conditions as 
Figure 1. 


exhaust velocity. The rotational temperature was esti- 
mated to be 1500 K based on the absorption measure- 
ments of Pollard (1993). Using these conditions, the 
spectrum shown in Figure 3 was calculated using the 
same setup parameters as for the reference calculation. 
The maximum estimated signal level is 0.6 photon s 1 
from the J - 3 transition. At these higher temperatures, 
the population is spread out over more rotational levels, 
with significant population in all 7 of the rotational states 
shown. The coupling of the rotational degeneracy to the 
nuclear spin degeneracy is evident from the alternating 
high-low intensity of the odd-even transitions. These very 
low signal estimates suggest that detection of the Raman 
scattering will be difficult, especially in the highly lumi- 
nous background of the plume. Long collection periods 
will likely be required in order to accumulate statistically 
significant data. 
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Table 1 Summary of Raman transition data at refer- 
ence conditions using an incident wavelength of 
532.1 nm, with the Q-branch transitions highlighted. 
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Experiment 

The arcjet thruster uses an electric arc to heat 
propellant to high temperatures, often 10 000 to 20 000 
K. The hot gases are then expanded through a Laval noz- 
zle. The arc passes from the cathode to the nozzle walls, 
which act as the anode. The arcjet nozzle is machined 
from a tungsten alloy containing 98% W and 2% Th. Both 
the converging and diverging sections are conical, with 
half angles of 30° and 20° respectively. The nozzle diverg- 
ing section has a 225:1 area ratio relative to the 0.64 mm 
diameter throat. The cathode is machined from a rod of 
the same alloy, with a tip angle matching that of the con- 
verging section. The cathode is installed by placing it in 
contact with the converging section, then pulling it 
0.7 mm back. A drawing of the nozzle is shown in 
Figure 4 . The anode is held by a stainless-steel anode 
housing which is bolted to the anode body. The cathode is 
insulated from the housing with a boron nitride sleeve. 
The hydrogen is injected tangentially through two small 
passages in the injector disk upstream of the cathode tip. 
The swirl imparted by this injection mechanism is 
believed to help stabilize the arc, primarily during star- 
tup. 



Figure 3 Predicted Raman intensity versus wavelength 
at the exit plane of the arc-ignited thruster. 

The arcjet used for these experiments was a lkW 
class NASA laboratory hydrogen thruster. The thruster 
was operated in a vacuum chamber 1.09 m long and 
0.53 m in diameter with optical access through 75 mm 
diameter ports.lt was mounted on a two-axis axial and 
radial translation stage. Twelve optical ports were located 
at 30“ intervals around the circumference of the tank, 
with additional ports in the endwall. The tank was evacu- 
ated by two Roots blowers backed by mechanical rough- 
ing pumps with a total capacity of 1.2 m 3 s' 1 . This system 
was capable of maintaining a back pressure of 43 ±2 Pa 
while flowing 13.3 mg s' 1 hydrogen. 

In the previous cold flow experiments (Beattie 
and Cappelli 1992, Boyd et aL 1994, Beattie 1995) a CW 
laser was used. However, in order to improve the signal- 
to-noise ratio (SNR), a pulsed laser excitation scheme 
was chosen here. The Nd:YAG laser used was a Spectra 
Physics DCR-1 flashlamp-pumped, solid-state model. The 
output was frequency doubled using a type II KDP crys- 
tal, and the 532 nm light was separated using two dich- 
roic filters. When operating at a 10 Hz repetition rate, the 
laser was capable of producing 200 mJ pulses at 532 nm 
with a 10 ns length. 

Initially, the laser beam was sent axially into the 
arcjet through a port in the chamber endwall, as in the 
cold-flow experiments. However, the beam intensity at 
nozzle was so high that each laser pulse created a tung- 
sten plasma and after a relatively short period of time, the 
nozzle was visibly eroded by the laser. In order to allow 
full laser energy to be utilized, the setup was modified to 
send the laser beam radially into the arcjet plume. The 
laser was sent in through a port 30° from vertical and 
directly out through the opposite port (Figure 5). 
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Figure 4 Drawing of the 1 kW-class arcjet nozzle and 
cathode, showing key dimensions in mm. 

The scattered light was collected at 9= 120° 
through a port located in the same circumferential plane 
as the laser at optical table level, and was focused onto a 
spectrometer entrance slit with an achromatic two-lens 
system as shown in Figure 6. The collection lens 
(f = +400 mm, d = 76 mm) was limited to // 5.7 by the 
port geometry. The focusing lens (f = +600 mm, 
d = 76 mm) was chosen to match the spectrometer /- 
number. The spatial resolution was 2 mm in the radial 
direction, determined by the slit height, and 0.25 mm in 
the axial direction, determined by the slit width. The 
spectrometer was a 1.0 m focal length //8.7 single-pass 
Czemy-Tumer design (Jarrell Ash 78-466). The grating 
was an 1800 line/mm holographic model with high effi- 
ciency from 400 to 1000 nm. A long-pass filter (Schott OG 
590) was placed before the entrance slit to further 
reduce the intensity of collected laser light to negligible 
levels (10 8 ). The spectrometer slits were 400 p m giving 
a spectral resolution of -0.2 nm. 

Light was detected at the exit slits with a 
Hamamatsu R928 photomultiplier tube (PMT) which had 
an efficiency of 6% at 680 nm and was operated at 1000 V. 
The output was amplified by an integral socket/preampli- 
fier assembly (Hamamatsu C716-01) and sent to a photon 
counter (SRS 400). The counter was triggered by the 
laser using a photodiode. The 15 ns gate width was as 
narrow as possible to maximize the background rejection 
while still being wide enough to capture the peak of the 
signal pulse including jitter. The gate delay was first set 
approximately by simultaneously monitoring the counter 
input signal and the gate output on a two-channel oscillo- 
scope. To set it more accurately, the count rate was 
recorded for a constant level of pulsed-light input to the 
detector. The delay was then scanned over the time- 
frame of interest and set at the position of the maximum 
count rate. 


Figure 5 Schematic diagram of laser and collection 
optics (front view). 



Figure 6 Schematic diagram of experimental apparatus 
(plan view) 

The counter was gated to the laser pulse using a 
photodiode. Counts are only accumulated when the 
counter is gated. The photon pulse height distribution 
peaked at -300 mV and the lower level discriminator was 
set to 100 mV. The spectrometer was scanned over the 
spectral region of interest and the photon counter output 
was recorded on a 486 computer using SRS 440 software. 

The system was calibrated for absolute efficiency 
by measuring the Raman scattering from stagnant H 2 at a 
known temperature and pressure. In order to account for 
the different efficiency at each transition, the relative 
efficiency versus wavelength was measured using a tung- 
sten strip lamp. This second calibration was necessary 
since the transitions from J =4 through 6 were too weak 
at room temperature to be accurately measured. 

Results and Analysis 
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Given the low signal levels expected in the arc- 
ignited flow', it was feasible to make measurements at 
only a few spatial positions. Since arcjet performance is 
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determined entirely by the exit plane properties, data 
was collected as close to the exit as possible. Measure- 
ments were taken at 5 radial locations each 2 mm apart. 
With spatial resolution of 2 mm radially, this mesh cov- 
ered the entire exit plane with no overlapping of seg- 
ments. The arcjet was positioned at z = +4.0 mm and the 
data acquisition program was started. The signal was 
recorded for 4000 laser pulses (-400 s) at each of the 5 
positions across the exit plane. 

Radial profiles in the arcjet plume were taken at 
5 power levels: 800, 900, 1000, 1200, and 1400 W. With 
the hydrogen mass flow at 13.3 mg s' 1 , the corresponding 
specific power levels were from 60 MJ kg -1 to 
106 MJ kg -1 . A plot of the Raman signal versus radial 
position for the 900 W case is shown in Figure 7. The sig- 
nal levels ranged from under 300 photons 
(0.7 photon s" 1 ) from state J = 1 to 5 photons 
(0.01 photon s -1 ) from J = 6. One notable feature is that 
the data is not symmetrical about z = 0. This asymmetry 
was observed at all power levels, although to varying 
degrees. In order to determine if this effect was system- 
atic, one data set was taken by scanning the arcjet up for 
some J values and down for the others. However, there 
was no correlation between the scan direction and the 
asymmetry. Furthermore, the plume emission provided 
confirmation that the arcjet was properly centered. The 
plume emission that was measured simultaneously with 
the Raman scattering using the CW gate is shown in 
Figure 8 for the 900 W case. Gaussian curves were fitted 
to the emission to check the centering and in all cases the 
center of the fit was within 50 pm of z - 0. The emission 
intensity varied with / depending on whether the Raman 
transition wavelength coincided with any H 2 emission 
lines. 
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Figure 7 Raman signal versus radial position at 900 W 
power. 
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The first step in the analysis was to correct the 
data for possible saturation effects, which arise when 
consecutive photons overlap within the temporal resolu- 
tion of the photon counting electronics (Beattie 1995). 
This correction was negligible except in the calibration 
scans, where the J = 1 peak was reduced -20% by satu- 
ration. Next, the contribution of the CW signal to the 
gated counts was subtracted, using the measured CW 
rejection ratio of 6.4x1 0 6 . This amount was no more than 
4 counts at the 1400 W power level and less than 1 count 
at 800 W. With signal levels typically 40 to 200 counts, the 
largest source of error was the statistical variation in N R , 
with the SNR equal to JN^ . The raw data was normalized 
by the laser power and relative efficiency at each transi- 
tion. Then a best fit was performed to the corrected data 
using a function which assumed that the flow was in local 
thermodynamic equilibrium (LTE) with rotational popu- 
lation distributions dictated by a rotational temperature. 
The two parameters in this fit (LTE fit) were the Ho tem- 


Figure 8 Plume emission signal measured simulta- 
neously with the Raman signal shown in Figure 7. 

perature and density. For the data presented here, states 
J = 0 through 3 were used for calibration and states J ~ 0 
through 6 were used in the plume. 

The corrected calibration data and LTE fit inten- 
sities are plotted as a function of wavelength in Figure 9 
and in a Boltzmann plot in Figure 10. The data in Figure 
9 is normalized by the laser pulse energy in Joules (-J). 
The excellent agreement is primarily due to the good 
SNR afforded by the relatively high density at 100 Pa. 
The LTE fit temperature of 303 ± 9 K is consistent with 
the measured room temperature of 300 K. The intensity 
of the transitions from J -4 through J = 6 are shown for 
the LTE fit, but were not measured because of their rela- 
tive weakness at this temperature. 


► 





Figure 9 Normalized Raman signal versus wavelength 
from calibration data. 



X (nm) 

Figure 11 Normalized Raman signal versus wavelength 
measured at x = 2 mm and 2 = 0 mm while operating at 
P = 900 W. 



Figure 10 Boltzmann plot of rotational population frac- 
tions from Figure 9. 



Figure 12 Boltzmann plot of data in Figure 11, showing 
the measured density and rotational temperature, along 
with the calculated pressure. 


Figure 11 shows a plot of the corrected data 
from Figure 7 at z = 0 mm along with the results of the 
LTE fit. As an indication of how weak the scattering pro- 
cess is, the normalized peak signal of 0.3 photon J -1 cor- 
responds to 1 detected photon for every 10 19 laser 
photons. The corresponding Boltzmann plot of the data is 
shown in Figure 12. This data was typical in that all but 2 
or 3 points fell within one standard error of the intensity 
fit. For data following a Poisson distribution, one would 
expect measurements to lie within one standard error 
68% of the time (for large AQ. Thus the observed agree- 
ment of the data with the fit is quite reasonable. 

The temperature profile across the arcjet exit 
plane is plotted in Figure 13 for the 900 W case. There is 
a relatively small variation in temperature across the exit 


plane, ranging from 750 to 1120 K. The asymmetry that 
was apparent in the raw data is also evident here, with 
the peak temperature located at z = -2 mm, rather than 
along the centerline as expected. The asymmetry 
appears again in the density profile shown in Figure 14, 
except that the peak density is located at z = +2 mm. The 
shape is more peaked than for temperature, with the 
maximum density slightly more than double the mini- 
mum. 

The temperature profiles measured at the 5 
power levels are plotted in Figure 15. Similar asymme- 
tries are evident to varying degrees in all cases. The max- 
imum temperature always occurs at either 2 = - 2 or 
2 = —4 mm. The temperatures generally rise across the 
entire exit plane as the power is increased. This behavior 
is illustrated in Figure 16, where the maximum measured 
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Figure 13 Rotational temperature versus radial position 
for the 900 W power level. 
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Figure 14 H 2 density versus radial position for the 
900 W power level. 

temperature is plotted as a function on arcjet power. The 
maximum temperature was chosen instead of the center- 
line temperature in order to allow a consistent quantity to 
be compared between profiles of varying shape. The 
results range from T m ^ x = 960 K at 800 W to 1850 K at 
1400 W, with a relatively smooth upward trend associated 
with increasing power. 

The corresponding density profiles are plotted in 
Figure 17. Again, the profiles were all asymmetric, with 
the maximum density at z - +2 mm except for the 800 W 
power level. The densities range from 2.8xl0 21 to 
8.8xl0 21 m -3 , all below the predicted value of 9xl0 21 m -3 . 
This difference is not surprising, since that prediction 
was based on the assumption that the exhaust was 100% 
Ho at a uniform velocity. 
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Figure 15 Rotational temperature versus radial posi- 
tions for 5 power levels. 
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Figure 16 Maximum rotational temperature versus arc- 
jet power level. 

Discussion 

The H 2 rotational temperatures were compared 
to previous experimental results. The H translational 
temperature was previously measured on the same facil- 
ity using LIF and emission of the Baimer a line. In the LIF 
study (Liebeskind et al. 1993), the mass flow rate was 
reported as 13.7 mg s _I (versus 13.3 mg s -1 here) and the 
power levels were 870, 1000, and 1500 W. This LIF data 
was re-analyzed to account for Stark broadening (Storm 
and Cappeili 1995) and the resulting temperatures are 
significantly lower than the originally published data. The 
emission spectra were collected with a flow of 13.3 mg s _1 
and P = 1500 W, and were Abel inverted to give spatially 
resolved data. Figure 18 shows a comparison of the exit 
plane temperatures from Raman (1400 W), LIF 
(1500 W), and emission. The re-analyzed LIF results 
show the most peaked temperature distribution, reach- 
ing a maximum of 3400 K on axis, and in agreement with 
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Figure 17 Density versus radial positions for 5 power 
levels. 


number of collisions an H 2 molecule would undergo while 
moving from the nozzle wall to the center is proportional 
to Kn~ l . The mean-free path was estimated for H 2 in H 2 
(since the flow is predominantly using the tabulated 
value at room temperature and scaling as n~ x . At the 
throat A = 0.6 pm, so Kn~ x = 500. Since Z rot = 60 at 
2000 K, then Kn~ l » Z nt and the flow is likely near rota- 
tional equilibrium. However, at the exit plane A = 500 pm 
and Kn~ 1 = 10, so KrT 1 « Z rot . Hence near the exit plane 
the density is sufficiently reduced that H 2 molecules 
could be transported to the center without undergoing 
enough collisions to maintain rotational-translational 
equilibrium. Molecules colliding with the wall are 
expected to equilibrate at the wall temperature which 
has been estimated as -1700 K (Butler et al. 1994). 
Hence the walls provide a source of H 2 with a rotational 
temperature comparable to that measured in the flow. 
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Figure 18 Comparison of temperature versus radial 
position from Raman scattering, LIF, and emission for 
P - 1500 W. 

the more recent emission measurements of temperature. 
The Raman temperatures are the lowest of all, except 
near the edges, where all methods give similar values. 
However, the Raman temperatures are consistent with 
the absorption measurements of Pollard (1993) which 
gave T toi = 1500 K averaged across the exit plane. 


Summary 

Raman scattering measurements of H 2 density 
and temperature were made at the exit plane of a 1 kW- 
class hydrogen arcjet thruster. Unlike in the cold-flow 
studies reported previously, pulsed-laser excitation was 
used in this study to improve the SNR in this high-back- 
ground environment. Quantum-limited detection was 
achieved through the use of gated photon counting and a 
high-power Nd:YAG laser. Radial profiles of rotational 
temperature and density at the exit plane were measured 
for 5 power levels. In all cases the profiles were asymmet- 
ric about the arcjet centerline. The temperature profiles 
were compared with the translational temperatures of H 
from LIF studies and the bulk temperature from contin- 
uum simulations. The rotational temperatures were sig- 
nificantly lower than the other results, implying the flow 
is not in translational-rotational equilibrium. This finding 
is not so surprising, since results from the cold-flow stud- 
ies and simple calculations indicate that rotational equi- 
librium is unlikely under cold-flow conditions. These 
results also suggest that the predictive capabilities of 
numerical simulations of such flows would be enhanced 
by incorporating a rotational energy transfer mechanism. 


These temperature measurements suggest the 
arc-ignited flow is not in rotational-translational equilib- 
rium. This result is not entirely unexpected, since the 
^ cold-flow also showed significant rotational non-equilib- 
rium. At the higher temperatures in the arc -ignited case, 
the rotational collisional number Z TOt decreases which 
. should bring the flow closer to equilibrium. However, the 
density is roughly an order of magnitude lower, which will 
have a counteracting effect by reducing the number of 
collisions. The lower rotational temperatures can poten- 
tially be explained by considering whether H 2 should be 
in rotational equilibrium as it diffuses to the centre. The 
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